REPORT  DOCUMENTATION  PAGE 


2 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  coilection  of  information  is  estimated  to  average  1  hour  per  response,  inciuding  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  compieting  and  reviewing  the  coiiection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  coilection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Ariington  VA,  22202-4302. 
Respondents  shouid  be  aware  that  notwithstanding  any  other  provision  of  iaw,  no  person  shaii  be  subject  to  any  oenaity  for  failing  to  comply  with  a  coiiection 
of  information  if  it  does  not  dispiay  a  currentiy  vaiid  OMB  controi  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 
Ph.D.  Dissertation 


1.  REPORT  DATE  (DD-MM-YYYY) 
31-07-2014 


4.  TITLE  AND  SUBTITLE 

Application  of  optical  forces  in  microhptonic  systems 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 
W911NF-09- 1-0473 


5b.  GRANT  NUMBER 


6.  AUTHORS 
ling  Ma 


5c.  PROGRAM  ELEMENT  NUMBER 
611102 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

University  of  Southern  California 
Contracts  and  Grants 
3720  S.  Flower  Street 

Los  Angeles,  CA  90089  -0701 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS 
(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

56801-MS-PCS.12 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contmed  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

Optical  forces  represent  an  exciting  new  approach  for  manipulating 

microphotonic  devices.  In  this  thesis,  the  overall  goal  is  to  invent  and  demonstrate 

novel  microphotonic  device  functionalities  based  on  optical  forces.  There  are  two 

major  tasks.  One  is  to  explore  how  optical  forces  can  be  used  to  achieve  highly 

tunable,  on-chip  photonic  devices.  The  other  task  is  to  utilize  optical  forces  for  light  assisted. 


15.  SUBJECT  TERMS 

optical  forces,  microphotonics,  optomechanics,  optical  trapping,  light-assisted  self  assembly 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

UU 

UU 

UU 

UU 

Michelle  Povinelli 


19b.  TELEPHONE  NUMBER 
213-740-8682 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 18 


Report  Title 

Application  of  optical  forces  in  microhptonic  systems 

ABSTRACT 

Optieal  forees  represent  an  exeiting  new  approaeh  for  manipulating 

mierophotonie  devices.  In  this  thesis,  the  overall  goal  is  to  invent  and  demonstrate 

novel  mierophotonie  deviee  funetionalities  based  on  optical  forces.  There  are  two 

major  tasks.  One  is  to  explore  how  optieal  forees  ean  be  used  to  aehieve  highly 

tunable,  on-ehip  photonie  deviees.  The  other  task  is  to  utilize  optieal  forees  for  light  assisted, 

template  self-assembly  of  nanoparticles. 

Optieal  forees  are  numerieally  investigated  in  different  eonfigurations. 

Attraetive  forees  exist  between  a  suspended  one-dimensional  periodie  photonie 
erystal  waveguide  and  underlying  substrate  in  a  silieon-on-insulator  platform.  It  is 
shown  that  the  optieal  foree  can  be  enhaneed  by  designing  the  waveguide  eross 
seetion  to  make  the  mode  approaeh  the  band  edge  or  substrate  light  line.  For  periodie 
waveguides,  the  optieal  foree  is  non-monotonie  with  waveguide-substrate  separation. 

This  effeet  may  enable  the  design  of  eompaet,  integrated  optieal  power  limiters. 

An  analytieal  method  is  proposed  to  calculate  optieal  forees  between  silieon 
waveguides  based  on  the  perturbation  of  effeetive  index  at  fixed  frequeney.  The 
method  is  used  to  investigate  the  meehanieal  Kerr  effeet  in  a  eoupled-waveguide 
system  with  bipolar  forees.  It  is  shown  that  positive  meehanieal  Kerr  eoeffieient 
results  from  either  an  attractive  or  repulsive  foree.  An  enhaneed  meehanieal  Kerr 
eoeffieient  several  orders  of  magnitude  larger  than  the  intrinsie  Kerr  eoeffieient  is 
obtained  in  waveguides  for  whieh  the  optieal  mode  approaehes  the  air  light  line, 
given  appropriate  design  of  the  waveguide  dimensions. 

Optieal  forees  are  proposed  to  tune  phase  and  group  birefringenee  in  parallel 
silieon  strip  waveguides.  Widely  tunable  phase  and  group  birefringenee  ean  be 
achieved  by  varying  the  pump  power,  with  maximum  values  of  0.026  and  0.13, 
respeetively.  The  giant  phase  birefringenee  allows  linear  to  circular  polarization 
conversion  within  30  ?m  for  a  pump  power  of  67  mW.  The  group  birefringenee 
gives  a  tunable  differential  group  delay  of  6  fs  between  orthogonal  polarizations. 

A  novel  photonie  erystal  lattiee  is  proposed  for  assembling  a  two-dimensional 
array  of  partieles  by  optieal  forees  with  low  power.  The  lattiee  is  ereated  by 
introdueing  a  reetangular  slot  in  eaeh  unit  eell  of  the  Suzuki-Phase  lattiee  to  enhanee 
the  light  eonfinement  of  guided  resonance  modes.  Large  quality  factors  on  the  order 
of  10^5  are  predieted  in  the  lattice.  A  signifieant  decrease  of  the  optieal  power 
required  for  optieal  trapping  ean  be  aehieved  eompared  to  the  previous  design  based 
on  square  lattiee.  Experiments  are  earried  out  to  optieally  eharaeterize  the  high-Q 
guided  resonanee  modes  with  slot  eonfinement.  The  evolution  of  the  measured 
wavelengths  and  quality  faetors  follows  the  trend  predieted  by  the  simulations. 

The  intrinsie,  radiation  loss  of  a  eoupled  resonator  optieal  waveguide  (CROW)  is 
studied  by  the  tight  binding  approximation  (TBA).  The  TBA  predicts  that  the  quality 
faetor  of  the  CROW  inereases  with  that  of  the  isolated  eavity.  The  results  provide  a 
method  to  design  CROWs  with  low  intrinsie  loss  aeross  the  entire  waveguide  band. 

The  method  may  faeilitate  the  design  of  large-area,  eoupled-eavity  modes  with  high 
quality  factor  that  nevertheless  eouple  to  normally-ineident  radiation  for  assembly. 
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List  of  Figures 

Figure  1-1  Optical  forces  between  parallel  waveguides,  (a)  schematic  structure; 

multicolored  arrow  indicates  the  direction  of  light  propagation,  (b) 
calculated,  normalized,  optical  force  as  a  function  of  waveguide 
separation.  Insets  show  modal  profiles,  with  red/blue  indicating 
positive/negative  values  of  the  electric  field  component  parallel  to 
the  air  gap .  3 

Figure  1-2  Optical  forces  (a)  between  waveguide  and  substrate,  and  (b)  between 
vertically-stacked  microrings.  Light  is  coupled  into  the  rings  by  a 
bus  waveguide.  For  both  (a)  and  (b),  the  structure  is  shown  on  top, 
with  a  multicolored  arrow  labeling  the  direction  of  light  propagation. 

The  modal  profiles  are  shown  on  the  bottom,  with  red/blue 
indicating  positive/negative  values  of  the  electric  field  component 
perpendicular  to  the  air  gap.  Yellow  arrows  show  the  direction  of  the 
optical  force.  The  system  in  (a)  has  an  attractive  mode  only,  whereas 
(b)  has  both  symmetric  (attractive)  and  anti-symmetric  (repulsive) 
modes .  4 

Figure  1-3  Schematic  of  photonic  crystal  cavities  showing  opto-mechanical 
coupling,  (a)  “zipper”  cavity,  (b)  nanobeam  cavity,  and  (c)  double¬ 
layer  photonic  crystal  slab  cavities.  In  (a)  and  (b),  the  lattice  constant 
near  the  center  of  the  device  is  decreased  slightly  (not  shown)  to 
form  a  microcavity .  ^ 

Figure  1-4  Guided-resonance  optomechanical  devices:  (a)  double-layer 
photonic  crystal  slab  supporting  guided  resonances,  and  (b) 
asymmetric  guided  resonant  structure  composed  of  photonic  crystal 
slab  and  silicon-silica  substrate .  7 

Figure  1-5  Optical  forces  generated  in  plasmonic  devices:  (a)  finite-thickness 
metal  waveguides  separated  by  a  dielectric  gap,  (b)  semi-infinite 
metal  plates  separated  by  a  dielectric  gap,  and  (c)  a  hybrid 
plasmonic  system  with  a  dielectric  waveguide  and  a  metallic 
substrate .  ^ 

Figure  1  -6  Optically-tunable  directional  coupler,  in  which  the  effective  index  of 


each  arm  is  tuned  via  an  adjustable  slot  width .  12 

Figure  1-7  Schematic  design  of  a  self-tuning  resonator .  14 


Figure  1-8  Schematic  of  light-assisted,  templated  self-assembly.  Incident  light 
from  below  excites  a  guided-resonance  mode  of  a  photonic-crystal 
slab,  giving  rise  to  optical  forces  on  nanoparticles  in  solution.  Under 
the  influence  of  the  forces,  the  nanoparticles  self  assemble  into 
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regular,  crystalline  patterns .  19 

Figure  1-9  Light-assisted,  templated  self-assembly  of  520  nm  diameter  particles 
above  a  photonic-crystal  slab,  (a-c)  Sequential  snapshots  taken  with 
the  light  beam  on.  (d)  Snapshot  taken  after  the  beam  is  turned 
off .  21 


Figure  2-1  (a)  Patterned,  ID-periodic  photonic  crystal  Si  waveguide  with  period 

a  is  separated  from  the  Si02  substrate  by  a  distance  c?  =100  nm.  r  is 
the  hole  radius,  (b)  Dispersion  relation  for  TM  (blue  solid  curve)  and 
TE  (red  dashed  curve)  modes  in  waveguide  with  r=0.3a,  and  h=w=a. 

The  gray  region  shows  the  substrate  light  cone.  The  point  where  the 
TE/TM  mode  crosses  the  light  line  is  indicated  by  an  orange/yellow 
circle .  37 

Figure  2-2  The  attractive  force  per  length  per  unit  power  as  a  function  of  period 
a  for  (a)  TM  and  (b)  TE  modes  at  wavelength  A,  =1550  nm.  One 
period  of  the  waveguide  (wv.)  is  a  cube  of  side  length  a  with  a 
circular  air  hole  of  radius  r=0.3a.  The  separation  d  between  the 
suspended  waveguide  and  substrate  is  100  nm . 

Figure  2-3  The  attractive  force  per  length  per  unit  power  as  a  function  of  40 

distance  d  for  (a)  TM  and  (b)  TE  modes  at  wavelength  A,  =1550  nm... 

Figure  3-1  (a)  Two  coupled  silicon  waveguides,  each  with  cross  section  w^h, 

separated  by  a  distance  d  and  resting  on  a  Si02  substrate  with 
suspended  length  L.  (b)  Dispersion  relation  for  modes  of  the 
suspended  section.  The  two  lowest  modes  with  electric  field  vectors 
primarily  in  the  z-direction  are  plotted  for  separations  c?=100  nm  and 
c/=1000  nm.  The  cross  section  of  each  waveguide  is  350x350  nm^. 

The  gray  region  shows  the  light  cone  for  air .  47 

Figure  3-2  Schematic  diagram  showing  a  shift  in  the  dispersion  relation  due  to  a 
change  in  separation  between  waveguides  d.  The  solid  and  dashed 
lines  represent  the  initial  and  shifted  dispersion  relation, 
respectively .  48 

Figure  3-3  For  fixed  wavelength  1550  nm  and  waveguide  dimensions  h=w=350 
nm:  (a)  Symmetric  mode:  ties  as  a  function  of  d  (black  curve); 
F/{LP)  calculated  by  the  MST  method  (red  curve)  and  by  Equation 
(3-8)  (red  triangles),  (b)  Antisymmetric  mode:  tisf  as  a  function  of  d 
(black  curve);  F/{LP)  calculated  by  the  MST  method  (red  curve)  and  ^  j 
by  Equation  (3-8)  (red  triangles) . 

Figure  3-4  Symmetric  mode:  (a)  Band  structure.  The  gray  region  represents  the 
air  light  cone.  The  red  line  shows  fixed  wavelength  1550  nm 
(frequency  193.55  TFIz).  (b)  Attractive  force  per  length  per  unit 
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power  F/(LP).  The  initial  separation  d=100  nm,  suspended  length 
L=20  pm,  optical  power  T*  =  1 0  mW,  and  the  signal  wavelength  is 
1550  nm.  (c)  Mechanical  Kerr  coefficient  . 

Figure  3-5  Anti-symmetric  mode:  (a)  Band  structure.  The  gray  region 
represents  the  air  light  cone.  The  red  line  shows  fixed  wavelength 
1550  nm  (frequency  193.55  TFIz).  (b)  Repulsive  force  per  length  per 
unit  power  F/(LP).  The  initial  separation  g?=100  nm,  suspended 
length  L  =  20  pm,  optical  power  R  =  10  mW,  and  the  signal 
wavelength  is  1550  nm.  (c)  Mechanical  Kerr  coefficient  . 


Figure  4-1  (a)  Two  coupled  Si  waveguides,  each  with  cross  section  wxh 

separated  by  a  distance  d,  rest  on  a  Si02  substrate  with  a  free¬ 
standing  section  of  length  L.  (b)  Dispersion  relation  for  the  lowest- 
frequency  TE  mode  (solid  lines)  and  TM  mode  (dashed  lines)  of  the 
coupled  waveguides  for  several  separations.  Insets  respectively 
show  the  Ey  field  distribution  of  the  TE  mode  and  the  E^  field 
distribution  of  the  TM  mode  with  d=0.2a  at  frequency  ftia/27ic=0.18 
(darker  shades  correspond  to  larger  magnitudes  of  the  electric  field 
at  a  snapshot  in  time).  The  yellow  region  shows  the  light 
cone .  66 

Figure  4-2  (a)  Normalized  force  per  unit  area  for  the  lowest-frequency  TE  (red) 

and  TM  (black)  modes  as  a  function  of  pump  light  frequency,  at 
fixed  separation  d=035a.  (b)  Force  per  unit  area  as  a  function  of 
separation  (red  triangles),  at  optimized  frequency 
cOpa/2’iic=a/Xp=0A65.  The  right  and  top  axes  are  in  physical  units 
with  incident  power  P=20  mW,  w=h=a=263.5  nm,  and  L=30  pm. 

The  blue  solid  line  is  the  second-order  polynomial  fit  of  the  force 
per  unit  area . 


Figure  4-3  Displacement  of  the  suspended  section  of  each  waveguide  as  a 
function  of  position  along  the  waveguide.  The  cross-sectional 
dimension  is  iv=/i=a=263.5nm,  the  initial  waveguide  separation  is 
d=035a=92.2  nm,  the  pump  frequency  is  cOpa/2Ti  c=allp=0.\65,  the 
incident  power  is  P=20  mW,  and  the  suspended  length  is  L=30  pm. 
Note  that  the  x  andy  axes  differ  in  scale . 

Figure  4-4  (a)  Phase  birefringence  as  a  function  of  signal  frequency  in  the 

coupled  waveguides  with  varying  separations.  The  arrow  shows  that 
at  a  frequency  a)sal2'Kc=0.\l,  the  absolute  value  of  increases  as 
the  waveguide  separation  d  decreases,  (b)  Phase  birefringence  An^  as 
a  function  of  position  along  the  waveguides.  The  initial  separation  is 
d=035a=92.2  nm.  The  attractive  force  is  induced  by  CW  pump  light 
at  frequency  cOpa/2’iic=a/lp=0A65  and  power  P=20  mW.  The 
difference  between  the  maximum  and  the  minimum  of  Aup  is  0.026.. 
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Figure  4-6 


Figure  4-6 


Figure  4-7 


Figure  5-1 


Figure  5-2 


Figure  5-3 


Figure  5-4 


Tuning  the  relative  phase  shift  by  increasing  CW  pump  power  from 
20  mW  to  70  mW  (black  squares).  The  waveguide  length  is  23.05 
pm.  The  red  line  is  a  2”‘*  -order  polynomial  fit  of  the  phase  shift.  A 
power  of  approximately  67  mW  yields  a  0.57t  phase  difference .  73 

Phase  birefringence  Anp  as  a  function  of  signal  frequency  in  the 
two-waveguide  system  with  varying  separations,  (a)  The  cross 
section  of  each  waveguide  is  ax2a.  The  arrow  shows  that  Anp 
decreases  with  decreasing  separation,  (b)  The  cross  section  of  each 
waveguide  is  ax0.5a.  The  absolute  value  of  Anp  increases  with 
decreasing  separation . 

(a)  Group  birefringence  Atig  as  a  function  of  signal  frequency,  in  the 
coupled  waveguides  with  varying  separations.  The  arrow  shows  that 
as  the  waveguide  separation  d  increases,  the  group  birefringence  Aug 
increases  fastest  at  a  frequency  a)sal2nc=0.\l .  (b)  Group 

birefringence  Arig  as  a  function  of  position  along  the  waveguides. 

The  waveguides  are  deformed  by  an  attractive  force  induced  by  a 
CW  pump  at  frequency  a)pal2nc=allp=0.\65,  with  a  power  P=20 
mW.  The  difference  between  the  maximum  and  the  minimum  of  Arig 
is  0.13 . 

(a)  Diagram  of  the  Suzuki-Phase  photonic  crystal  lattice,  (b) 
Normalized  transmission  spectra  calculated  by  the  transfer-matrix 
method.  Red  line  for  x-polarization,  blue  for  y-polarization.  (c)  H^- 
field  profile  (left)  and  (right)  of  ol  resonance,  d)  profile  84 

(left)  and  (right)  of  e2  resonance . 

(a)  Diagram  of  the  Slot-Suzuki-phase  hybrid  lattice,  (b)  Normalized 
transmission  spectra.  Red  line  for  x  -polarization,  blue  for  y- 
polarization.  (c)  F/z-field  profile  of  se2  resonance,  d)  E^  field  profile 
of  se2  resonance .  86 

(a)  Evolution  of  the  Q  and  wavelength  Ir  as  a  function  of  slot  width 
Wy.  b)  Evolution  of  the  Q  and  wavelength  1r  as  a  function  of  slot 
length  w^.  Red  color  indicates  Q;  blue  color  indicates  wavelength. 

Dots  correspond  to  calculated  values;  lines  represent  a  guide  for  the  88 
eye . 

(a)  Optical  force  as  a  function  of  slot  width  Wy  for  four  different 
particle  radii.  The  slot  length  is  fixed  to  464  nm.  (b)  Optical  force 
E)  as  a  function  of  particle  radius.  Dots  represent  the  calculated 
values.  Eines  represent  a  guide  for  the  eye.  In  both  (a)  and  (b),  the 
particle  is  at  (x  =  0,  y  =  0)  and  has  its  bottom  edge  45  nm  above  the  ^  j 
top  surface  of  the  slab . 
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Figure  5-5  (a)  Diagram  of  the  Slot  Suzuki-phase  (SSP)  lattice  with  a  particle 

above,  (b)  Vertical  force  as  a  function  of  particle  position  in  the 
XY  plane,  (c)  In-plane  force  Fjcy  as  a  function  of  position  in  the  XY 
plane.  The  force  magnitude  is  indicated  by  the  colormap,  and  the 
force  direction  is  shown  by  the  blue  arrows,  (d)  Potential  map  in  the 
XY  plane.  To  obtain  the  results  shown  in  (b),  (c),  and  (d),  we  assume 
that  the  particle  has  a  radius  of  25  nm  and  is  placed  such  that  its 
bottom  edge  is  45  nm  above  the  top  surface  of  the 
slab . 


Figure  5-6  (a)  In-plane  optical  force  F^z  for  a  particle  with  radius  of  25  nm.  (b) 

Optical  potential  map  for  the  particle .  94 

Figure  6-1  (a)  Diagram  of  the  graphite  photonic  crystal  lattice,  (b)  Diagram  of 

the  Slot-graphite  photonic  crystal  lattice.  The  arrows  represent  the 
primitive  vectors  of  the  Bravais  lattice .  101 


Figure  6-2  (a)  Photonic  band  structure  of  the  Graphite  photonic  crystal  lattice. 

(b)  Photonic  band  structure  of  the  Slot-graphite  photonic  crystal 
lattice.  The  highlighted  region  represents  the  targeted  bands,  (c) 
field  profde  (left)  and  (right)  profde  of  the  fifth  band  (down)  and 
fourth  band  (up)  of  the  graphite  photonic  crystal  lattice,  (d)  F/z-field 
profile  (left)  and  (right)  profile  of  the  fourth  band  (down)  and 
fifth  band  (up)  of  the  Slot-graphite  photonic  crystal  lattice .  102 

Figure  6-3  (a):  Evolution  of  the  Q  and  wavelength  as  function  of  Wy.  (b): 

Evolution  of  the  Q  and  wavelength  as  function  of  Wx.  Red  color  for 
Q.  Black  for  wavelength.  Dots:  Corresponds  to  the  calculate  values. 


Eines  represents  a  guide  for  the  eye .  104 

Figure  6-4  SEM  image  of  one  fabricated  Slot-graphite  device .  105 


Figure  6-5  (a)  Diagram  of  experimental  setup  used  to  characterize  the  optical 

properties  of  the  guided  resonant  device.  Two  cross  polarizers  (PC) 
are  applied  before  and  after  the  device  to  cancel  out  Fabry-Perot 
noise.  TE:  tunable  laser;  MO:  microscope  objective;  ES:  lens;  TS: 
translation  stage;  PD:  photodiode,  (b)  One  measured  spectrum.  The 
right  transmission  peak  correspond  to  the  slot  mode,  and  the  left 
peak  shows  the  dipole  mode  with  lower  Q-factor .  j 

Figure  6-6  (a):  Evolution  of  the  Q  factor  as  function  of  Wy.  (b):  Evolution  of  the 

wavelength  as  function  Wy.  Black  spheres  for  experimental  data. 

Gray  dashed  curves  for  simulation.  Blue  stars  for  experimental  data 
when  Wy  is  increased  by  10%  by  considering  SEM  system  errors .  108 

Figure  7-1  (a)  Schematic  diagram  of  CROW  structure  E2S3.  (b)  Dispersion 
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curve  of  CROW  mode  in  L2S2  (red),  L2S3  (blue)  and  L2S4 
(magenta).  Dashed  line  shows  the  fundamental  mode  of  the  isolated 
L2  cavity . 

Figure  7-2  (a)  Schematic  picture  of  tight-binding  prediction  for  the  power 

spectrum  of  a  CROW  mode,  (b)  CROW  power  spectrum  giving  a 
large  value  of  Q{k)  compared  to  the  isolated  cavity,  (c)  CROW 
Power  spectrum  giving  a  small  value  of  Q{k)  compared  to  the 
isolated  cavity . 


Figure  7-3  (a)  Mode  distribution  of  the  real  part  of  Ey  on  the  x-y  plane  0.2a 

above  the  surface  of  the  isolated  L2  cavity,  (b)  Power  spectrum  of 
L2  cavity.  The  blue  circle  is  the  light  cone,  (c)  Power  spectrum  of 
CROW  mode  in  L2S3  structure  with  A:  =  0.  (d)  Visual  test  of  tight- 
binding  approximation  via  comparison  of  CROW  power  spectra  to  ^22 
that  of  isolated  L2  cavity.  The  blue  dashed  lines  indicate  light  lines... 

Figure  7-4  Q{k)  for  L2  CROWs  with  varying  cavity  separations:  (a)  L2S2,  (b) 

L2S3,  and  (c)  L2S4.  Dots  represent  calculated  values.  Lines 
represent  guides  for  the  eye .  125 


Figure  7-5  (a)  Quality  factor  for  CROWs  based  on  different  constituent 

cavities.  Triangles  indicate  Q(0);  circles  indicate  Qiksz)-  Dashed 
lines  show  Q  for  single  cavity.  Calculations  are  performed  using  the 
TEA  method,  (b)  Q{k)  for  CROWs  designed  with  a  number  of 
separation  holes  equal  to  the  first  node  in  (a) . 
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Abstract 

Optical  forces  represent  an  exeiting  new  approaeh  for  manipulating 
mierophotonic  deviees.  In  this  thesis,  the  overall  goal  is  to  invent  and  demonstrate 
novel  microphotonic  device  functionalities  based  on  optical  forces.  There  are  two 
major  tasks.  One  is  to  explore  how  optical  forces  can  be  used  to  achieve  highly 
tunable,  on-ehip  photonie  devices.  The  other  task  is  to  utilize  optical  forces  for  light- 
assisted,  template  self-assembly  of  nanopartieles. 

Optical  forces  are  numerieally  investigated  in  different  configurations. 
Attractive  forces  exist  between  a  suspended  one-dimensional  periodic  photonic 
crystal  waveguide  and  underlying  substrate  in  a  silicon-on-insulator  platform.  It  is 
shown  that  the  optical  force  ean  be  enhanced  by  designing  the  waveguide  cross 
section  to  make  the  mode  approaeh  the  band  edge  or  substrate  light  line.  For  periodic 
waveguides,  the  optical  force  is  non-monotonic  with  waveguide-substrate  separation. 
This  effect  may  enable  the  design  of  compact,  integrated  optical  power  limiters. 

An  analytieal  method  is  proposed  to  caleulate  optical  forces  between  silieon 
waveguides  based  on  the  perturbation  of  effective  index  at  fixed  frequeney.  The 
method  is  used  to  investigate  the  meehanical  Kerr  effect  in  a  coupled-waveguide 
system  with  bipolar  forces.  It  is  shown  that  positive  meehanical  Kerr  eoefficient 
results  from  either  an  attraetive  or  repulsive  foree.  An  enhanced  mechanieal  Kerr 
coefficient  several  orders  of  magnitude  larger  than  the  intrinsie  Kerr  coefficient  is 
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obtained  in  waveguides  for  whieh  the  optieal  mode  approaches  the  air  light  line, 
given  appropriate  design  of  the  waveguide  dimensions. 

Optical  forces  are  proposed  to  tune  phase  and  group  birefringence  in  parallel 
silicon  strip  waveguides.  Widely  tunable  phase  and  group  birefringence  can  be 
achieved  by  varying  the  pump  power,  with  maximum  values  of  0.026  and  0.13, 
respectively.  The  giant  phase  birefringence  allows  linear  to  circular  polarization 
conversion  within  30  pm  for  a  pump  power  of  67  mW.  The  group  birefringence 
gives  a  tunable  differential  group  delay  of  6  fs  between  orthogonal  polarizations. 

A  novel  photonic  crystal  lattice  is  proposed  for  assembling  a  two-dimensional 
array  of  particles  by  optical  forces  with  low  power.  The  lattice  is  created  by 
introducing  a  rectangular  slot  in  each  unit  cell  of  the  Suzuki-Phase  lattice  to  enhance 
the  light  confinement  of  guided  resonance  modes.  Large  quality  factors  on  the  order 
of  10^  are  predicted  in  the  lattice.  A  significant  decrease  of  the  optical  power 
required  for  optical  trapping  can  be  achieved  compared  to  the  previous  design  based 
on  square  lattice.  Experiments  are  carried  out  to  optically  characterize  the  high-Q 
guided  resonance  modes  with  slot  confinement.  The  evolution  of  the  measured 
wavelengths  and  quality  factors  follows  the  trend  predicted  by  the  simulations. 

The  intrinsic,  radiation  loss  of  a  coupled  resonator  optical  waveguide  (CROW)  is 
studied  by  the  tight  binding  approximation  (TBA).  The  TBA  predicts  that  the  quality 
factor  of  the  CROW  increases  with  that  of  the  isolated  cavity.  The  results  provide  a 
method  to  design  CROWs  with  low  intrinsic  loss  across  the  entire  waveguide  band. 
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The  method  may  faeilitate  the  design  of  large-area,  eoupled-eavity  modes  with  high 
quality  faetor  that  nevertheless  eouple  to  normally-ineident  radiation  for  assembly. 
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Chapter  1:  Introduction 

1.1  Background 

It  has  long  been  known  that  light  can  exert  a  force  on  objects,  an  effect  known  as  radiation 
pressure.  Back  in  the  16th  century,  it  was  suggested  that  mechanical  effects  from  solar  radiation 
cause  comet  tails  to  point  away  from  the  sun.  In  2010,  Japanese  Aerospace  Exploration  Agency 
successfully  launched  a  spacecraft  with  radiation  pressure  from  sunlight  as  propulsion  force. 

Under  most  everyday  circumstances,  optical  forces  are  too  weak  to  cause  noticeable  effects. 
For  nano-  and  micro-scale  objects,  however,  optical  forces  can  have  significant  effects.  One 
well-known  application  of  optical  forces  is  optical  trapping  [1,  2],  a  technique  that  uses  laser 
light  to  trap  and  position  small  particles.  The  objective  of  this  thesis  is  to  explore  the  applications 
of  optical  forces  in  nano-  and  micro-scale  region. 

Dr.  Povinelli  [3]  has  introduced  the  idea  of  using  optical  forces  to  precisely  control  the 
positions  of  microscale  optical  waveguides  and  resonators.  Optical  waveguides,  which  act  like 
“photonic  wires,”  and  microresonators,  devices  that  confine  light  to  the  scale  of  the  wavelength, 
are  important  elements  for  developing  photonic  integrated  circuits  [4-6].  Ultimately,  photonic 
integrated  circuits  may  manipulate  photons  in  the  way  that  electronic  circuits  manipulate 
electrons.  The  development  of  photonic  circuits  that  can  filter,  buffer,  and  reroute  optical  signals 
is  expected  to  transform  optical  communications  [7],  enable  high-speed  optical  interconnects 
between  multi-chip  computer  processors  [8],  and  provide  a  platform  for  quantum  communication 
[9].  However,  current  photonic  circuits  are  largely  static;  methods  for  reconfiguring  their 
response  are  limited.  Optical  forces  in  integrated  photonics,  or  optomechanical  effects,  represent 
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an  exciting  new  approach  for  achieving  flexible,  all-optical  tuning  and  reconfiguration  of  on-chip 
microphotonic  devices. 

1.2  Optical  Forces  in  Microphotonic  Systems 

In  2005,  Povinelli  et  al.  theoretically  investigated  the  optical  forces  between  parallel 
waveguides  [3],  shown  in  Figure  1-1  (a).  A  mechanical  force  arises  from  the  overlap  of  the 
guided  waves  propagating  in  the  two  waveguides.  The  sign  of  the  force  can  be  either  attractive  or 
repulsive,  depending  on  the  relative  phase  of  light  in  the  waveguides.  The  “bonding”  mode, 
which  is  symmetric  with  respect  to  the  two  waveguides,  gives  an  attractive  optical  force,  whereas 
the  “anti-bonding,”  or  anti-symmetric  mode,  is  repulsive.  By  controlling  which  mode  is  input 
into  the  waveguides,  input  light  can  be  used  to  pull  the  waveguides  together  or  push  them  apart. 
Physically,  the  force  originates  from  the  interaction  of  dipoles  induced  in  the  dielectric 
waveguides  by  the  electromagnetic  field  of  the  light  wave.  The  optical  force  scales  linearly  with 
input  power. 

The  calculated  optical  force  in  a  typical  system  is  shown  in  Figure  1-1  (b).  The  waveguides 
each  have  cross  sections  of  310  nm  by  310  nm,  and  the  wavelength  in  air  is  1550  nm.  In  the 
figure,  the  force  per  unit  length  of  the  waveguide  is  divided  by  the  input  power  and  plotted  as  a 
function  of  waveguide  separation.  For  waveguide  separations  larger  than  -60  nm,  the  symmetric 
and  antisymmetric  modes  correspond  to  attractive  (positive  sign)  and  repulsive  (negative  sign) 
forces,  respectively.  The  magnitude  of  the  optical  force  decays  monotonically  with  increasing 
separation.  The  decay  is  due  to  reduced  optical  coupling  between  the  waveguides,  that  is. 


reduced  overlap  between  the  evanescent  tails  of  the  individual  waveguide  modes.  For  smaller 
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waveguide  separations,  strong  coupling  between  the  waveguides  can  significantly  perturb  the 
individual  waveguide  modes.  In  this  regime,  the  optical  force  of  the  antisymmetric  mode 
changes  sign. 


Figure  1-1  Optical  forces  between  parallel  waveguides,  (a)  schematic  structure; 
multicolored  arrow  indicates  the  direction  of  light  propagation,  (b)  calculated, 
normalized,  optical  force  as  a  function  of  waveguide  separation.  Insets  show  modal 
profiles,  with  red/blue  indicating  positive/negative  values  of  the  electric  field 
component  parallel  to  the  air  gap. 


In  Figure  I-I  (b),  it  can  be  seen  that  for  the  symmetric  (attractive)  mode,  the  optical  force 
increases  with  decreasing  gap.  However,  if  the  gap  is  too  small,  the  Casimir  force  becomes 
comparable  to,  or  larger  than,  the  optical  force  [10].  Devices  with  gaps  as  small  as  ~80  nm  have 
been  used  in  experiments  [11]. 

Interestingly,  the  optical  force  between  waveguides  is  in  a  direction  perpendicular  to  the 
direction  of  light  propagation.  This  is  different  than  the  traditional  case  of  radiation  pressure, 
where  the  optical  force  on  a  surface  is  parallel  to  the  direction  of  light  propagation.  The  radiation 
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pressure  ean  be  simply  derived  from  momentum  eonservation  arguments:  refleetion  of  a  photon 
from  a  surface  must  correspond  to  transfer  of  mechanical  momentum  to  the  surface.  Similarly, 
the  optical  force  on  the  mirrors  of  a  Fabry-Perot  cavity  acts  perpendicular  to  the  direction  of  light 
propagation  [12].  For  parallel  waveguides,  the  situation  is  notably  different:  the  momentum 
acquired  by  each  individual  waveguide  is  perpendicular  to  the  propagation  direction,  and  the 
momenta  of  the  two  waveguides  are  equal  and  opposite. 


1,2,1  Early  Experiments 

Initial  theoretical  work  on  optical  forces  was  followed  rapidly  by  experimental 
demonstrations. 


Figure  1-2  Optical  forces  (a)  between  waveguide  and  substrate,  and  (b)  between 
vertically-stacked  microrings.  Light  is  coupled  into  the  rings  by  a  bus  waveguide.  For 
both  (a)  and  (b),  the  structure  is  shown  on  top,  with  a  multicolored  arrow  labeling  the 
direction  of  light  propagation.  The  modal  profiles  are  shown  on  the  bottom,  with 
red/blue  indicating  positive/negative  values  of  the  electric  field  component 
perpendicular  to  the  air  gap.  Yellow  arrows  show  the  direction  of  the  optical  force.  The 
system  in  (a)  has  an  attractive  mode  only,  whereas  (b)  has  both  symmetric  (attractive) 

and  anti-symmetric  (repulsive)  modes. 
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In  2008,  a  suspended  waveguide  eoupled  to  an  underlying  substrate  was  designed  [13]  and 
fabricated  [14]  by  Li  et  al.  A  schematic  is  shown  in  Figure  1-2  (a).  The  authors  demonstrated  an 
attractive  optical  force  that  is  induced  by  the  overlap  between  the  evanescent  tail  of  waveguide 
mode  and  the  underlying  oxide  substrate.  Waveguide  displacement  was  inferred  from  the  optical 
transmission  using  an  on-chip  interferometer  scheme.  Further  experiments  demonstrated 
attractive  and  repulsive  forces  in  coupled  microresonators,  shown  in  Figure  1-2  (b). 


1,2,2  Other  Configurations 

Following  initial  work  on  waveguide  and  microresonator  systems,  a  variety  of  alternate 
microphotonic  device  geometries  have  been  investigated.  Various  mechanisms  have  been  studied 
for  enhancing  the  optical  force.  These  include  tailoring  the  overlap  between  photonic  and 
mechanical  modes,  exploiting  guided  resonance  modes,  and  incorporating  plasmonic  materials. 


Figure  1-3  Schematic  of  photonic  crystal  cavities  showing  opto-mechanical  coupling,  (a) 
“zipper”  cavity,  (b)  nanobeam  cavity,  and  (c)  double-layer  photonic  crystal  slab  cavities. 
In  (a)  and  (b),  the  lattice  constant  near  the  center  of  the  device  is  decreased  slightly  (not 

shown)  to  form  a  microcavity. 
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Photonic-crystal  microcavities  can  be  designed  to  eonfine  light  in  ultra-small  mode  volumes 
and  exhibit  high  optieal  quality  faetors,  Q  [15-18].  Reeent  work  has  exploited  these  attributes  to 
aehieve  strong  optomeehanical  effeets.  Painter  and  co-workers  [19,  20]  have  designed  and 
experimentally  demonstrated  a  “zipper”  system  with  large  optomechanieal  eoupling.  The  zipper 
is  made  up  of  two  parallel,  patterned  nanobeams  (shown  in  Figure  1-3  (a))  eaeh  of  whieh  is 
designed  to  support  an  optical  microcavity  mode.  Coupling  between  the  modes  gives  rise  to 
attraetive  and  repulsive  forces.  A  relevant  figure  of  merit  in  the  strueture  is  the  optomeehanical 
coupling  length,  =  co  I  ,  where  co  is  the  optical  resonance  frequency.  Experimental 
values  lower  than  3  microns  were  obtained. 

In  this  system,  optomeehanieal  coupling  results  in  an  “optieal  spring  effect:”  due  to  the 
presenee  of  optical  forces,  the  mechanical  resonance  frequency  and  spring  constant  are  shifted. 
In  experiments,  a  modified  stiffness  of  five  times  the  intrinsie  meehanical  stiffness  of  the 
nanobeam  was  observed. 

In  Reference  [21],  Painter  and  eo- workers  designed  and  experimentally  demonstrated  a 
nanobeam  strueture  (in  Figure  1-3  (b))  that  confines  both  the  optical  modes  and  the  mechanical 
vibration  modes  to  the  seale  of  the  optical  wavelength.  Localized  meehanieal  modes  are  created 
using  a  phononic  band  gap  microcavity,  analogous  to  the  optieal  case.  The  authors  suggest  that 
the  strong  optomeehanieal  eoupling  between  optical  modes  (200  THz  optieal  resonance 
frequeney)  and  high-frequeney  meehanical  modes  (2  GHz  mechanical  resonance  frequeney)  will 
allow  extremely  sensitive  mass  detection  via  optical  readout. 

Notomi  and  coworkers  have  theoretically  analyzed  microeavities  in  parallel  photonie -crystal 
slabs  [22].  A  schematic  double-layer  cavity  is  shown  in  Figure  1-3  (e).  Intriguingly,  the  authors 
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suggest  that  if  the  meehanieal  displacement  between  the  structures  occurs  faster  than  the 
response  time  of  the  optical  cavity  (proportional  to  the  optical  period  times  the  optical  Q), 
wavelength  conversion  will  occur  [23].  This  phenomenon  is  analogous  to  the  frequency  shifts 
predicted  and  measured  in  microphotonic  systems  with  time-varying  linear  refractive  index  [24, 
25], 


Figure  1-4  Guided-resonance  optomechanical  devices:  (a)  double-layer  photonic  crystal 
slab  supporting  guided  resonances,  and  (b)  asymmetric  guided  resonant  structure 
composed  of  photonic  crystal  slab  and  silicon-silica  substrate. 

Other  than  microcavities,  optical  forces  have  been  demonstrated  in  guided  resonance  systems. 
When  light  is  normally  incident  upon  a  ID-  or  2D-periodically  patterned  slab,  it  can  excite 
guided  resonance  modes,  which  propagate  in  the  plane  of  the  slab  while  leaking  partially  to  the 
surrounding  air  [26].  Early  work  examined  forces  between  identical,  parallel  photonic -crystal 
slabs  (Figure  1-4  (a)).  Liu  et  al.  [27]  showed  theoretically  that  optical  forces  are  enhanced  near 
guided  resonances.  In  an  ideal  system,  infmite-g  resonances,  or  “dark  states”  can  be  designed, 
and  the  force  will  diverge  as  1/g  near  the  resonant  frequency.  Experiments  on  a  similar  system 
were  reported  by  Roh  et  al.  [28].  Changes  in  separation  between  the  slabs  due  to  optical  forces 
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were  inferred  from  the  optical  reflection  spectrum.  Displacements  of  approximately  3.6  nm  were 
obtained  for  -'20  mW  of  power.  The  force  per  photon,  0.051  pN/photon,  was  similar  to  values 
obtained  in  double-layer  disk  cavities  [29]  and  zipper  cavities  [20]. 

Rodriguez  et  al.  [30]  have  theoretically  studied  asymmetric  systems  consisting  of  a  photonic- 
crystal  slab  above  an  unpatterned,  layered  substrate,  shown  in  Figure  1-4  (b).  This  system  may 
offer  advantages  in  terms  of  fabrication.  The  asymmetric  system  supports  both  quasi-symmetric 
and  quasi-antisymmetric  modes  with  opposite  signs  of  the  force.  As  for  symmetric  structures,  the 
sign  of  the  force  can  be  changed  by  tuning  the  frequency  of  the  incident  light.  Assuming  the 
photonic-crystal  slab  can  move,  changing  the  separation  between  slab  and  substrate,  bistable 
behavior  is  predicted  [31].  The  precise  conditions  for  bistability  are  affected  by  the  Casimir  force, 
suggesting  that  observations  of  bistability  could  be  used  to  infer  the  magnitude  of  the  Casmir 
force.  In  addition,  the  optical  force  might  be  used  as  a  way  of  controlling  or  preventing  stiction, 
which  results  from  Casimir  attraction. 

One  advantage  of  such  guided-resonant  devices  is  that  light  can  be  coupled  into  the  mode  at 
normal  incidence  using  free-space  optics,  potentially  allowing  ease  of  alignment.  However,  the 
mass  of  guided-resonant  slabs  is  usually  bigger  than  that  of  photonic  crystal  cavity  devices, 
which  may  reduce  the  mechanical  resonance  frequency  and  hence  increase  the  characteristic 
response  time  for  an  optical  force. 

Surface  Plasmon  Polaritons  (SPP’s)  are  propagating  electromagnetic  modes  formed  at 
metallic-dielectric  interfaces  [32]  that  offer  the  ability  to  confine  light  to  the  deep  subwavelength 
scale.  There  has  been  interest  in  determining  whether  the  enhanced  light  intensity  in  SPP  can  be 
exploited  to  yield  large  optical  forces.  Woolf  et  al.  [33]  have  theoretically  analyzed  the  forces 
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generated  by  SPP’s  between  metal  slabs.  They  determined  that  semi-infinite  metal  waveguides 
separated  by  an  air  gap  (MIM  geometry  as  shown  by  Figure  1-5  (a))  give  rise  to  an  attraetive 
foree,  while  finite -thiekness  metal  waveguides  separated  by  an  air  gap  (IMIMI  geometry  as 
shown  by  Figure  1-5  (b))  yield  either  attraetive  or  repulsive  forees,  depending  on  the  mode.  The 
authors  foeus  eonsideration  on  modes  with  long-range  propagation  eharaeteristies,  for  whieh  the 
field  overlap  with  the  metal,  and  eonsequently  the  loss,  is  relatively  low. 


Figure  1-5  Optical  forces  generated  in  plasmonic  devices:  (a)  finite -thickness  metal 
waveguides  separated  by  a  dielectric  gap,  (b)  semi-infinite  metal  plates  separated  by  a 
dielectric  gap,  and  (c)  a  hybrid  plasmonic  system  with  a  dielectric  waveguide  and  a 

metallic  substrate. 


Two  papers  have  calculated  force  enhancements  in  hybrid  plasmonic  waveguides  consisting 
of  dielectric  waveguides  coupled  to  metallic  substrates  (Figure  1-5  (c)).  References  [34]  and  [35] 
show  that  the  force  between  a  silicon  waveguide  and  a  metal  substrate  is  up  to  an  order  of 
magnitude  larger  than  the  force  between  the  same  waveguide  and  a  dielectric  substrate,  given  the 
same  cross-sectional  power.  However,  this  conclusion  applies  only  for  relatively  small 
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waveguide-substrate  gaps.  The  plasmonic  enhaneement  is  non-resonant,  and  thus  can  be  used  for 
a  broad  range  of  wavelengths.  However,  due  to  ohmic  loss,  the  hybrid  plasmonic  mode  decays  in 
the  propagation  direction.  As  a  result,  the  enhancement  will  only  be  maintained  for  distances 
shorter  than  the  plasmon  decay  length,  typically  on  the  order  of  10 ’s  of  microns. 

In  Chapter  2,  the  effect  of  slow  light  on  optical  force  enhancement  will  be  examined.  The 
optical  force  is  proportional  to  the  light  intensity  in  the  propagating  mode.  Under  conditions  of 
fixed  power  input  to  a  guided-wave  mode,  the  light  intensity  scales  as  the  input  power  over  the 
group  velocity.  Consequently,  reducing  the  group  velocity  of  a  mode  increases  the  optical  force 
for  fixed  input  power.  In  the  extreme  case,  as  the  light  group  velocity  approaches  zero,  the 
optical  force  diverges. 

1.3  Applications  for  On-chip  Manipulation  of  Light  Signals 

As  discussed  above,  optical  forces  have  been  calculated  for  a  range  of  microphotonic  device 
geometries,  and  mechanical  motion  arising  from  optical  forces  has  been  measured  in  a  variety  of 
experiments.  This  section  will  discuss  some  of  the  applications  of  these  effects  in  on-chip 
information  processing,  such  as  filtering,  switching,  and  other  operations. 

Optomechanical  techniques  offer  particular  advantages  and  disadvantages  relative  to  other 
methods  for  tuning  the  response  of  a  microphotonic  system.  Mechanical  motion  results  in  a  large 
change  in  the  effective  index  of  a  microphotonic  structure,  as  compared  to  other  tuning 
mechanisms  such  as  thermal,  electrooptic,  and  carrier  injection.  Thermal  tuning  through  10 
degrees  kelvin,  for  example,  changes  the  index  of  silicon  by  less  than  0.07%.  Mechanical  motion 
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changes  the  refraetive  index  by  order  unity,  due  to  the  large  index  differenee  between  silieon  and 
air.  The  eharaeteristie  time  seale  for  optomeehanieal  effeets  is  determined  by  the  meehanieal 
motion.  If  the  deviee  is  operated  in  steady  state,  the  time  seale  is  determined  by  the  meehanieal 
period  (proportional  to  the  inverse  meehanieal  resonanee  frequeney).  This  number  sets  a  lower 
limit  on  the  time  seale  for  response.  If  the  optieal  signal  is  switehed  on  abruptly,  the  meehanieal 
displaeement  will  equilibrate  on  a  time  seale  determined  by  the  produet  of  the  meehanieal  period 
and  the  meehanieal  quality  faetor. 

With  these  general  features,  optomeehanieal  effeets  ean  be  used  to  eonstruet  novel,  tunable, 
on-ehip  optieal  eomponents.  These  inelude  deviees  for  wavelength  filtering,  self-adaptive 
mieroeavities,  slow-light  effeets,  and  tunable  lasers. 


1,3,1.  Tunable  Directional  Coupler 

One  example  of  a  waveguide  switeh  is  a  direetional  eoupler  [36].  The  standard  direetional 
eoupler  design  eonsists  of  two  parallel  waveguides.  Light  input  to  one  waveguide  transfers  to  the 
other  waveguide  over  a  eharaeteristie  length  seale  ealled  the  eoupling  length.  For  proper 
adjustment  of  the  deviee  length,  L,  all  input  light  is  transferred  to  the  drop  port.  Tuning  the 
effeetive  index  of  the  waveguides  ehanges  the  eoupling  length,  also  ehanging  the  fraetion  of 
power  that  is  output  at  the  through  and  drop  ports. 
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in  port 


through  port 


Figure  1-6;  Optically-tunable  directional  coupler,  in  which  the  effective  index  of  each 
arm  is  tuned  via  an  adjustable  slot  width.  Figure  adapted  from  [37]. 

Fong  et  al  [37]  have  used  optical  forces  to  demonstrate  an  all-optically  tunable  directional 
coupler.  Each  of  the  waveguides  is  a  slot  waveguide  (Figure  1-6).  An  optical  pump  signal  is  used 
to  induce  an  optical  force,  which  reduces  the  slot  width.  The  change  in  slot  width  causes  a 
change  in  effective  index  of  each  slot  waveguide,  changing  the  behavior  of  the  coupler.  The 
authors  demonstrated  that  the  output  power  of  a  probe  signal  at  the  drop  port  could  be  tuned 
from  0.1%  to  70%  using  1.3  mW  pump  power.  Dynamic  modulation  was  also  demonstrated. 


1,3,2,  Wavelength  Filtering 

Optical  forces  provide  a  method  to  achieve  wideband  tuning  of  microcavity  resonances,  with 
potential  applications  to  filtering  and  routing. 

In  an  early  paper,  Eichenfield  et  al.  [38]  demonstrated  an  all-optical  tunable  filter  based  on  a 
high-Q  silicon  nitride  microresonator  coupled  to  a  movable  silica  waveguide.  The  optical  force 
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due  to  a  strong  pump  beam  pulls  the  waveguide  toward  the  mieroresonator,  ehanging  the 
resonanee  frequeney  at  a  probe  wavelength.  The  pump  wavelength  was  tuned  to  change  the 
waveguide-resonator  separation,  shifting  the  probe  transmission  spectrum. 

Two  papers  have  demonstrated  tunable  fdters  based  on  power  tuning  of  a  pump  beam. 
Conceptually,  the  devices  used  are  similar  to  Figure  1-2  (b);  stacked  microrings  support  coupled 
optical  modes.  The  rings  are  mechanically  supported  by  hub-and-spoke  structures  (not  shown), 
which  can  be  tailored  to  control  the  effective  spring  constant  of  the  mechanical  modes. 

Rosenberg  et  al.  [29]  demonstrated  the  use  of  optical  forces  to  tune  the  microcavity 
resonance  in  a  double -ring,  silica  “spiderweb”  structure.  The  optical  force  changes  the  separation 
between  the  rings,  shifting  the  resonance.  Static  tuning  was  demonstrated  with  an  efficiency  of 
2.5nm/mW.  Dynamic  tuning  was  demonstrated  with  a  switching  time  less  than  200  ns. 

Wiederhecker  et  al.  [39]  later  used  optical  forces  to  tune  the  optical  resonance  of  a  silicon 
nitride  “double  wheel”  microcavity.  Tuning  of  30  nm  was  achieved  with  an  efficiency  of 
2.3nm/mW.  The  authors  argue  that  the  tuning  range  is  ultimately  limited  by  the  onset  of 
mechanical  regeneration. 

In  both  these  papers,  the  transmission  line  shape  is  a  dip,  or  rejection  filter.  It  has  been 
suggested  that  modified  devices  and  measurement  configurations  could  be  used  to  demonstrate 
related  functions,  such  as  routing,  switching,  buffering,  dispersion  compensation,  pulse  trapping 
and  release,  and  tunable  lasing  [29].  The  wide  tunability  suggests  these  functions  could  be 
achieved  over  the  full  telecom  C-  or  L-bands  [39]. 


13 


32 


1,3,3  Self-adaptive  Microcavity 

Rakich  et  al.  [40]  have  proposed  theoretieally  that  optieal  forces  can  be  used  to  design  self- 
adaptive  microcavities,  cavities  which  tune  themselves  to  the  incident  laser  frequency.  This 
behavior  is  achieved  by  careful  design  of  the  optomechanical  potential  in  a  modified,  double-ring 
system,  shown  in  Figure  1-7.  The  waveguide  on  the  right  is  used  to  couple  light  into  the  rings. 
The  force  between  rings  is  attractive  at  large  separations  and  repulsive  at  small  ones.  As  a  result, 
the  ring  separation  will  adjust  to  a  stable,  intermediate  separation  value.  Using  the  optical  forces 
between  an  additional  pair  of  coupled  waveguides  (shown  on  the  left)  to  further  shape  the 
optomechanical  potential,  the  system  can  be  designed  so  that  the  stable  ring  separation  always 
corresponds  to  a  resonance  frequency  equal  to  the  incident  laser  frequency. 


Figure  1-7  Schematic  design  of  a  self-tuning  resonator,  adapted  from  Reference  [40]. 


The  concept  of  a  self-aligning  microcavity  may  provide  a  route  to  more  robust  on-chip 
microcavity  devices  [41].  Temperature  variations,  fabrication  error,  and  other  factors  all 
influence  the  resonant  frequencies  of  microcavities.  Using  optomechanical  effects,  a  reference 
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laser  could  be  used  to  lock  multiple  microcavity  devices  on  a  chip  (such  as  fdters)  to  a  standard 
frequency. 

The  main  challenge  in  demonstrating  self-adaptive  behavior  is  fabrication  feasibility.  The 
design  of  Reference  [40]  requires  two  device  layers.  Each  layer  contains  at  least  one  waveguide 
and  microring,  and  all  devices  within  the  layer  must  move  together.  Creative  approaches  and/or 
modified  designs  may  be  required  to  demonstrate  self-adaptive  concepts  experimentally. 

1,3,4  Slow  Light 

Recent  papers  by  Painter  and  coworkers  have  suggested  the  use  of  optomechanical  effects  to 
achieve  slowing  of  light  [42,  43].  Motivated  by  the  dramatic  demonstration  of  stopped  light  in 
atomic  gasses  [44],  resulting  from  electromagnetically-induced  transparency  (BIT),  much  recent 
work  has  attempted  to  demonstrate  slow  light  in  practical,  on-chip  microphotonic  systems  [45- 
47].  Applications  of  slow  light,  such  as  buffering  and  optical  memories,  are  found  in  both 
classical  and  quantum  information  processing. 

The  maximum  pulse  delay  that  can  be  achieved  in  a  slow  light  system  is  fundamentally 
limited  by  the  lifetime  of  a  “dark  state,”  a  resonance  of  the  system  that  is  weakly  coupled  to  the 
outside  world.  The  dark-state  lifetime  scales  as  the  product  of  the  resonance  period  times  the 
quality  factor.  Painter  and  coworkers  have  suggested  the  use  of  a  mechanical  resonance  as  the 
dark  state;  the  mechanical  period  is  typically  much  longer  than  the  optical  period.  In  Reference 
[42],  the  authors  experimentally  demonstrated  a  reflection  spectrum  characteristic  of  BIT  in  a 
optomechanical  device.  The  structure  is  similar  to  Bigure  1-3  (b),  however,  the  microcavity  is 
formed  by  changing  the  hole  size  and  shape  near  the  center  of  the  device.  At  low  temperature 
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(8.7  kelvin),  an  optically  tunable  delay  of  50  nanoseconds  with  near-unity  optieal  transpareney 
was  aehieved.  A  method  for  tunable  eontrol  of  the  optomeehanieal  eonstant  in  time  would 
further  allow  trapping  and  release  of  finite-bandwidth  pulses. 


1,3,5  Tunable  Lasers 

The  ability  to  tune  the  operating  wavelength  of  a  semieonduetor  laser  by  mechanically 
modifying  the  optieal  eavity  is  an  intriguing  applieation  of  optomeehanieal  deviees.  Alegre  et  al. 
[48]  have  analyzed  laser  tuning  in  photonic  crystal  “zipper”  eavities  (similar  to  Figure  1-3  (a)). 
They  propose  to  use  two  optieal  resonances,  one  eorresponding  to  a  master  (or  pump)  frequeney, 
and  the  other  to  the  slave  (or  lasing)  frequeney.  The  optieal  foree  resulting  from  pump  light 
would  be  used  to  adjust  the  separation  between  two  nanobeams  that  make  up  the  zipper.  The 
change  in  separation  will  shift  the  frequeney  of  the  lasing  mode.  A  partieular  advantage  of  the 
zipper  strueture  is  the  ability  to  obtain  large  optomeehanieal  eoupling  eonstants  [48];  the  shift  in 
the  lasing  frequeney  seales  as  the  product  of  the  optomechanical  coupling  constants  for  the  pump 
and  lasing  modes.  However,  wavelength-tunable  lasers  eould  also  be  realized  through  other 
mieroseale  resonant  photonie  topologies  with  optomeehanieally  adjustable  air  gaps.  Tradeoffs 
between  eonsiderations  sueh  as  output  power,  wavelength  tunability,  frequeney  stability,  and 
mode  profile  will  affeet  the  ehoiee  of  design. 
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1,3,6,  Our  Contributions 

In  Chapter  2,  3  and  4,  more  applieations  of  optical  forces  will  be  discussed.  In  Chapter  2,  it 
will  be  proposed  that  optomechanical  response  also  offers  the  opportunity  to  control  or  regulate 
transmitted  power  levels.  Optical  forces  have  been  calculated  between  a  periodically-patterned, 
movable  silicon  waveguide  and  a  substrate.  For  a  given  wavelength  of  operation,  the  separation 
between  waveguide  and  substrate  will  determine  whether  light  propagates  or  falls  in  the  photonic 
bandgap  of  the  device.  For  proper  design,  increasing  the  optical  power  will  pull  the  waveguide 
toward  the  substrate,  moving  the  waveguide  into  the  bandgap  and  prohibiting  propagation.  This 
effect  should  limit  the  total  amount  of  power  that  can  be  transmitted  through  the  device.  The 
design  of  novel  devices  that  perform  the  reverse  function,  only  allowing  transmission  above 
certain  power  levels,  will  allow  the  regulation  of  power  levels  on  chip.  In  combination  with  gain 
elements,  optomechanical  elements  could  thus  provide  a  way  to  adaptively  shape  the  average 
power  of  signal  bit  streams. 

In  Chapter  3,  an  effective  Kerr  nonlinearity  will  be  studied  in  coupled  coupled-waveguide 
systems,  where  the  force  between  waveguides  can  be  either  attractive  or  repulsive.  The  essence 
of  an  optical  nonlinearity  is  that  the  propagation  characteristics  of  a  device  change  with  incident 
power.  This  is  the  case  for  optomechanical  effects.  Light  results  in  a  force  that  physically  moves 
the  device,  affecting  light  propagation.  The  effective  Kerr  nonlinearity  is  named  as  “mechanical 
Kerr  nonlinearity”  [49].  It  will  be  shown  that  either  sign  of  the  force  results  in  a  positive,  giant 
mechanical  Kerr  coefficient  and  how  the  magnitude  can  be  optimized  with  appropriate  device 
parameters.  The  mechanical  Kerr  coefficient  is  several  orders  of  magnitude  larger  than  the 
ordinary  Kerr  coefficient  of  silicon. 
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In  Chapter  4,  it  will  be  theoretieally  proposed  that  optomeehanical  effects  could  be  used  to 
achieve  all-optically  tunable  birefringence.  The  birefringence  in  a  system  of  two  coupled 
waveguides  (Figure  1-1  (a))  will  be  analyzed.  By  tuning  the  pump  power  injected  in  the 
waveguide,  the  separation  between  the  waveguides  can  be  adjusted.  The  change  in  separation 
affects  the  relative  phase  between  TE  and  TM  polarized  modes.  For  appropriate  values  of  the 
pump  power,  the  device  acts  as  a  polarization  converter,  changing  linear  to  circularly  polarized 
light. 

The  ability  to  achieve  mechanical  Kerr  effects  and  optically-tunable  polarization  rotation 
suggest  that  optomechanical  effects  can  be  used  in  a  variety  of  switching  configurations.  The  use 
of  non-resonant  devices  should  allow  wide  bandwidth  response. 

1.4  Applications  for  Assembly  of  Nanoparticles 

The  force  of  light  on  objects  provides  tremendous  flexibility  for  microscale  or  nanoscale 
manipulation.  As  discussed  above,  optical  forces  can  be  used  to  construct  novel,  tunable,  on-chip 
optical  components  for  information  processing.  Other  than  repositioning  elements  of  integrated 
photonics,  optical  forces  can  be  utilized  to  manipulate  particles.  One  famous  example  is  optical 
tweezer  [50-52],  which  use  the  optical  gradient  force  of  a  focused  laser. 

Recent  work  has  leveraged  the  strong  field  gradients  near  microphotonic  devices  for  particle 
trapping  [53-64].  However,  such  work  has  focused  on  trapping  single  or  few  particles.  Our  group 
have  proposed  to  use  optical  forces  near  microphotonic  devices  for  a  fundamentally  different 
purpose:  to  assemble  periodic  arrays  of  nanoparticles  resembling  synthetic,  reconfigurable  two 


dimensional  (2D)  crystals  [65].  This  approach,  called  light- assisted,  templated  self-assembly 
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(LATS),  exploits  photonic-crystal  slabs  to  create  resonantly  enhanced  optical  forces  orders  of 
magnitude  larger  than  radiation  pressure. 

The  LATS  process  is  shown  schematically  in  Figure  1-8.  Light  is  incident  from  below  on  a 
photonic-crystal  slab,  which  consists  of  a  silicon  device  layer  patterned  with  a  periodic  array  of 
air  holes.  The  slab  is  designed  to  support  guided-resonance  modes,  electromagnetic  modes  for 
which  the  light  intensity  near  the  slab  is  strongly  enhanced  [26].  The  previous  work  has 
theoretically  predicted  [65]  that  when  the  incident  laser  is  tuned  to  the  wavelength  of  a  guided- 
resonance  mode,  nanoparticles  will  be  attracted  toward  the  slab.  The  attractive,  optical  force 
arises  from  a  strong  electric-field  gradient  just  above  the  slab  surface.  In  addition,  the 
nanoparticles  will  experience  lateral  optical  forces  due  to  the  electromagnetic  field  structure  of 
the  guided  resonance  mode,  resulting  in  the  assembly  of  a  nanoparticle  array. 


Figure  1-8  Schematic  of  light-assisted,  templated  self-assembly  (LATS).  Incident  light 
from  below  excites  a  guided-resonance  mode  of  a  photonic-crystal  slab,  giving  rise  to 
optical  forces  on  nanoparticles  in  solution.  Under  the  influence  of  the  forces,  the 
nanoparticles  self  assemble  into  regular,  crystalline  patterns  [66]. 
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In  the  first  experiment  of  LATS  [66],  a  square  array  of  over  100  polystyrene  particles  was 
assembled  near  a  silicon  photonic-crystal  slab.  The  photonic  crystal  was  fabricated  on  silicon-on- 
insulator  wafer.  A  square  lattice  of  300-nm-diamter  holes  with  periodicity  of  860  nm  were 
etched  through  250-nm-deep  silicon  device  layer.  The  slab  was  immersed  in  a  suspension  of 
520-nm-diameter  polystyrene  particles  and  illuminated  from  below  with  a  laser.  The  particles 
floating  above  the  crystal  were  then  moved  into  the  holes  by  optical  forces,  forming  a  square 
crystal  lattice  measuring  13  pm  on  each  side. 

Figure  1-9  shows  snapshots  of  the  LATS  process.  When  the  laser  beam  is  turned  on, 
nanoparticles  are  attracted  to  the  slab  and  begin  to  occupy  sites  of  the  square  lattice  (Figure  1-9 
(a)).  As  time  progresses,  additional  particles  diffuse  into  the  region  where  the  beam  intensity  is 
high  and  begin  to  form  a  cluster  (Figure  1-9  (b)).  Eventually,  a  regular  array  of  particles  is 
formed  (Figure  1-9  (c)).  The  square  symmetry  of  the  assembled  particles  is  evident  from  the 
picture.  When  the  laser  beam  is  turned  off,  the  particles  immediately  begin  to  disperse  and 
diffuse  away  from  the  slab  (Figure  1-9  (d)). 
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a)  b)  c)  d) 


■> 


assembly  release 


Figure  1-9  Light-assisted,  templated  self-assembly  of  520  nm  diameter  partieles  above  a 
photonic -erystal  slab,  (a-c)  Sequential  snapshots  taken  with  the  light  beam  on.  (d) 
Snapshot  taken  after  the  beam  is  turned  off. 


Nanoparticle  arrays  assembled  using  LATS  can  be  viewed  as  “programmable  optical  matter” 
[67];  turning  the  laser  on  and  off  will  reversibly  assemble  or  disassemble  the  structure. 
Moreover,  exciting  different  resonance  modes  of  the  photonic  crystal,  by  adjusting  the 
wavelength  of  the  input  laser,  should  allow  different  crystalline  structures  to  be  formed  [65]. 

Light-driven  assembly  of  multiparticle  patterns  has  previously  been  achieved  using 
structured  light  fields  generated  by  interference  fringes,  holography,  spatial  light  modulators,  or 
other  methods  [68].  The  LATS  approach  differs  crucially  from  previous  work  in  that  it  exploits 
near-field,  rather  than  farfield  effects.  Rather  than  generating  a  structured  light  beam  via  firee- 
space  optics,  a  simple,  Gaussian  input  beam  is  used.  The  structured  light  field  responsible  for 
trapping  is  generated  by  the  interaction  of  light  with  the  photonic-crystal  device.  Ultimately, 
however,  LATS  could  be  carried  out  using  a  photonic-crystal  laser,  allowing  the  integration  of 
the  light  source  with  the  trapping  device  and  making  our  approach  highly  suitable  for  on-chip 
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integration.  Thus  a  wide  range  of  applieations  can  be  expected  from  all-optically  tunable 
photonic  devices,  to  materials  assembly,  to  biological  trapping  and  manipulation. 

In  the  previous  work,  a  simple,  square-lattice  design  for  the  photonic  crystal  is  considered.  It 
is  envisioned  that  the  LATS  technique  can  be  extended  to  assemble  larger  numbers  of  particles 
by  designing  the  photonic -crystal  slab  to  reduce  the  input  power  per  area  required  for  trapping. 
The  input  beam  can  then  be  spread  over  a  larger  area,  resulting  in  a  larger  cluster.  One  approach 
to  reducing  the  required  power  per  area  is  to  use  a  mode  with  a  higher  quality  factor.  Another 
approach  is  to  strongly  localize  the  field  in  the  trapping  regions.  With  these  two  principles,  novel 
lattices  will  be  proposed  in  Chapter  5  and  6.  By  judicious  design,  guided  resonance  modes  with 
high  quality  factors  and  mode  localization  can  be  achieved  in  the  novel  photonic  crystal  lattice.  It 
is  expected  that  the  new  design  will  allow  assembly  of  larger  numbers  of  particles. 

1.5  Force  Calculation 

Accurate  calculations  of  optical  forces  in  realistic,  3D  structures  can  be  performed  using 
electromagnetic  modeling  techniques. 


1.5.1.  Maxwell  Stress  Tensor 

One  method  for  force  calculation  is  to  evaluate  the  Maxwell  stress  tensor  (MST)  numerically 
[69].  This  general,  flexible  technique  allows  one  to  directly  calculate  the  force  based  on  the  full 
electromagnetic  field  distributions,  which  can  be  obtained  from  an  appropriate  full-vectorial 
electromagnetic  solver.  The  MST  method  has  been  used  to  evaluate  optical  forces  in  many 
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contexts  [22,  70-73].  To  compute  the  optical  forces  acting  on  a  movable  component,  the 
Maxwell  stress  tensor  is  numerically  integrated  over  a  closed  surface  surrounding  the  component: 

where  a  and  (5  indicate  direction  x,  y  or  z,  n  is  the  outward  normal  to  the  closed  surface  S,  and 
the  stress  tensor, 

■  E  +  -  H) . 

Typically,  the  time  scale  for  mechanical  response  is  longer  than  the  optical  period,  and  the  time- 
averaged  force  over  the  optical  period  may  be  considered  [73].  The  disadvantage  of  the  MST 
method  is  that  since  it  requires  the  knowledge  of  the  full  electromagnetic  fields,  it  is 
computationally  intensive  and  can  become  prohibitive  for  large  system  size. 


1,5,2  Derivative  Method 


It  has  been  shown  previously  [3]  that  for  infinitely  long  waveguides  with  no  radiation  loss, 
the  force  between  two  waveguides  can  be  calculated  from  the  derivative  of  the  eigenmode 
frequency  with  respect  to  separation.  This  relationship  follows  from  expressing  the  force  as  a 
derivative  of  the  electromagnetic  energy  U  with  respect  to  waveguide  separation,  g: 


dU 

d(Nh(o) 

--Nh— 

U  da 

dg 

k  dg 

k  dg 

,  ®  dg 

(1-2) 


where  N  is  the  number  of  photons,  co  is  the  eigenmode  frequency,  and  k  is  the  wave  vector. 
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An  optomechanical  coupling  constant  can  be  defined  as  .  The  value  of  goM  is 

proportional  to  the  force  per  photon  {FIN)  on  the  mechanical  system.  Larger  values  of  goM 
indicate  a  larger  change  in  the  optical  properties  of  the  system  for  a  given  mechanical 
displacement. 

It  is  often  convenient  to  consider  the  force  per  unit  length  of  the  waveguide,  normalized  to 
the  input  power  P.  Under  conditions  of  fixed  power  input  to  a  guided-wave  mode,  the 
electromagnetic  field  energy  U  scales  as  the  input  power  over  the  group  velocity:  U=PL/vg.  The 
normalized  force  is: 


F  _  1  dco 

PL  cov^  dg  ^ 


(1-3) 


Consequently,  reducing  the  group  velocity  of  a  mode  tends  to  increase  the  optical  force  per 
unit  waveguide  length  at  fixed  input  power.  The  force  may  alternately  be  written  as  a  derivative 
of  the  effective  index  Ueff  with  respect  to  separation  [13]: 


F  _  1  ^^eff 

PL  dg 


In  Equations  (1-3)  and  (1-4),  the  derivatives  are  taken  at  fixed  wave  vector. 

Derivative  methods  generally  assume  that  the  optomechanical  system  is  closed  and  does  not 
exchange  energy  with  the  surrounding  media.  The  force  between  two  microcavities  can  also  be 
calculated  as  a  derivative  of  the  cavity  mode  frequency  with  respect  to  separation,  provided  that 
the  cavity  quality  factor  is  large  [22,  74].  For  microresonator  systems,  the  optical  force  is 
proportional  to  the  electromagnetic  field  energy  stored  in  the  resonators. 
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Rakich  et  al.  have  shown  that  for  general  linear,  lossless  systems,  the  force  can  be  written  in 
terms  of  the  scattering  matrix  [75].  This  method  is  particularly  useful  for  systems  that  can  be 
accurately  described  by  coupled-mode  theory  [76],  such  as  systems  of  multiple  microresonators 
coupled  to  waveguides,  in  which  full-field  calculations  can  prove  prohibitive. 

Given  the  value  of  the  optical  force  at  a  given  optical  power  level,  a  resulting  static 
displacement  can  be  calculated  by  solving  the  corresponding  mechanical  problem.  A  typical 
approach  is  to  use  finite  element  software,  such  as  COMSOL,  for  this  purpose. 

In  the  following  chapters,  the  force  calculation  will  be  discussed  in  more  detail  in  specific 
systems. 

1.6  Thesis  Overview 

Since  the  initial  theoretical  studies  on  optical  forces  between  microphotonic  devices, 
experimental  progress  has  been  rapid.  The  goal  of  this  thesis  is  to  explore  applications  of  optical 
forces  in  micro-photonics  systems.  There  are  two  major  tasks.  One  is  to  exploit  effects  of  optical 
forces  in  highly  tunable,  on-chip  devices,  and  the  corresponding  discuss  will  be  in  Chapter  2,  3 
and  4.  The  other  task  is  to  utilize  optical  forces  for  light-assisted,  template  self-assembly  (LATS) 
of  nanoparticles.  A  novel  photonic  crystal  template  is  designed  in  Chapter  5  for  lower-power 
assembly.  In  Chapter  6,  experiments  are  carried  out  to  characterize  optical  performance  of  a  new 
photonic  crystal  lattice  potential  for  lower-power  LATS.  In  Chapter  7,  a  method  is  proposed  to 
predict  the  radiation  loss  of  coupled  resonator  optical  waveguides  in  photonic  crystal  slabs.  The 
method  may  facilitate  the  design  of  large-area,  coupled-cavity  modes  with  high  quality  factor 
that  nevertheless  couple  to  normally-incident  radiation  for  LATS. 
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In  Chapter  2,  the  attractive  optical  force  is  numerically  investigated  between  a  suspended 
one-dimensional  periodic  photonic  crystal  waveguide  and  underlying  substrate  in  a  silicon-on- 
insulator  platform.  It  is  shown  that  the  optical  force  can  be  enhanced  by  designing  the  waveguide 
cross  section  to  make  the  mode  approach  the  band  edge  or  substrate  light  line.  For  periodic 
waveguides,  the  optical  force  is  non-monotonic  with  waveguide-substrate  separation.  This  effect 
may  enable  the  design  of  compact,  integrated  optical  power  limiters. 

In  Chapter  3,  an  analytical  method  is  proposed  to  calculate  optical  forces  between  silicon 
waveguides  based  on  the  perturbation  of  effective  index  at  fixed  frequency.  The  method  is  used 
to  investigate  the  mechanical  Kerr  effect  in  a  coupled-waveguide  system  with  bipolar  forces.  It  is 
shown  that  positive  mechanical  Kerr  coefficient  results  from  either  an  attractive  or  repulsive 
force.  An  enhanced  mechanical  Kerr  coefficient  several  orders  of  magnitude  larger  than  the 
intrinsic  Kerr  coefficient  is  obtained  in  waveguides  for  which  the  optical  mode  approaches  the 
air  light  line,  given  appropriate  design  of  the  waveguide  dimensions. 

In  Chapter  4,  an  optomechanical  method  is  presented  to  tune  phase  and  group  birefringence 
in  parallel  silicon  strip  waveguides.  The  deformation  of  suspended,  parallel  strip  waveguides  due 
to  optical  forces  is  calculated.  The  frequency  and  polarization  of  the  pump  light  is  optimized  to 
obtain  a  9nm  deformation  for  an  optical  power  of  20  mW.  Widely  tunable  phase  and  group 
birefringence  can  be  achieved  by  varying  the  pump  power,  with  maximum  values  of  0.026  and 
0.13,  respectively.  The  giant  phase  birefringence  allows  linear  to  circular  polarization  conversion 
within  30pm  for  a  pump  power  of  67mW.  The  group  birefringence  gives  a  tunable  differential 
group  delay  of  6fs  between  orthogonal  polarizations.  The  tuning  performance  of  waveguides  is 
evaluated  with  different  cross  sections. 
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In  Chapter  5,  a  novel  photonie  erystal  lattiee  is  proposed  for  trapping  a  two-dimensional 
array  of  partieles.  The  lattiee  is  ereated  by  introdueing  a  reetangular  slot  in  eaeh  unit  eell  of  the 
Suzuki-Phase  lattiee  to  enhance  the  light  confinement  of  guided  resonance  modes.  Large  quality 
factors  on  the  order  of  10^  are  predicted  in  the  lattice.  A  significant  decrease  of  the  optical  power 
required  for  optical  trapping  can  be  achieved  compared  to  the  previous  design  based  on  square 
lattice. 

In  Chapter  6,  experiments  are  carried  out  to  optically  characterize  a  new  photonic  crystal 
structure.  This  photonic  crystal  device  is  generated  by  using  regular  graphite  lattice  as  the  base 
and  adding  a  slot  into  the  center  of  each  unit  cell  for  enhancing  field  confinement.  The  lattice  has 
been  optimized  to  obtain  guided  resonant  modes  with  high  quality  factor  over  4x10^.  The 
evolution  of  the  measured  wavelengths  and  quality  factors  follows  the  trend  predicted  by  the 
simulations. 

In  Chapter  7,  the  tight  binding  approximation  (TBA)  is  used  to  relate  the  intrinsic,  radiation 
loss  of  a  coupled  resonator  optical  waveguide  (CROW)  to  that  of  a  single  constituent  resonator 
within  a  light  cone  picture.  The  validity  of  the  TBA  is  verified  via  direct,  full-field  simulation  of 
CROWs  based  on  the  L2  photonic  crystal  cavity.  The  TBA  predicts  that  the  quality  factor  of  the 
CROW  increases  with  that  of  the  isolated  cavity.  Moreover,  the  results  provide  a  method  to 
design  CROWs  with  low  intrinsic  loss  across  the  entire  waveguide  band. 
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Chapter  2:  Effect  of  Periodicity  on  Optical  Forces 

2.1  Background 

In  the  last  few  years,  intense  researeh  has  been  earried  out  on  optieal  forees  induced  by  the 
strongly  enhanced  gradient  of  the  electromagnetic  field  close  to  micro-  and  nanophotonic  devices 
[1-11].  The  forces  provide  a  way  to  reconfigure  microphotonic  elements  with  experimentally 
measurable  deflection  [12-14],  This  mechanical  actuation  allows  a  number  of  important  device 
functionalities,  such  as  all-optically  controlled  filters  [15],  wavelength  routers  [16],  and 
polarization  rotators  [17],  as  well  as  optomechanical  Kerr  effects  [18].  One  structure  that  is 
particularly  amenable  to  fabrication  is  a  suspended,  microbridge  waveguide  over  a  substrate  [19]. 
For  appropriately  designed  devices,  it  has  been  predicted  and  verified  experimentally  that  the 
optical  force  is  dominant  over  other  effects,  such  as  thermo-optical  effects,  Casimir  and  van  der 
Waals  forces  [19,  20].  Previous  theoretical  and  experimental  work  has  demonstrated  that  light 
traveling  through  an  unpatterned  waveguide  induces  an  attractive  optical  force  between  the 
waveguide  and  substrate  [19,  21].  Introducing  a  periodic  pattern  into  the  waveguide  to  form  a 
one-dimensional  (ID)  periodic  photonic  crystal  should  introduce  additional  flexibility  in  tuning 
the  behavior  of  the  optical  force. 

Here  we  numerically  investigate  the  effects  of  periodicity  on  optical  forces  between  a 
microphotonic  waveguide  and  substrate.  While  previous  work  has  analyzed  ID-periodic 
photonic  crystal  microcavities  [22,  23],  we  consider  light  propagation  in  the  pass  band  of  a  ID- 
periodic  periodic  crystal.  We  show  that  the  optical  force  on  a  ID-periodic  photonic  crystal 
waveguide  is  enhanced  due  to  several  different  physical  effects:  delocalization  near  the  substrate 


light  line,  effective  reduction  of  the  waveguide  refractive  index,  and  slow-light  enhancement 
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near  the  band  gap.  We  find  that,  in  eontrast  to  unpatterned  waveguides,  the  optieal  foree  shows  a 
non-monotonie  dependenee  with  waveguide-substrate  separation.  The  foree  behavior  near  the 
photonic  band  edge  may  enable  the  design  of  optical  switches  that  transmit  only  below  a  certain 
power  threshold,  i.e.  power  limiters. 

2.2  Band  Edge  Effects 

We  investigate  the  optical  force  between  a  patterned,  ID-periodic  photonic  crystal  silicon 
waveguide  and  an  underlying  silica  substrate,  as  shown  in  Figure  2-1  (a).  In  the  paper,  the 
physical  wavelength  is  fixed  to  1550  nm.  We  assume  one  period  of  the  waveguide  is  a  cube  of 
side  length  a  with  a  circular  air  hole  of  radius  r=0.3a.  The  modes  of  the  waveguide-substrate 
system  can  be  characterized  as  even  or  odd  with  respect  to  the  mirror  plane  at  y=0,  according  to 
the  symmetry  of  the  vector  field  [24].  We  refer  to  the  y-even  mode  as  TM;  the  electric  field  is 
primarily  in  the  z-direction.  We  refer  to  the  y-odd  mode  as  TE;  the  electric  field  is  primarily  in 
the  y-direction. 

In  Figure  2-1  (b),  we  plot  the  dispersion  curves  for  the  two  lowest  TM  bands  and  the  lowest 

TE  band.  The  two  lowest  TM  bands  are  separated  by  a  band  gap.  For  TE  polarization,  the  second 

lowest  band  does  not  couple  strongly  to  a  linearly  polarized  source  and  is  not  shown  in  the  figure. 

We  use  the  MIT  Photonic  Bands  package  [25]  to  calculate  the  full  vectorial  electromagnetic 

fields  for  one  period  of  the  infinite  structure.  From  the  fields,  we  calculate  the  force  using  the 

Maxwell  Stress  Tensor  [26].  We  plot  the  force  per  unit  length  per  unit  power  (F/L/P),  where  we 

assume  P=  Uvg /L.  Vg  is  the  group  velocity,  and  U  is  the  electromagnetic  field  energy  [2].  We 

confine  our  calculation  to  the  values  of  the  period  a  for  which  the  modes  are  guided  (lie  below 
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the  gray  region  in  Figure  2- 1(b)).  All  three  of  the  modes  shown  give  rise  to  an  attraetive  optieal 
foree. 


(a)  (b) 


Figure  2-1  (a)  Patterned,  ID-periodie  photonie  erystal  Si  waveguide  with  period  a 
is  separated  from  the  Si02  substrate  by  a  distance  d  =100  nm.  r  is  the  hole  radius, 
(b)  Dispersion  relation  for  TM  (blue  solid  curve)  and  TE  (red  dashed  curve) 
modes  in  waveguide  with  r=0.3a,  and  h=w=a.  The  gray  region  shows  the 
substrate  light  cone.  The  point  where  the  TE/TM  mode  crosses  the  light  line  is 
indicated  by  an  orange/yellow  circle. 


Eigure  2-2  (a)  shows  force  per  unit  length  per  unit  power  as  a  function  of  the  side  length  a  for 
the  TM  modes  (blue  circles).  The  waveguide-substrate  separation  d=  \00  nm.  The  yellow  circle 
indicates  the  value  of  a  (~365nm)  for  which  the  TM  mode  enters  the  light  line.  The  blue  box 
indicates  the  range  of  a  for  which  1550nm  light  falls  within  the  band  gap  of  the  TM  modes.  In 
this  range,  no  light  can  propagate  along  the  waveguide,  and  thus,  no  optical  force  is  induced.  We 
observe  from  the  figure  that  the  force  increases  near  the  light  line  and  near  the  band  edge.  To 
understand  these  trends,  we  plot  the  force  for  an  unpattemed  strip  waveguide  with  the  same 


dimensions  but  no  air  holes  (black  squares).  Eor  the  strip  waveguide,  the  mode  crosses  the  light 
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line  at  a  value  of  a~280nm,  shown  by  the  yellow  eirele.  Like  the  patterned  waveguide,  the  foree 
inereases  as  the  mode  approaehes  the  light  line,  due  to  stronger  coupling  between  the  waveguide 
and  substrate.  Overall,  the  force  curve  for  the  patterned  waveguide  is  shifted  to  the  right  with 
respect  to  the  strip  waveguide.  This  effect  can  be  interpreted  as  a  reduction  of  the  effective  index 
due  to  the  air  holes.  We  plot  the  force  curve  for  a  strip  waveguide  with  reduced  index  equal  to 
2.79  with  pink  triangles.  Away  from  the  band  gap,  the  force  curve  for  the  patterned  waveguide 
follows  that  of  the  reduced-index  waveguide.  Near  the  band  gap,  the  enhancement  of  force  can 
be  interpreted  as  a  slow-light  effect  [1].  For  fixed  power,  the  energy  in  the  mode  increases 
proportionally  to  1/vg,  where  Vg  is  the  group  velocity.  Since  the  force  is  proportional  to  mode 
energy,  the  force  per  length  per  unit  power  goes  to  infinity  for  an  ideal  slow-light  waveguide. 


(a) 


Figure  2-2  The  attractive  force  per  length  per  unit  power  as  a  function  of  period  a 
for  (a)  TM  and  (b)  TE  modes  at  wavelength  k  =1550  nm.  One  period  of  the 
waveguide  (wv.)  is  a  cube  of  side  length  a  with  a  circular  air  hole  of  radius  r=0.3a. 
The  separation  d  between  the  suspended  waveguide  and  substrate  is  100  nm. 
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We  plot  the  foree  for  TE  modes  for  values  of  a  between  340  nm  and  375  nm  in  Figure  2-2  (b). 
The  light  wavelength  is  1550  nm.  As  for  TM  light,  the  foree  is  enhaneed  near  the  light  line  and 
band  edge.  Away  from  the  slow-light  region,  the  force  curve  follows  that  of  a  strip  waveguide 
with  reduced  refractive  index  equal  to  2.97,  (green  triangles). 

We  next  consider  the  dependence  of  the  force  on  the  waveguide-substrate  separation,  d.  In  a 
conventional  strip  waveguide,  the  force  increases  monotonically  with  decreasing  distance  [19]. 
In  the  patterned  waveguide,  however,  a  change  in  distance  can  move  a  mode  in  or  out  of  the 
photonic  band  gap,  yielding  quite  different  force  behavior. 

We  consider  a  waveguide  with  a  =  428nm.  From  Figure  2-2  (a),  we  can  infer  that  the  TM 
mode  lies  just  above  the  band  gap  for  d  =  lOOnm.  We  plot  the  force  as  a  function  of  separation  in 
Figure  2-3  (a)  (blue  triangles).  The  force  initially  decreases  with  increasing  distance,  then 
increases  as  the  mode  approaches  the  band  edge.  For  d  >  120nm,  the  mode  lies  in  the  band  gap. 
The  transmission,  and  thus  force,  are  zero. 

For  a  waveguide  with  a  =  408nm,  we  can  see  from  Figure  2-2  (a)  that  the  mode  lies  just 
below  the  band  gap  for  d  =  lOOnm.  The  force  (red  circles)  is  zero  for  distances  below  ~80nm, 
where  the  mode  falls  into  the  band  gap,  and  decreases  as  a  function  of  distance  for  d  >  80nm. 
The  behavior  should  allow  the  design  of  an  optical  power  limiter.  If  the  device  is  designed  with 
an  initial  waveguide-substrate  distance  >  80nm,  light  can  be  transmitted  through  the  waveguide. 
As  the  optical  power  is  increased,  the  optical  force  will  attract  the  waveguide  to  the  substrate 
more  strongly.  In  the  center  of  the  double-clamped  waveguide,  where  the  displacement  is  largest, 
the  local  mode  will  approach  or  enter  the  band  gap,  resulting  in  reflection,  and  reducing  the 
transmitted  power. 
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Figure  2-3  (b)  gives  the  force  induced  by  near  band  edge  TE  light  for  a  waveguide  with 
a=376  nm.  When  a  is  less  than  60  nm,  the  fundamental  mode  is  no  longer  guided,  and  we  set  the 
force  to  zero.  While  non-negligible  force  values  may  be  obtained  for  leaky  modes  with  long 
propagation  lengths,  the  force  will  decay  as  a  function  of  length  along  the  waveguide. 


Figure  2-3  The  attractive  force  per  length  per  unit  power  as  a  function  of  distance  d 
for  (a)  TM  and  (b)  TE  modes  at  wavelength  k  =1550  nm. 


The  displacement  can  be  obtained  from  the  calculated  force  values  in  Eigure  2-3  using  the 
analytical  expression  u^=F]^  l32Ewh^  for  a  doubly-clamped  beam  [27],  where  E  is  Young's 
modulus  for  bulk  silicon.  Eor  a  60-pm  long  beam  with  cross  section  376x376  nm  ,  if  the 
waveguide-substrate  separation  is  100  nm,  TE -polarized  light  with  power  P=50  mW  results  in  a 
E/E/P  of  2.8x10'  ,  and  the  displacement  at  the  center  of  the  waveguide  is  estimated  to  be  14.3  nm. 
We  have  verified  using  COMSOE  finite  element  simulations  that  the  displacement  of  the 
waveguide  with  hole  radius  r=  112.8  nm  is  15  nm.  The  discrepancy  is  likely  due  to  the  fact  that 
the  holes  modify  the  beam  flexibility  [28]. 
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2.3  Potential  Application 

In  conclusion,  we  have  studied  the  attraetive  optieal  foree  between  a  ID-periodie  photonie 
erystal  waveguide  and  underlying  silica  substrate.  The  optieal  foree  inereases  as  the  mode 
approaehes  the  light  line,  as  for  an  unpatterned  strip  waveguide.  However,  the  periodicity 
introduces  two  additional  meehanisms  for  enhancing  the  foree;  reduction  of  the  waveguide 
effective  index  and  slow  light  enhaneement  near  the  band  edge.  Moreover,  the  shift  of  the 
waveguide  mode  in  and  out  of  the  band  gap  with  ehanges  in  waveguide-substrate  separation  give 
rise  to  non-monotonie  foree  behavior  and  may  enable  eompaet  optieal  power  limiters.  It  is 
intriguing  to  eonsider  whether  optomeehanieal  eoupling  ean  be  used  more  widely  to  tailor  optieal 
power  response.  The  eoupling  of  dieleetrie  waveguides  to  substrates  with  loss  or  gain, 
ineorporating  plasmonie  or  quantum-confined  struetures,  may  enable  a  variety  of  devices  to 
tailor  power  response,  sueh  as  threshold-aetivated  switehes,  whieh  only  transmit  above  a  eertain 
power,  or  power  regulators,  which  self-adjust  to  maintain  the  output  transmission  within  a 
speeified  range. 
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Chapter  3:  Mechanical  Kerr  Nonlinearity 

3.1  Background 

In  the  last  few  years,  intense  researeh  has  been  earned  out  on  optieal  forees  indueed  by  the 
strongly  enhaneed  gradient  of  the  eleetromagnetie  field  elose  to  miero-  and  nanophotonie  deviees 
[1-14],  Optieal  forees  provide  a  novel  way  to  tune  the  properties  of  mierophotonie  deviees  by 
ehanging  the  separation  between  their  eomponents.  This  meehanieal  aetuation  allows  a  number 
of  important  deviee  funetionalities,  sueh  as  all-op tieally  eontrolled  filters  [15],  wavelength 
routers  [16],  polarization  rotators  [17],  and  power  limiters[18]. 

It  has  previously  been  suggested  that  optomeehanieal  eoupling  in  maeroseale  Fabry-Perot 
eavities  ean  be  viewed  as  an  effeetive  nonlinearity  [19].  The  eoneept  ean  be  generalized  to  other 
types  of  deviees  that  move  in  response  to  optieal  forees.  A  meehanieal  Kerr  effeet  has  been 
predieted  in  a  mierophotonie  Si  waveguide  with  a  suspended,  movable  seetion  [20].  The 
meehanieal  Kerr  eoeffieient  is  several  orders  of  magnitude  larger  than  the  intrinsie  Kerr 

eoeffieient  of  silieon.  The  origin  of  the  meehanieal  Kerr  effeet  lies  in  optieal  eoupling  between 
the  waveguide  and  substrate.  Coupling  results  in  an  attraetive  optieal  foree  that  pulls  the 
waveguide  eloser  to  the  substrate,  ehanging  the  effeetive  index  of  the  waveguide  mode.  Sinee  the 
optieal  foree  is  proportional  to  the  optieal  power  in  the  waveguide,  the  shift  in  effeetive  index 
depends  on  intensity;  for  this  system,  the  meehanieal  Kerr  eoeffieient  is  positive. 

It  is  intriguing  to  explore  whether  mierophotonie  deviees  ean  be  designed  to  tailor  the  value 
and  sign  of  the  meehanieal  Kerr  eoeffieient,  yielding  “eustomizable”  nonlinear  materials.  In  this 
paper,  we  examine  a  eoupled,  two-waveguide  system  in  whieh  the  forees  ean  be  either  attraetive 


or  repulsive.  We  show  that  the  meehanieal  Kerr  eoeffieient  seales  as  the  optieal  foree  squared. 
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As  a  result,  the  meehanieal  Kerr  eoeffieient  is  positive  for  either  an  attraetive  or  repulsive  foree. 
The  foree  ean  be  expressed  as  a  derivative  of  the  effeetive  index  with  respeet  to  separation  at 
fixed  frequeney.  We  ealeulate  the  meehanieal  Kerr  eoeffieient  numerieally  for  varying 
waveguide  eross  seetions.  We  find  that  w"  is  greatest  when  the  modes  are  neither  too  deloealized 
nor  too  eonfined  in  deformable  waveguides.  Under  these  eonditions,  the  foree  whieh  is  enhaneed 
by  strong  evaneseent-wave  eoupling  and  the  inereased  displaeement  give  rise  to  the  largest  shift 
in  effeetive  index. 

3.2  Coupled  Waveguides  Exhibiting  Bipolar  Optical  Forces 

We  eonsider  the  eoupled  waveguides  shown  in  Figure  3-1  (a).  Two  parallel,  silieon  strip 
waveguides  (refraetive  index  n  =  3.45)  are  separated  by  a  distanee  d.  Both  ends  of  the 
waveguides  rest  on  a  Si02  (refraetive  index  n=  1.5)  substrate.  The  suspended  seetion  has  length 
L.  Eaeh  waveguide  has  a  eross  seetion  of  dimensions  w^h.  From  previous  work  [2],  it  is  known 
that  optieal  eoupling  between  the  waveguides  gives  rise  to  a  meohanieal  foree.  The  magnitude  of 
the  foree  is  proportional  to  the  optieal  power,  and  the  sign  is  either  attraetive  or  repulsive, 
depending  on  the  relative  phase  of  light  in  the  two  waveguides.  The  foree  per  unit  length  in  the 
suspended  seetion  ean  be  ealeulated  from  the  full-veetorial  eigenmodes  of  the  eoupled 
waveguides  using  the  Maxwell  Stress  Tensor  method,  as  in  Referenee  [2]. 
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Figure  3-1  (a)  Two  coupled  silicon  waveguides,  each  with  cross  section  w^h,  separated  by 
a  distance  d  and  resting  on  a  Si02  substrate  with  suspended  length  L.  (b)  Dispersion 
relation  for  modes  of  the  suspended  section.  The  two  lowest  modes  with  electric  field 
vectors  primarily  in  the  z-direction  are  plotted  for  separations  J=100  nm  and  J=1000  nm. 
The  cross  section  of  each  waveguide  is  350x350  nm  .  The  gray  region  shows  the  light 

cone  for  air. 


Figure  3-1  (b)  shows  the  dispersion  relation  for  the  two  lowest-order  modes  that  have  an 
electric  field  primarily  in  the  z-direction.  Calculations  were  performed  using  the  MIT  Photonic 
Bands  (MPB)  package  [21].  To  illustrate  the  dependence  on  waveguide  separation,  the 
dispersion  relation  is  plotted  for  a  smaller  (lOOnm)  and  larger  (lOOOnm)  value  of  separation  d. 
For  the  lowest  mode,  the  electric  field  component  Ez  is  symmetric  in  the  two  waveguides.  The 
symmetric  mode  corresponds  to  an  attractive  force.  Decreasing  separation  shifts  the  dispersion 
relation  to  a  larger  k  value,  in  the  direction  shown  by  the  blue,  solid  arrow.  For  the  second  lowest 
mode,  E2  is  antisymmetric  and  the  force  is  repulsive.  Increasing  separation  also  shifts  the 
dispersion  relation  to  larger  k  value,  in  the  direction  shown  by  the  blue,  dashed  arrow.  Below,  we 
relate  the  shift  in  dispersion  relation  at  fixed  frequency  to  the  optical  force. 
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Figure  3-2  introduces  notation  to  describe  the  shift  in  dispersion  relation  co{k)  resulting  from  a 
small  change  in  waveguide  separation  dd.  The  solid  and  dashed  lines  are  the  initial  and  shifted 
bands,  respectively.  For  constant  frequency  operation  at  frequency  coq,  the  magnitude  of  the 
wave  vector  shifts  from  ko  to  k2.  If  we  consider  a  constant  wave  vector  ko,  the  dispersion  relation 
shifts  from  co  o  to  coj. 


Figure  3-2  Schematic  diagram  showing  a  shift  in  the  dispersion  relation  due  to  a  change  in 
separation  between  waveguides  d.  The  solid  and  dashed  lines  represent  the  initial  and 

shifted  dispersion  relation,  respectively. 


3.3  Calculation  of  Optical  Force 


For  infinitely  long  waveguides  with  no  radiation  loss,  the  optical  force  may  be  calculated  as  a 
derivative  of  the  eigenmode  frequency  with  respect  to  separation  [2]; 


F  

1  do) 

LP 

c  (0  dd 

(3-1) 
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where  F/{LF)  represents  the  foree  per  length  per  unit  power,  Ug  is  the  group  index,  and  c  is  the 
speed  of  light  in  vaeuum.  Note  that  the  derivative  is  taken  at  eonstant  k.  Equation  (3-1)  has  been 
shown  to  agree  well  with  direct  calculations  via  the  Maxwell  Stress  Tensor  (MST)  [22]  for  the 
coupled  waveguide  system  [2], 

By  substituting  (x>=ck/nsn  into  Equation  (3-1),  the  force  may  alternately  be  expressed  as  a 
function  of  the  effective  index  «eff  [23]; 


F 

LP 


1  "g 


c  dd 


(3-2) 


We  may  recast  the  force  in  terms  of  a  derivative  at  fixed  frequency  as  follows.  As  shown  by 
the  dashed  arrow  in  Figure  3-2,  for  fixed  wave  vector  ko,  a  small  change  in  separation  6d  moves 
the  eigenmode  from  (ko,  cod)  to  (ko,  co\),  corresponding  to  a  shift  of  Hen- 
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^1  1  ^ 


ckc 


{o)Q  +  6o))cOf^ 


(3-3) 


where  Sco=coi-coo.  For  fixed  frequency  coq,  the  shift  of  the  dispersion  relation  changes  the  wave 
vector  from  ko  to  k2.  The  change  of  iZeff  for  fixed  co  is: 


L„  =  —  iK- K)  ^ 


0)a 


CO, 


(3-4) 


where  Sk  =  ki-  ko.  By  comparing  Equation  (3-3)  and  (3-4),  we  find: 
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(3-5) 


If  the  band  is  shifted  by  a  small  amount,  the  ratio  Sco  ISk  can  be  estimated  by  -c/ug.  (The 
negative  sign  results  from  the  fact  that  Sco=  co\-coq  is  defined  to  be  negative  in  Figure  3-2. 
Equation  (3-5)  can  then  be  rewritten  as: 
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Substituting  Equation  (3-6)  into  Equation  (3-7)  and  taking  the  limit  as  Sd  ^0  yields  the  result 


F  ^  1  3^eff 
LP  c  dd 


(3-8) 


showing  that  the  optical  force  can  be  calculated  in  terms  of  the  effective  index  risi^d)  at  a  fixed 
optical  frequency.  Negative  values  here  correspond  to  attractive  forces.  Reference  [13]  used  a 
different  derivation,  based  on  energy  and  photon-number  conservation,  to  arrive  at  the  same 


expression. 

We  compared  force  values  calculated  with  Equation  (3-8)  to  those  calculated  using  the  MST 
method.  The  wavelength  was  fixed  to  1550  nm  and  the  waveguide  dimensions  to  h=w=350  nm. 
The  derivative  in  Equation  (3-8)  was  approximated  as  a  finite  difference,  using  values  of  neff(^0 
calculated  with  the  MPB  package.  Eigure  3-3  (a)  shows  neff(^0  for  the  symmetric  mode  (black 
line).  The  force  calculated  using  Equation  (3-8)  is  plotted  in  Eigure  3-3  (a)  (red  triangles).  It 
agrees  well  with  the  force  calculated  using  the  MST  method  (red  curve).  The  MST  was  evaluated 
using  the  full  vectorial  eigenmodes  obtained  from  MPB.  As  predicted  by  Equation  (3-8),  the 
force  is  attractive  (negative),  corresponding  to  a  decreasing 

Eigure  3-3  (b)  shows  ngff  (black  curve)  and  F/{LP)  for  the  antisymmetric  mode.  Eorces 
calculated  by  Equation  (3-8)  (red  triangles)  agree  well  with  values  calculated  using  the  MST 
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method  (red  eurve).  Note  that  the  value  of  /leff  first  deereases  and  then  inereases  as  a  function  of 
d.  Correspondingly,  the  derivative  turns  from  negative  to  positive,  and  the  optical 

force  changes  from  negative  (attractive)  to  positive  (repulsive)  as  a  function  of  separation.  We 
emphasize  that  the  sign  of  the  force  can  be  easily  inferred  from  the  slope  of  n^n. 


separation  d  [nm] 
(b) 


Figure  3-3  For  fixed  wavelength  1550  nm  and  waveguide  dimensions  h=w=350  nm;  (a) 
Symmetric  mode;  ngff  as  a  function  of  d  (black  curve);  F/(LP)  calculated  by  the  MSI 
method  (red  curve)  and  by  Equation  (3-8)  (red  triangles),  (b)  Antisymmetric  mode;  neff  as  a 
function  of  d  (black  curve);  F/{LP)  calculated  by  the  MST  method  (red  curve)  and  by 

Equation  (3-8)  (red  triangles). 


3.4  Calculation  of  Mechanical  Kerr  Coefficient 

The  mechanical  Kerr  effect  arises  from  the  change  in  effective  index  due  to  waveguide 
deformation.  We  consider  static  deformations.  The  optical  force  between  the  waveguides  is 
proportional  to  the  power,  P,  and  causes  the  waveguides  to  deform.  Eor  the  symmetric  mode,  the 
force  is  attractive  and  the  waveguides  are  pulled  together.  For  the  antisymmetric  mode,  the  force 
is  repulsive  and  the  waveguides  are  pushed  apart.  The  maximum  deformation  will  be  at  the 
center  of  the  waveguide,  since  the  ends  of  the  waveguides  are  fixed  to  the  silica  substrate.  We 
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consider  the  case  where  the  deformation  is  small  eompared  to  the  initial  separation.  The 
magnitude  of  the  deformation  is  also  mueh  smaller  than  the  length  of  the  suspended  region,  L. 
Under  deformation,  the  effeetive  index  ngff  varies  along  the  length  of  the  waveguides,  aeeording 
to  the  loeal  separation.  Following  Referenee  [10],  we  define  the  meehanieal  Kerr  eoeffieient  ri^ 

through  the  relation  An^j.f.(/)  =  ,  where  An^ff.(7)  is  the  ehange  in  the  effeetive  index  at  the 

eenter  of  the  suspended  region,  and  /  is  the  light  intensity.  Speeifieally,  we  take  I=P/A,  where  A 
is  the  mode  area,  equal  to  (j^y.  z) dyd^j max  ^  | 

Using  Equation  (3-8),  whieh  relates  the  force  to  the  derivative  of  the  effective  index  with 
respeet  to  separation,  we  ean  show  that  either  an  attraetive  or  repulsive  foree  gives  rise  to  a 
positive  meehanieal  Kerr  eoeffieient.  Figure  3-3  (a)  shows  that  for  the  symmetrie  mode,  neff 
deereases  with  d.  From  Equation  (3-8),  the  foree  is  negative  (attraetive)  and  pulls  the  waveguides 
closer  together.  The  separation  at  the  eenter  of  the  suspended  length  deereases,  resulting  in  an 
inerease  in  neff-  The  meehanieal  Kerr  eoeffieient  n"  is  thus  positive.  For  the  antisymmetrie 


mode,  in  the  region  where  neff  deereases  with  d,  the  same  argument  shows  that  is  positive.  In 

the  region  where  ^eff  inereases  with  d,  the  foree  is  positive  (repulsive),  whieh  increases  the 
separation  at  the  center  of  the  suspended  region  and  inereases  neff-  Again,  a  positive  meehanieal 
Kerr  eoeffieient  n”  is  obtained. 

The  meehanieal  Kerr  eoeffieient  ean  be  evaluated  numerieally  from  the  value  of  the  optieal 
force.  When  the  deformation  is  small  with  respeet  to  the  initial  separation,  the  ehange  of  the 
effective  index  for  a  signal  with  fixed  frequency  can  be  approximated  by: 


An  =  Adx^^^ 
dd 


(3-9) 
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where  AJ  is  the  ehange  in  separation,  equal  to  twiee  the  deformation  of  a  single  waveguide, 
which  can  be  modeled  as  a  double-clamped  beam  [24]; 

AJ  =  2x(Ff?/32£W),  (3-10) 

where  E  is  the  Young's  modulus  of  silicon.  By  manipulating  Equations  (3-8),  (3-9)  and  (3-10), 
we  get 

PL^ 


An^j.f 


1 

^^eff 

c 

dd 

m) 

\6Ehw 


.3  ’ 


(3-11) 


or 


An^ff  =  c 


The  mechanical  Kerr  coefficient  is 
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16EW 

Note  that  F/{LP)  depends  on  the  initial  separation  d  and  the  waveguide  cross  section  parameters, 
h  and  w.  The  parameter  n”  also  explicitly  depends  on  the  waveguide  geometry  parameters  h,  w, 
and  A. 

Note  that  Equation  (3-13)  provides  another  means  of  seeing  that  both  attractive  and  repulsive 


forces  give  positive  mechanical  Kerr  coefficient:  depends  on  {FI{LP))  . 
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3.5  Symmetric  Mode 


In  this  section,  we  discuss  how  to  design  the  waveguide  cross  section  parameters  h  and  w  to 


enhance  F/(LP)  and  for  the  symmetric  mode.  We  set  X  =  1550  nm,  L  =  20  pm,  initial 
separation  J  =  100  nm,  and  P  =  10  mW. 


Figure  3-4  Symmetric  mode:  (a)  Band  structure.  The  gray  region  represents  the  air  light 
cone.  The  red  line  shows  fixed  wavelength  1550  nm  (frequency  193.55  THz).  (b) 
Attractive  force  per  length  per  unit  power  F/{LP).  The  initial  separation  d=100  nm, 
suspended  length  L=20  pm,  optical  power  f  =  10  mW,  and  the  signal  wavelength  is  1550 

nm.  (c)  Mechanical  Kerr  coefficient  . 


Figure  3-4  (b)  shows  F/{LP)  for  the  symmetric  mode  in  waveguides  with  various  heights  and 
widths. 

The  strongest  attractive  force  is  shown  by  the  darkest  region,  which  corresponds  to 
waveguides  with  large  heights  and  small  widths.  To  understand  the  mechanism  behind  the  force 
enhancement,  we  plot  the  dispersion  relation  of  the  mode  in  Figure  3-4  (a).  For  a  waveguide  with 
h  =  260  nm  and  w  =  350  nm,  marked  by  the  circle  in  Figure  3-4  (b),  F/{LP)  has  a  relatively  large, 
negative  value  of  -4.1x10'  N/m/W.  The  corresponding  point  in  the  dispersion  relation  is  below 
and  close  to  the  air  light  line,  as  indicated  by  the  circle  in  Figure  3-4  (a). 
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Figure  3-4  (b)  shows  that  if  the  waveguide  cross  section  is  too  large,  the  attractive  force  is 
decreased  (yellow  region).  As  marked  by  the  square  in  Figure  3-4  (b),  when  h  =  540  nm  and  w  = 
350  nm,  F/{LP)  is  reduced  to  -2.7x10'  N/m/W.  The  mode  in  such  a  system  is  well  below  the  air 
light  line.  The  light  is  confined  in  the  two  waveguides,  with  minimal  optical  coupling  between 
them,  and  the  magnitude  of  the  force  is  decreased. 

The  white  region  of  Figure  3-4  (b)  shows  that  when  the  waveguide  cross  section  is  too  small, 
the  force  magnitude  is  less  than  10'  N/m/W.  For  example,  F/{LP)  on  the  waveguide  with  h  = 
110  nm  and  w  =  350  nm  (marked  by  the  star)  is  on  the  order  of  10'^  N/m/W.  Figure  3-4  (a)  shows 
that  the  mode  is  on  the  asymptote  of  the  light  line.  The  mode  energy  spreads  into  the  surrounding 
air,  giving  rise  to  a  small  force.  In  order  to  achieve  a  larger  attractive  force,  we  can  design  the 
waveguide  cross  section  to  make  the  optical  mode  approach  the  air  light  line  but  not  overlap  with 
the  light  line,  as  indicated  by  the  circle  in  Figure  3-4  (a). 

The  mechanical  Kerr  coefficient  can  be  determined  from  the  force  by  using  Equation  (3- 

13)  and  is  shown  in  Figure  3-4  (c).  The  biggest  in  Figure  3-4  (c)  is  over  40,000  times  bigger 
than  the  intrinsic  Kerr  coefficient  of  silicon,  equal  to  4.5x10'  m  /W.  The  value  of  is  largest 

when  the  optical  force  is  strong  and  the  waveguide  width  is  small.  A  thin  waveguide  is  more 
deformable  to  the  optical  force  along  the  y  direction.  For  example,  F/{LP)  for  a  waveguide  oih  = 
400  nm  and  w  =  200  nm  is  -6.7x10'  N/m/W,  which  has  the  same  order  of  magnitude  as  F/{LP)  = 
-7.2x  10'  N/m/W  for  a  waveguide  oih  =  400  nm  and  w  =  100  nm.  However,  by  reducing  w  from 
200  nm  to  100  nm,  we  can  increase  n”  from  2.4x1 0'^"^  m^/W  to  1.7x10''^  m^/W. 
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3.6  Antisymmetric  Mode 


We  calculated  F/{LP)  and  n”  for  the  antisymmetric  mode  for  waveguides  with  various 
heights  and  widths.  We  only  consider  guided  (non-leaky)  modes,  lying  underneath  the  light  cone. 

Figure  3-5  (b)  shows  F/{LP).  The  force  is  positive  (repulsive).  The  bright  region  indicates  the 
combination  of  waveguide  heights  and  widths  that  result  in  relatively  large  optical  force  for 
1550-nm  light.  The  blue  curve  indicates  the  conditions  for  mode  cutoff:  a  guided  mode  exists 
only  above  and  to  the  right  of  the  curve.  Large  forces  can  be  obtained  by  choosing  cross- 
sectional  parameters  on  the  upper,  right  boundary  of  the  curve.  For  example,  for  h  =  320  nm  and 
w  =  250  nm,  marked  by  the  circle  in  Figure  3-5  (b),  F/{LP)  is  9.5x10'^  N/m/W.  The 
corresponding  mode  in  the  dispersion  relation  lies  close  to  and  below  the  light  line,  as  indicated 
by  the  circle  in  Figure  3-5  (a). 

The  repulsive  force  decreases  for  waveguides  with  cross-sectional  dimensions  much  larger 
than  cutoff.  As  marked  by  the  square  in  Figure  3-5  (b),  when  h  =  540  nm  and  w  =  250  nm, 
F/(LP)  is  reduced  to  5.7x10'  N/m/W.  The  mode  energy  is  confined  in  the  waveguides,  and  the 
light  intensity  in  the  air  gap  is  reduced,  decreasing  the  optical  force.  If  the  waveguide  cross- 
sectional  dimension  is  too  small,  as  marked  by  the  star  in  Figure  3-5  (b)  for  /?  =  110  nm  and  w  = 
250  nm,  the  mode  becomes  leaky. 

Figure  3-5  (c)  shows  that  the  sign  of  n"  induced  by  a  repulsive  force  is  positive  for  various 
cross-sectional  dimensions.  In  order  to  increase  n"  ,  we  should  choose  dimensions  on  the  right 
boundary  of  the  blue  curve  and  reduce  w  to  make  the  waveguides  more  deformable.  For 
example,  if  we  choose  h=540  nm  and  w=140  nm,  is  approximately  8.1x10'  m  /W,  which  is 
more  than  150,000  times  larger  than  the  intrinsic  Kerr  coefficient  of  silicon. 
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Figure  3-5  Anti-symmetric  mode;  (a)  Band  structure.  The  gray  region  represents  the  air 
light  cone.  The  red  line  shows  fixed  wavelength  1550  nm  (frequency  193.55  THz).  (b) 
Repulsive  force  per  length  per  unit  power  F/{LP).  The  initial  separation  J=100  nm, 
suspended  length  L  =  20  pm,  optical  power  f  =  10  mW,  and  the  signal  wavelength  is  1550 

nm.  (c)  Mechanical  Kerr  coefficient  n”  . 


3.7  Discussion 

Above,  we  have  shown  that  for  fixed  waveguide  length  and  initial  separation,  we  can 
increase  the  mechanical  Kerr  coefficient  by  adjusting  the  waveguide  width  and  height.  From 

Equations  (3-8)  and  (3-9),  we  may  infer  that  is  proportional  to  the  product  of  the  optical 

force  F  and  the  waveguide  displacement  Ad.  The  mechanical  Kerr  nonlinearity  can  thus  be 

enhanced  by  i),  increasing  the  optical  force  and  ii),  making  the  waveguide  more  deformable. 

Figures  3-4  (b)  and  3-5  (b)  show  that  the  force  depends  on  waveguide  cross-section 

dimensions.  For  either  the  symmetric  or  anti-symmetric  modes,  decreasing  the  waveguide  height 

at  fixed  width  shifts  the  mode  closer  to  the  light  line.  For  the  symmetric  mode,  the  largest  force 

is  obtained  when  the  mode  is  close  to,  but  not  too  close  to,  the  light  line.  For  the  anti-symmetric 

mode,  the  largest  force  is  obtained  when  the  mode  approaches  the  cut-off  geometry  (blue  line), 

where  the  mode  crosses  the  light  line.  We  may  interpret  both  these  conditions  as  meaning  that 

57 


76 


the  modes  should  be  neither  too  deloealized  in  the  surrounding  air,  nor  too  eonfined  inside  the 
waveguides,  to  optimize  eoupling  between  the  waveguides  and  produee  a  large  force.  Meanwhile, 
the  waveguide  can  also  be  made  more  deformable  by  adjusting  waveguide  dimensions.  Equation 
(3-10)  indicates  that  Ar/  is  inversely  proportional  to  hv/  .  Overall,  the  largest  ri^  values  are 

obtained  for  small  values  of  the  width  w  (Figures  3-4  (c)  and  3-5  (c)). 

We  have  presented  results  for  both  the  symmetric  and  anti-symmetric  mode.  Both  modes 
give  positive  mechanical  Kerr  coefficients  n”  ,  and  the  maximum  values  are  both  on  the  order  of 
10'  m  /W.  Comparing  Figures  3-4  (c)  and  3-5  (c),  we  observe  that  the  symmetric  mode  has  a 
bigger,  brighter  region  representing  n"  greater  than  10“'"^  m^/W.  In  other  words,  to  realize  the 

same  value  of  n” ,  there  is  greater  flexibility  in  waveguide  dimensions  using  the  symmetric 

mode.  Another  reason  to  choose  the  symmetric  mode  is  that  exciting  the  antisymmetric  mode 
increases  device  complexity;  for  example,  it  has  been  demonstrated  that  a  Mach-Zender 
interferometer  can  be  used  to  control  the  phase  difference  between  light  incident  in  the  two 
waveguides  [10,  14]. 

Note  that  the  optical  frequency  is  many  orders  of  magnitude  larger  than  the  mechanical 
resonance  frequency,  which  is  in  the  MHz  range.  Because  the  instantaneous  optical  force 
oscillates  much  faster  than  the  mechanical  motion,  we  consider  the  effect  of  the  time-averaged, 
or  “smoothed,”  optical  force  on  the  waveguides,  where  the  average  is  taken  over  the  optical 
cycle.  Assume  the  incident  power  is  switched  on  and  remains  fixed.  Driven  by  the  optical  force, 
the  waveguides  will  start  to  deform.  A  steady-state  displacement  will  be  reached  after  a  period  of 
time  proportional  to  the  mechanical  quality  factor  divided  by  the  mechanical  resonance 
frequency.  In  this  paper,  we  focus  on  the  static  mechanical  Kerr  coefficient  n”  caused  by  the 
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steady-state  displacement.  However,  for  fixed  power,  larger  displacements  can  be  obtained  by 
modulating  the  incident  optical  power  at  the  mechanical  resonance  frequency  [20].  The 
displacement  will  vary  as  a  function  of  time  throughout  the  period  of  mechanical  vibration.  At 
the  peak  displacement  value,  the  instantaneous  reaches  its  maximum.  The  coefficient 
will  be  increased  by  a  factor  of  with  respect  to  the  case  of  static  displacement,  where 
Qmech  mechanical  quality  factor.  However,  if  the  average  mechanical  displacement  is  zero, 

the  average  value  of  will  also  be  zero. 

3.8  Summary 

We  have  investigated  mechanical  Kerr  effects  in  a  coupled-waveguide  system  exhibiting  both 
attractive  and  repulsive  forces.  The  optical  force  can  be  related  to  the  change  in  effective  index 
with  separation  at  frequency.  This  formulation  makes  it  clear  that  either  an  attractive  or  repulsive 
force  gives  rise  to  a  positive  mechanical  Kerr  coefficient.  Values  several  orders  of  magnitude 
larger  than  the  intrinsic  Kerr  coefficient  are  obtained  for  waveguides  in  which  the  optical  mode 
approaches  the  air  light  line. 

Since  the  mechanical  Kerr  effect  results  from  physical  motion,  the  time  scale  for  response  will 
be  slower  than  for  the  intrinsic  Kerr  nonlinearity.  At  the  same  time,  however,  the  magnitude  is 
orders  of  magnitude  larger.  This  suggests  the  potential  for  the  design  of  ultra-low  threshold, 
integrated  optical  devices  such  as  all-optical  transistors  and  isolators  [25]  that  use  the  mechanical 
Kerr  effect.  Such  devices  might  respond  to  the  average  power  in  an  optical  data  stream,  where 
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the  time-average  is  related  to  the  time  seale  for  mechanieal  response,  and  ultimately  find  use  in 
the  regulation  of  on-ehip  optieal  networks. 
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Chapter  4:  Large  Tuning  of  Birefringence 

4.1  Background 

Optical  forces  have  reeently  been  investigated  as  a  way  of  repositioning  mierophotonic 
elements  sueh  as  waveguides  [1-4]  and  mieroeavities  [5-8].  In  eoupled  waveguides  or 
mieroeavities,  the  foree  ean  show  either  attraetive  or  repulsive  behavior  depending  on  the 
relative  phase  of  the  ineident  light  in  eaeh  waveguide  or  mieroeavity.  As  a  result,  waveguide  (or 
eavity)  separation  ean  be  adjusted  by  tuning  the  power  and/or  polarization  of  input  light. 

Experiments  have  now  demonstrated  motion  indueed  by  optieal  forees  in  a  variety  of 
mierophotonie  systems.  Eiehenfield  et  al.  [9]  aehieved  mieron-seale  displaeement  of  a 
mierometer-scale  waveguide  evaneseently  eoupled  to  a  high-Q  optieal  mieroresonator  at 
milliwatt  optieal  powers.  Ei  et  al.  demonstrated  attraetive  and  repulsive  optieal  forees  between  a 
suspended  waveguide  and  dieleetric  substrate,  actuated  by  a  pulsed  light  souree  [10-12]. 
Attraetive  and  repulsive  forees  have  also  been  observed  in  staeked  mieroring  struetures  [13], 
where  Wiederheeker  et  al.  aehieved  statie  displaeements  up  to  12  nm.  Optomeehanieal  modes 
have  been  analyzed  [14]  and  measured  [15]  in  doubly-clamped  beam  structures  with  linear 
arrays  of  etehed  air  holes  forming  optieally  resonant  cavities.  As  the  field  of  optieal  forees  in 
integrated  photonie  devices  matures,  it  is  important  to  map  out  the  novel  functionalities  that  this 
relatively  new  form  of  meehanieal  aetuation  allows. 

In  this  seetion,  we  study  the  use  of  optieal  forees  to  aehieve  highly  tunable  birefringenee  in  a 
coupled  waveguide  system.  Birefringence  refers  to  differenees  in  propagation  between 
orthogonally  polarized  modes.  Phase  birefringenee,  the  differenee  between  the  phase  indiees. 
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induces  a  relative  phase  shift  between  orthogonal  eomponents  of  ineident  light,  whieh  may  be 
used  for  polarization  rotators  in  optieal  eommunieation  systems.  Group  birefringenee,  defined  as 
the  differenee  between  the  group  indiees  of  different  polarizations,  leads  to  a  polarization- 
dependent  time  delay  and  ean  be  exploited  to  realize  polarization-sensitive  delay  lines,  splitters, 
and  multiplexers.  It  has  previously  been  observed  that  silieon  waveguides  with  multiple  etehed 
slots  exhibit  giant  birefringenee  behavior,  with  a  group  index  differenee  as  large  as  1.5  between 
TE  and  TM  polarized  light  [16].  The  magnitude  of  birefringenee  depends  on  the  slot  width,  as 
well  as  other  design  parameters.  However,  after  the  deviee  is  fabrieated,  the  slot  width  is  fixed, 
as  is  the  value  of  the  birefringenee.  Here,  we  present  an  optomeehanieally-eontrolled  system 
whieh  resembles  a  suspended  one-slot  waveguide  with  an  adjustable  slot  width  that  depends  on 
the  optieal  foree  and  exhibits  highly  tunable  birefringenee. 

A  variety  of  methods  have  been  used  to  aehieve  tunable  birefringenee,  ineluding  thermal  and 
stress  tuning.  However,  the  aehievable  tuning  range  is  generally  small.  For  example,  thermal 
tuning  of  polymer-fdled  mierostruetured  fibers  [17,  18]  and  photonie  liquid  erystal  fibers  [19] 
yields  ehanges  in  birefringenee  on  the  order  of  10'"^.  For  miorophotonie  deviees,  Tsia  et  al.  [20] 
demonstrated  dynamie  eontrol  of  birefringenee  on  the  order  of  10'"^  in  a  silieon  waveguide  by 
eleetrieally  tuning  the  stress  on  the  waveguide  eore  using  an  integrated  piezoeleetrie  film. 
Reeently,  Kumar  et  al.  [21]  varied  the  air-eore  thickness  in  a  3D  hollow  waveguide  via  MEMS 
aeutation  to  offer  a  large  birefringenee  of  0.012  with  a  tuning  range  of  0.01.  This  result  points  to 
the  relatively  large  ehanges  in  birefringenee  that  ean  be  aehieved  using  meehanieal  motion  to 
physieally  ehange  the  waveguide  geometry. 
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In  this  work,  we  calculate  the  tunable  birefringence  that  can  be  achieved  by  using  an  optical 
force  to  adjust  the  separation  between  parallel,  suspended,  strip  silicon  waveguides.  The  force 
arising  from  the  evanescent  coupling  between  the  modes  of  the  two  waveguides  deforms  both 
waveguides,  changing  the  air-slot  width.  Such  behavior  results  in  all-optically  adjustable  giant 
birefringence.  In  Subsection  4.2,  we  characterize  the  two-waveguide  system,  and  optimize  the 
frequency  and  polarization  of  input  light  to  achieve  a  considerable  optical  force  at  a  moderate 
laser  power.  We  calculate  the  resultant  displacement  of  both  waveguides  based  on  the  finite- 
element  method  (FEM).  In  Subsections  4.3  and  4.4,  we  respectively  evaluate  the  performance  of 
the  optimized  optical  force  in  tuning  phase  and  group  birefringence.  We  show  that  the 
optomechanical  approach  is  characterized  by  widely  tunable  birefringence.  We  explore  the 
applications  of  the  two-waveguide  system  for  tunable  linear-to-circular  polarization  conversion 
and  polarization-dependent  delay. 

4.2  Design  Method 

Figure  4-1  (a)  shows  two  parallel  silicon  strip  waveguides  (refractive  index  n=3.45) 
separated  by  a  distance  d.  Both  ends  of  the  waveguides  rest  on  a  Si02  substrate,  with  a  free¬ 
standing  section  of  length  L.  Each  waveguide  has  a  cross  section  of  dimensions  w^h.  We 
consider  the  transverse  electric  (TE)  mode  with  electric  field  vector  primarily  parallel  to  the  air 
slot  (Ej;),  and  the  transverse  magnetic  (TM)  mode  with  electric  field  vector  primarily 
perpendicular  to  the  slot  (E^).  We  first  consider  the  case  where  the  cross  section  of  each 
waveguide  is  square  {w=h=a).  Figure  4-1  (b)  shows  the  dispersion  relation  calculated  by  the 
MIT  Photonic  Bands  package  for  the  lowest  TE  and  TM  mode  for  several  different  separations. 
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We  only  consider  the  guided  modes  lying  below  the  light  cone  (yellow).  The  insets  illustrate  the 
Ey  field  distribution  of  the  TE  mode  and  the  field  distribution  of  the  TM  mode,  respectively. 
TM-polarized  light  is  concentrated  in  the  air  region,  satisfying  the  continuity  condition  on  the 
normal  component  of  the  displacement  field  at  the  silicon-air  interface  [22,  23]. 


Figure  4-1  (a)  Two  coupled  Si  waveguides,  each  with  cross  section  w^h  separated  by  a 
distance  d,  rest  on  a  Si02  substrate  with  a  free-standing  section  of  length  L.  (b)  Dispersion 
relation  for  the  lowest-frequency  TE  mode  (solid  lines)  and  TM  mode  (dashed  lines)  of  the 
coupled  waveguides  for  several  separations.  Insets  respectively  show  the  Ey  field 
distribution  of  the  TE  mode  and  the  Ez  field  distribution  of  the  TM  mode  with  d=0.2a  at 
frequency  coa/2nc=0.\S  (darker  shades  correspond  to  larger  magnitudes  of  the  electric  field 
at  a  snapshot  in  time).  The  yellow  region  shows  the  light  cone. 


TE  and  TM  modes  have  different  propagation  characteristics.  The  group  birefringence  is 
obtained  by  taking  the  difference  between  the  group  indices  An  =  -  n™ ,  where  and  n™ 

are  the  group  indices  of  the  TE  and  TM  modes,  respectively,  at  a  fixed  frequency  co.  The  group 
index  is  calculated  as  =  c/  (d col  dkf^’™  and  is  proportional  to  the  reciprocal  of  the  slope 

of  the  dispersion  curve  in  Figure  4-1  (b).  Phase  birefringence  is  An^  -n^  ^  where  and 
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n™  respectively  denote  the  phase  indices  of  the  TE  and  TM  modes  at  a  fixed  frequency  co.  The 

phase  index  is  calculated  by  -  ck^^’™  /  0).  For  small  separation  {d=0.2a),  the  optical  fields 

in  the  two  strips  are  strongly  coupled.  The  bands  for  TE  and  TM  light  are  significantly  different. 
It  can  be  inferred  from  the  figure  that  group  and  phase  birefringence  are  relatively  large  in  some 
frequency  ranges.  If  the  separation  is  increased  {d=0Aa  or  d=0.6a),  the  coupling  between  the 
strips  becomes  weaker,  and  the  dispersion  curves  of  TE  and  TM  modes  are  closer  to  one  another. 
Group  and  phase  birefringence  are  reduced.  For  infinite  waveguide  separation  {d=cc),  no 
coupling  occurs.  Modes  propagate  independently  in  each  single-strip  waveguide.  The  bands  of 
the  TE  and  TM  modes  overlap,  and  the  birefringence  vanishes. 

Coupling  between  the  strips  gives  rise  to  an  optical  force  on  the  waveguides.  The  force 
depends  on  the  frequency  of  light  propagating  in  the  waveguides  (pump  frequency),  the  incident 
power,  and  the  waveguide  separation.  For  fixed  incident  power  and  waveguide  separation,  we 
optimize  the  force  magnitude  by  selecting  an  appropriate  frequency  and  polarization  of  pump 
light.  As  in  previous  work  [1],  we  calculate  the  force  as  a  function  of  frequency  using  the 
Maxwell  Stress  Tensor  formulation  [24],  which  gives  the  force  as  an  integral  expression  over  the 
electric  and  magnetic  fields  of  the  waveguides.  We  solve  for  the  full  vectorial  fields  using  the 
MIT  Photonic  Bands  package.  Figure  4-2  (a)  shows  the  force  magnitude  as  a  function  of  the 
pump  frequency  for  a  waveguide  separation  d  =  0.35a.  We  focus  on  the  dimensionless  frequency 
{(oa/lnc)  range  from  0.15  to  0.21.  For  lower  frequencies,  the  bands  are  closer  to  the  light  line, 
and  the  waveguide  modes  spread  out  in  the  surrounding  air.  For  frequencies  larger  than  0.21, 
there  are  at  least  two  symmetric  TE  (or  TM)  modes  for  a  single  frequency,  resulting  in 
multimodal  coupling.  We  compare  the  force  induced  by  TM  (red)  and  TE  (black)  modes.  The 


67 


strongest  force  is  obtained  for  TM -polarized  light  at  frequency  (Upa/27rc=0.165.  According  to 
Equation  (4-1), 


F  =  - 


1  dco 


03 


U 


(4-1) 


which  indicates  the  mechanical  force  F  is  proportional  to  the  derivative  of  frequency  co  with 
respect  to  waveguide  separation  ^  for  conserved  wave  vector  k  and  energy  U.  In  Figure  4-1  (b), 
we  see  that  for  a  fixed  wave  vector  k,  if  the  waveguide  separation  is  reduced,  the  frequency  of 
the  TM  mode  changes  more  than  that  of  the  TE  mode.  Therefore,  the  TM  mode  induces  a  bigger 
force.  In  the  remainder  of  the  paper,  we  will  assume  that  TM  pump  light  is  used  to  exert  force  on 
the  waveguides  and  deform  their  shape.  We  will  consider  the  effects  of  the  deformation  on  the 
propagation  of  both  TE  and  TM  signal  light. 


Figure  4-2  (a)  Normalized  force  per  unit  area  for  the  lowest-frequency  TE  (red)  and  TM 
(black)  modes  as  a  function  of  pump  light  frequency,  at  fixed  separation  d=Q35a.  (b) 
Force  per  unit  area  as  a  function  of  separation  (red  triangles),  at  optimized  frequency 
03pal2n  c=a/2p=0.165.  The  right  and  top  axes  are  in  physical  units  with  incident  power 
E’=20  mW,  w=h=a=263.5  nm,  andT=30  pm.  The  blue  solid  line  is  the  second-order 
polynomial  fit  of  the  force  per  unit  area. 


68 


87 


On  the  left  and  bottom  axes  of  Figure  4-2  (b),  we  plot  the  dimensionless  foree  per  unit  area 
(F/Lh)(a  c/P)  as  a  function  of  the  normalized  separation  dia  at  the  optimized  frequency 
(Upa/27rc=0.165.  The  symbol  c  represents  the  speed  of  light  in  vacuum,  andP  denotes  the  incident 
power.  Substituting  the  values  F’=20mW,  w=/z=a=263.5nm,  and  pump  wavelength  /lp=1597nm, 
we  plot  the  force  per  unit  area  versus  separation  on  the  right  and  top  axes  in  Figure  4-2  (b). 

We  perform  FEM  numerical  simulations  to  calculate  the  deformation  of  waveguides  caused 
by  the  attractive  force  using  COMSOL  Multiphysics  software.  We  simulate  the  full  three- 
dimensional  structure  consisting  of  two  silicon  waveguides,  each  with  a  square  cross  section  of 
263.5x263.5  nm  and  a  length  of  30pm.  The  ends  of  each  waveguide  are  fixed.  We  assume  an 
initial  waveguide  separation  J=0.35a=92.2  nm.  The  waveguide  ends  are  fixed  at  this  separation. 
The  optical  force  pulls  the  waveguides  toward  one  another.  As  the  waveguides  deform,  the 
magnitude  of  the  force  increases.  We  use  a  second-order  polynomial  to  approximate  the  force 
density  for  waveguide  separation  d  between  0.2a  and  0.35a,  as  shown  by  the  solid  blue  line  in 

•5 

Figure  4-2  (b).  For  the  mechanical  parameters  of  Si,  we  take  the  mass  density  p=2330  kg/m  and 
Young’s  modulus  £'=131  GPa.  Figure  4-3  shows  the  waveguide  deformation  as  a  function  of 
position  along  the  waveguides  with  £=30  pm,  w=h=a=263.5  nm,  d=92.2  nm,  and  £=20  mW. 
The  biggest  displacement  (9  nm)  is  obtained  at  the  waveguide  center.  The  smallest  distance 
between  the  two  waveguides  is  74.2  nm.  Following  [25],  we  estimate  the  van  der  Waals  force 
per  unit  area  to  be  approximately  25  N/m  between  two  parallel  infinitely  extended  Si  slabs  of 
thickness  263.5  nm  and  separation  74.2  nm.  For  this  separation,  the  van  der  Waals  force  is 
several  orders  of  magnitude  less  than  the  optically-induced  force  and  can  therefore  be  neglected. 
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It  is  interesting  to  note  that  experiments  have  been  done  in  a  similar  eoupled  system,  showing 
that  pulsed  light  ean  obtain  a  eomparable  displacement  at  a  lower  power  of  a  few  milliwatts  [10]. 


Figure  4-3  Displacement  of  the  suspended  section  of  each  waveguide  as  a  function  of 
position  along  the  waveguide.  The  cross-sectional  dimension  is  >v=/z=a=263.5nm,  the 
initial  waveguide  separation  is  r/=0.35a=92.2  nm,  the  pump  frequency  is 
(Upa/27ic=a/2p=0.165,  the  incident  power  is  P=20  mW,  and  the  suspended  length  is  T=30 
pm.  Note  that  the  x  andy  axes  differ  in  scale. 


4.3  Tunable  Phase  Birefringence 

The  phase  birefringence  is  a  function  of  the  separation  between  the  strips.  We  consider  the 
birefringence  experienced  by  light  at  a  signal  frequency  co.  We  assume  that  the  signal  power  is 
much  weaker  than  the  pump  power,  such  that  the  optical  force  due  to  the  signal  can  be  neglected. 
We  calculate  the  phase  birefringence  Aup  as  a  function  of  dimensionless  signal  frequency  roa/27ic 
for  different  separations  from  0.2a  to  0.35a,  as  shown  in  Figure  4-4  (a).  The  arrow  indicates  that 
at  signal  frequency  cOsal27i:c=0.n  (a=263.5  nm,  2.5=1550  nm),  the  absolute  value  of  increases 
as  the  waveguide  separation  decreases. 

When  the  waveguides  are  deformed  by  an  attractive  optical  force,  the  separation  varies  along 

the  waveguide  length  as  shown  in  Figure  4-3.  The  propagating  signal  experiences  a  spatially- 

varying  birefringence.  Figure  4-4  (b)  shows  the  phase  birefringence  as  a  function  of  position 
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along  waveguides  for  a  CW  pump  at  frequeney  (Upa/27rc=0.165  with  power  P=20  mW.  For  zero 
pump  power  (zero  force),  the  waveguides  are  parallel  and  the  birefringence  is  0.067  throughout 
the  whole  free-standing  section.  With  the  pump  on,  the  birefringence  is  increased.  The  maximum 
value  is  approximately  0.093  at  the  waveguide  center.  The  difference  between  the  maximum  and 
the  minimum  of  I I  is  0.026. 
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Figure  4-4  (a)  Phase  birefringence  as  a  function  of  signal  frequency  in  the  coupled 
waveguides  with  varying  separations.  The  arrow  shows  that  at  a  frequency  (X)sal2nc=Q.\l , 
the  absolute  value  of  An^  increases  as  the  waveguide  separation  d  decreases,  (b)  Phase 
birefringence  An^  as  a  function  of  position  along  the  waveguides.  The  initial  separation  is 
t/=0.35a=92.2  nm.  The  attractive  force  is  induced  by  CW  pump  light  at  frequency 
(X)pal2Tic=alXp=0.\65  and  power  P=20  mW.  The  difference  between  the  maximum  and  the 

minimum  of  An^  is  0.026. 


4,3,1  Tunable  Relative  Phase-shift 

A  potential  application  of  the  system  is  to  dynamically  adjust  the  relative  phase  shift  between 
TE  and  TM  modes.  The  total  phase  shift  between  TE  and  TM  modes  is  calculated  as; 

A0  =  \k^'^ (x)  -k™ {x)\ix  =  -  I nf  -  n™ \dx  =  —  | A |c/x  (4-2) 
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The  phase  shift  in  two  parallel  waveguides  with  a  separation  of  92.2  nm  is  estimated  to  be 
271x30  pmxQ. 067/1. 55  pm  =2.5971,  equal  to  0.5971  modulo  27i.  When  we  input  CW  pump  light 
with  a  power  of  20  mW,  the  waveguides  deform  and  the  phase  shift  is  ehanged  to  3.1  Iti,  or  1.1  Iti 
modulo  271.  The  tuning  range  is  about  0.57i.  If  we  reduee  the  waveguide  length  to  23.05  pm,  the 
total  relative  phase  shift  of  a  linear  signal  is  27i  at  the  output  of  the  un-deformed  deviee.  We  plot 
the  relative  phase  differenee  (phase  differenee  modulo  2n)  as  a  funetion  of  input  pump  power,  as 
shown  by  blaek  squares  in  Figure  4-5.  By  inereasing  the  pump  power,  the  phase  differenee  ean 
be  adjusted  from  0.127i  to  0.5 Iti.  Taking  the  pump  power  to  be  about  67  mW,  we  obtain  a  total 
relative  phase  shift  of  0.57i.  Thus,  an  ineoming  signal  of  linear  light  whose  £'-field  is  at  0.2571  to 
the  y  (or  z)  axes  will  emerge  from  the  deviee  eireularly  polarized. 

We  note  that  a  variety  of  approaehes  exist  for  polarization  eontrol  in  fiber  [26]  and  on-ehip 
mierophotonie  systems  [27-30].  In  mierophotonie  systems  in  partieular,  polarization  rotators 
have  been  designed  based  on  tapered  waveguide  eores  [27]  and  asymmetrie  waveguides  [28-30]. 
Referenee  [31]  presented  a  linear-to-eireular  polarization  eonverter  based  on  a  high-Q  dieleetrie 
mieroring  resonator.  In  eontrast  to  these  approaehes,  we  take  advantage  of  meehanieal  motion  to 
yield  a  widely  tunable  phase  shift  that  is  also  broadband  (non-resonant)  in  response.  For  a  deviee 
with  fixed  length,  the  phase  birefringenee  and  relative  phase  shift  ean  be  dynamieally  adjusted 
by  varying  pump  light  power.  Unlike  other  polarization  eontrol  sehemes  using  semieonduetor 
optieal  amplifiers  [32,  33]  or  eleetrooptie  materials,  our  approaeh  is  fully  eompatible  with  SOI 
materials  systems. 
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Figure  4-5  Tuning  the  relative  phase  shift  by  inereasing  CW  pump  power  from  20  mW  to 
70  mW  (blaek  squares).  The  waveguide  length  is  23.05  pm.  The  red  line  is  a  2nd  -order 
polynomial  fit  of  the  phase  shift.  A  power  of  approximately  67  mW  yields  a  0.57i  phase 

difference. 


4,3,2  Other  tuning  configurations 

The  dependence  of  phase  birefringence  on  separation  can  be  changed  by  altering  the 
waveguide  cross  section.  Figure  4-6  (a)  depicts  the  birefringence  versus  signal  frequency  in  two 
a^2a  waveguides.  The  arrow  indicates  that  the  value  of  Arip  decreases  with  decreasing  separation 
for  frequencies  from  0.1 1  to  0.16.  The  tuning  trend  is  opposite  from  that  in  the  system  with  two 
square  waveguides.  We  can  understand  the  trend  by  considering  two  extreme  cases.  For  zero 
separation,  light  propagates  in  a  square  waveguide  of  cross  section  2a^2a,  which  has  no 
birefringence  due  to  spatial  symmetry.  For  infinite  separation,  light  propagates  in  a  single 
rectangular  waveguide  asymmetric  with  respect  to  the  y  and  z  axes,  which  has  large 
birefringence.  Figure  4-6  (b)  shows  another  waveguide  structure  with  a  cross  section  ax0.5a.  The 
effective  index  of  the  TE  mode  is  less  than  that  of  the  TM  mode.  The  absolute  value  of  Aup 
increases  with  decreasing  separation  for  frequencies  from  0.20  to  0.25. 


73 


92 


008  0  09  010  011  012  013  014  015  016  020  021  022  023  024  025 

t,)a/2iic  .■.a/24C 

(a)  (b) 

Figure  4-6  Phase  birefringenee  Aup  as  a  function  of  signal  frequency  in  the  two-waveguide 
system  with  varying  separations,  (a)  The  cross  section  of  each  waveguide  is  ax2a.  The 
arrow  shows  that  Aup  decreases  with  decreasing  separation,  (b)  The  cross  section  of  each 
waveguide  is  ax 0.5 a.  The  absolute  value  of  Aup  increases  with  decreasing  separation. 

4.4  Tunable  Group  Birefringence 

The  group  birefringence  also  depends  on  waveguide  separation  and  can  be  tuned  by 
deforming  the  waveguides  via  an  optical  force.  We  plot  the  group  birefringence  Aug  as  a  function 
of  dimensionless  signal  frequency  coallnc,  for  different  separations  from  0.2a  to  0.35a  in  Figure 
4-7  (a).  The  birefringence  increases  with  decreasing  separation.  The  arrow  indicates  that  the 
group  birefringence  Aug  increases  fastest  at  a  frequency  (Usa/27rc=0.17  (a=263.5  nm,  2^=1550 
nm).  We  consider  deformation  due  to  an  attractive  optical  force  for  CW  pump  light  at  frequency 
(Upa/27rc=0.165  and  power  P=2Q  mW,  as  above.  Figure  4-7  (b)  shows  the  group  birefringence  of 
signal  light  as  a  function  of  length  along  the  waveguides  for  T=30  pm.  The  initial  value  1.09 
corresponds  to  the  birefringence  of  the  undeformed  waveguides.  The  peak  of  the  birefringence  is 
approximately  1.22  at  the  waveguide  center.  The  difference  between  the  maximum  and  the 
minimum  of  Aug  is  0.13. 
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Figure  4-7  (a)  Group  birefringence  Ang  as  a  function  of  signal  frequency,  in  the 
coupled  waveguides  with  varying  separations.  The  arrow  shows  that  as  the 
waveguide  separation  d  increases,  the  group  birefringence  Arig  increases  fastest  at 
a  frequency  cOsa/27j:c=0.\7 .  (b)  Group  birefringence  Ang  as  a  function  of  position 
along  the  waveguides.  The  waveguides  are  deformed  by  an  attractive  force 
induced  by  a  CW  pump  at  frequency  a>pa/27rc=a/2p=0.l65,  with  a  power  ^=20 
mW.  The  difference  between  the  maximum  and  the  minimum  of  Aug  is  0.13. 


The  system  can  be  used  to  tune  the  polarization-selective  delay  between  TE  and  TM  pulses. 
The  delay  is  tuned  by  adjusting  the  pump  power.  The  differential  time  delay  between  TE  and  TM 
pulses  is  calculated  as: 


(4-3) 


g 


g 


Eor  two  parallel  waveguides  without  deformation,  the  differential  time  delay  is 

o 

1.09x30pm/3xl0  =109  fs.  By  inputting  CW  pump  light  at  a  power  of  20  mW,  we  increase  the 
differential  time  delay  to  1 15  fs. 

If  we  consider  applications  where  the  pulse  width  is  comparable  to  the  tunable  delay  range, 
we  can  take  the  pulse  width  7=6  fs.  We  can  estimate  the  dispersion  length  as  L^=T^ !  [34], 
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where  the  dispersion  eoefficient  =  d^k  I  dco^  is  2.2x10'^^  s^/m.  The  parameter  T/)  is  around  164 

mierometers,  mueh  larger  than  the  waveguide  length.  We  ean  thus  assume  that  minimal 
distortion  oeeurs  due  to  propagation  through  the  30-pm  long  suspended  seetion. 

The  tuning  range  ean  be  inereased  by  inereasing  the  pump  power.  The  tuning  range  ean  also 
be  inereased  by  fabrieating  longer  suspended  seetions,  whieh  deform  more  easily  for  fixed 
power,  or  by  fabrieating  multiple  suspended  seetions  in  series. 

Note  that  this  analysis  also  applies  to  vibrating  waveguides  aetuated  by  pulsed  pump  light. 
The  dynamie  system  experienees  similar  displaeement  as  the  statie  system  for  mueh  smaller 
pump  powers.  While  the  waveguide  beam  vibrates  in  a  period  of  a  few  hundred  nanoseeonds,  the 
transmission  time  of  a  signal  pulse  through  the  30pm  waveguide  is  on  the  order  of  a 
femtosecond.  The  beam  can  thus  be  treated  as  static  in  the  calculation  of  time  delay,  provided 
that  the  signal  pulse  is  much  shorter  than  the  vibration  period  in  time. 

4.5  Thermal  Effects 

In  addition  to  optical  forces,  thermal  effects  may  also  contribute  to  waveguide  deformation. 
The  waveguides  absorb  propagating  light,  experience  thermal  expansion,  and  are  mechanically 
deformed.  We  performed  FEM  (finite  element  method)  simulations  to  calculate  the  displacement 
caused  by  the  photo-thermo-mechanical  effect  using  COMSOL  Multiphysics  software.  We 
considered  two  silicon  waveguides  with  identical  cross  section  263.5x263.5  nm  and  separation 
92.2  nm  sitting  on  top  of  a  silicon  dioxide  substrate,  as  shown  in  Figure  4-1  (a).  The  free-section 
length  of  the  beams  is  30  pm.  We  used  a  similar  simulation  method  to  that  introduced  in  the 
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supplementary  methods  of  Reference  [10].  The  beams  are  assumed  to  be  thermally  isolated, 
except  where  they  contact  the  substrate.  The  boundaries  of  the  silicon  dioxide  substrate  are 
assumed  to  be  of  constant  temperature  300  K.  The  absorption  constant  of  silicon  a  is  estimated  to 
be  0.0079  cm''  at  1.597  pm  pump  wavelength  by  using  the  equation  a=‘\Tik/X  and  the  value  of 
the  extinction  coefficient  k.  With  20mW  incident  optical  power,  the  absorbed  power  is 
approximately  0.474  pW.  After  using  the  RF  module  to  solve  for  the  pump  mode  profde,  we 
define  a  resistive  heating  power  of  0.474  pW  as  the  heating  source  in  the  thermal  module.  We 
obtain  a  calculated  temperature  profile  that  shows  a  peak  temperature  T-3 00=0. 097  K  near  the 
center  of  the  beams.  The  thermal  simulation  is  followed  by  a  mechanical  simulation  to  solve  for 
the  mechanical  stress  that  results  from  the  thermal  load.  We  are  primarily  concerned  with  the 
mechanical  displacement  in  the  z  direction,  which  may  introduce  inaccuracy  into  the 
birefringence  tuning.  A  peak  z-displacement  of  0.392  pm  is  reached  at  the  waveguide  center.  The 
peak  displacement  in  the  y-direction  was  <5  pm.  The  thermally  induced  displacement  is  thus 
negligible  compared  with  the  displacement  due  to  the  optical  force  of  9  nm. 

Considering  that  the  thermo-optical  coefficient  1.86x10'"'K''  for  silicon,  we  estimate  that  the 
refractive  index  change  due  to  our  calculated  temperature  increase  is  on  the  order  of  10'^,  which 
is  far  smaller  than  the  index  change  caused  by  optomechanical  deformation. 

We  conclude  that  thermal  effects  can  be  neglected  in  the  system  studied  here,  in  agreement 
with  the  conclusions  found  in  Reference  [10]  for  a  similar  experimental  system. 
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4.6  Summary 

In  this  work,  we  present  an  optomeehanieally-eontrolled  system  whieh  resembles  a 
suspended  one-slot  waveguide  with  adjustable  slot  width  depending  on  the  optieal  foree.  The 
ehange  in  slot  width  eorresponds  to  a  ehange  in  the  symmetry  of  the  eross-seetion  of  the  deviee, 
and  therefore  tunes  both  phase  and  group  birefringenee.  We  optimize  the  polarization  and 
frequeney  of  eontinuous  pump  light  to  obtain  a  eonsiderable  foree  at  low  operating  power.  The 
deformed  30  pm  waveguide  shows  a  maximum  displacement  of  about  9  nm  at  the  center  for  a 
power  of  20  mW.  Correspondingly,  the  phase  (group)  birefringence  at  the  waveguide  center  is 
different  from  that  at  the  clamped  ends  by  0.026  (0.13),  which  is  larger  than  previous  records. 
When  used  to  tune  the  relative  phase  of  two  orthogonally  polarized  light  signals,  the  device  has  a 
tuning  range  of  around  7i/2.  If  the  orthogonally  polarized  components  have  identical  amplitudes, 
it  is  possible  to  rotate  the  polarization  state  from  linear  to  circular  in  a  distance  of  30  pm.  We 
show  further  that  the  dependence  of  phase  birefringence  on  separation  can  be  changed  by 
altering  the  waveguide  cross  section.  When  used  as  a  polarization-dependent  tunable  delay 
element,  the  two-waveguide  system  has  a  tuning  range  of  6fs.  We  note  that  the  birefringence 
tuning  trends  studied  here  apply  not  only  to  a  two-beam  system  deformed  by  an  optical  force,  but 
also  by  other  actuation  methods,  including  static  electrostatic  forces  used  in  micro-electro- 
mechanical  systems  (MEMS). 

In  the  future,  the  potential  to  tune  birefringence  via  an  optical  force  may  be  applied  to  design 
novel  microring  resonator  systems.  Waveguide-coupled  microrings  exhibit  a  strong  (static) 
birefringence,  which  can  be  used  to  realize  polarization  converters  and  polarization-sensitive 
delay  devices  [35-37].  Recent  experiments  have  shown  that  the  separation  between  vertically- 
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stacked  microrings  can  be  adjusted  via  optical  forces  [13].  This  capability  may  in  the  future  be 
exploited  to  design  highly-tunable  birefringent  devices  along  similar  lines  as  studied  here. 
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Chapter  5:  Trapping  in  a  Slot-Suzuki-phase  Lattice 

5.1  Background 

Optical  trapping  of  objects  with  micro-  and  nano-scale  dimensions  has  opened  up  novel 
opportunities  in  areas  from  physios  to  biology  [1-4],  While  optioal  tweezers  make  use  of  firee- 
spaoe  beams,  reoent  researoh  has  explored  the  use  of  integrated  optioal  devioes  for  partiole 
trapping  [5].  The  strong  evanesoent  fields  in  suoh  mioron-sized  struotures  enable  trapping  on  the 
sub-wavelength  soale.  Designs  based  on  plasmonios  [6],  dieleotrio  waveguides  [7],  and  oavities 
[8-11]  have  been  proposed  and  demonstrated  for  trapping  of  single  partioles.  In  a  reoent  work 
[12],  we  have  proposed  a  method  to  assemble  2D  arrays  of  partioles  with  partioular,  designable 
patterns.  This  prooess,  whioh  we  oall  light  assisted,  templated  self  assembly,  relies  on  the 
struotured  light  fields  above  a  photonio  orystal  slab  to  oreate  an  array  of  partiole  trapping 
looations.  By  ohanging  the  wavelength  or  polarization  inoident  on  the  photonio  orystal,  the  fields, 
and  henoe  the  trapping  looations,  oan  be  reoonfigured.  This  prooess  oan  be  used  either  to  oreate 
tunable  photonio  filters  or  to  fabrioate  ordered  patterns  of  nanopartioles.  In  our  previous  work,  as 
a  proof  of  oonoept,  we  oonsidered  a  simple,  square-lattioe  design  for  the  photonio  orystal.  The 
estimated  power  for  trapping  was  ImW/unit  oell,  limiting  the  feasibility  of  the  method.  However, 
it  is  known  that  in  photonio-orystal  miorooavity  struotures,  the  introduotion  of  a  slot  into  the 
oavity  oan  provide  high  field  oonfmement,  reduoing  the  power  required  for  single-partiole 
trapping  [10,  11].  In  this  work,  we  propose  a  novel  2D  photonio  orystal  that  exploits  slot 
oonfinement  to  reduoe  the  optioal  power  for  trapping  of  partiole  arrays  by  two  orders  of 
magnitude. 
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Below,  we  present  our  novel  photonie  erystal  strueture,  designed  for  low-power  optieal 
trapping  of  partiele  arrays.  The  2D  periodie  strueture  is  ereated  by  introdueing  a  reetangular  slot 
at  the  eenter  of  eaeh  unit  eell  of  the  Suzuki-phase  lattiee  [13].  When  the  ineident  light 
wavelength  is  tuned  to  exeite  a  guided  resonanee  mode  (GRM),  the  optieal  field  is  eonfmed  and 
enhaneed  within  the  slot.  Meanwhile,  the  quality  (Q)  faetor  is  increased  by  orders  of  magnitude. 
Our  numerical  simulations  show  that  for  fixed  incident  power,  the  optical  force  on  a  particle 
above  the  lattice  is  also  enhanced  by  orders  of  magnitude,  significantly  reducing  the  power 
required  for  stable  trapping.  For  an  ideal  device,  stable  trapping  of  25-nm  radius  particles 
requires  a  power  as  low  as  27  pW  per  unit  cell,  around  40  times  smaller  than  for  the  simple 
square-lattice  structures  studied  in  our  previous  work  [12].  We  predict  that  a  particle  with 
dimensions  smaller  than  the  slot  will  be  pulled  inside  and  stably  trapped,  with  a  required  power 
of  only  3  pW  per  unit  cell,  300  times  smaller  than  in  our  previous  work.  The  novel  lattice 
structure  and  corresponding  low  power  requirements  also  open  up  an  opportunity  for 
experimental  implementation  in  active  structures.  Using  a  2D  photonic  crystal  surface-emitting 
laser  based  on  the  slot  Suzuki-phase  design,  we  expect  that  self-adaptive  trapping  of  2D  particle 
arrays  may  be  demonstrated. 

5.2  Structure  Design 

The  Slot-Suzuki-phase  photonic  crystal  lattice  we  propose  is  based  on  the  conventional 
Suzuki-phase  (SP)  lattice  shown  in  Figure  5-1  (a).  The  Suzuki-phase  lattice  is  obtained  by 
starting  with  a  triangular  lattice  of  holes  and  removing  selected  holes  to  generate  a  rectangular 
lattice  of  H\  cavities  [13,  14].  The  Suzuki-phase  lattice  has  different  periodicities  in  the  x  andy 
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directions,  equal  to  Sx  =  2a  and  Sy  =4^  a  respectively,  where  a  is  the  lattice  constant  of  the 
reference  triangular  lattice. 


d) 


H, 


•• 


Figure  5-1  (a)  Diagram  of  the  Suzuki-Phase  photonic  crystal  lattice,  (b)  Normalized 
transmission  spectra  calculated  by  the  transfer-matrix  method.  Red  line  for  x-polarization, 
blue  for  y-polarization.  (c)  //z-field  profile  (left)  and  (right)  of  ol  resonance,  d)  //z-field 

profile  (left)  and  (right)  of  e2  resonance. 


Figure  5-1  (b)  illustrates  the  normalized  transmission  spectrum  for  vertically  incident  light, 
calculated  by  the  Transfer  Matrix  Method  (TMM)  [15].  We  assume  a  high-refractive  index 
Suzuki-phase  lattice  with  relative  dielectric  constant  s  =  11.9,  hole  radius  r/a  =0.3  and  slab 
thickness  t/a  =  0.5  resting  on  an  oxide  substrate  (s  =  2.1),  and  immersed  in  a  fluid  with  dielectric 
constant  s  /  =  1.7.  Due  to  the  periodic  modulation  of  the  photonic  crystal  (PhC)  slab,  incident 
light  can  couple  to  the  GRM’s.  Guided  resonance  modes  are  strongly  confined  to  the  slab  and 
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appear  in  the  transmission  speetrum  as  Fano  line  shapes  superimposed  on  a  Fabry-Perot 
baekground  [16]  .We  foeus  in  partieular  on  F-point  modes  that  are  not  symmetry-forbidden  [17, 
18].  Figure  5-1  (b)  shows  the  four  GRM’s  we  are  interested  in.  We  label  the  four  modes  by  ol, 
o2,  el,  and  e2.  Figures  5-1  (e)  and  (d)  show  the  magnetie  field  eomponent  and  eleetric  field 
intensity  E  of  the  ol  and  e2  resonanee,  respeetively,  in  the  z  =  0  plane.  These  two,  dipole-like 
modes  are  eharaeteristie  of  the  set.  Modes  were  ealeulated  by  the  three-dimensional  (3D)  finite 
difference  time-domain  (FDTD)  method  [19]  and  normalized  to  the  maximum  value  in  the  z  =  0 
plane.  The  ol  mode  exhibits  odd  vector  symmetry  [20]  with  respect  to  the  x=  0  mirror  plane  {Hz 
is  even)  and  couples  to  a  plane  wave  with  electric  field  polarized  along  the  x  direction;  the  el 
mode  exhibits  even  vector  symmetry  with  respect  to  the  x  =  0  mirror  plane  {Hz  is  odd)  and 
couples  to  a  plane  wave  with  electric  field  polarized  in  the  y  direction,  as  has  been  shown 
experimentally  [21]. 

We  calculate  Q  factor  for  the  four  modes  in  Figure  5-1  (b)  using  3D  FDTD  calculations.  For 
a  real  photonic  device  with  finite  lateral  size,  the  total  Q  factor  depends  on  both  vertical  and 
lateral  losses.  However,  for  large  enough  structures  and  laser  excitation  spots,  it  is  reasonable  to 
consider  only  the  effect  of  vertical  loss  on  the  Q  factor  [22,  23].  We  model  the  PhC  structure  as 
infinitely  periodic  in  the  lateral  direction  by  imposing  periodic  boundary  conditions  along  the  x 
and  y  directions  of  the  unit  cell  in  the  FDTD  simulation.  The  corresponding  Q  factors  of  the  ol 
and  e2  modes  are  118  and  262,  respectively.  In  order  to  trap  a  particle  in  the  near  field  of  the 
PhC,  it  is  desirable  to  concentrate  the  optical  power  in  a  small  surface  area  and  to  have  a  large  Q 
factor  [5],  so  as  to  enhance  the  trapping  force  for  fixed  input  power.  However,  Figures  5-1  (c) 
and  (d)  show  that  the  electric  field  intensity  in  the  z  =  0  plane  of  both  the  x  -polarized  and  y- 
polarized  dipole  modes  is  relatively  spread  out  across  the  unit  cell.  In  order  to 
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increase  the  Q  factor  and  reduce  the  mode  area,  we  propose  a  Slot-Suzuki-phase  (SSP)  hybrid 
lattice,  shown  in  Figure  5-2  (a).  Rectangular  slots  with  cross-sectional  dimensions  of  Wx^Wy  are 
positioned  in  the  middle  of  eaeh  unit  cell  of  the  conventional  SP  lattice.  The  symbol  Wx  and  Wy 
represent  the  slot  length  along  the  x  and  y  axes,  respectively.  In  the  slot,  the  component  of 
electrie  field  normal  to  the  slot  boundary  is  enhanced  due  to  the  Maxwell  continuity  law  [24,  25]. 
If  the  smaller  dimension  of  the  reetangular  slot  is  along  the  y  direction  {wy  <  Wx),  the  electric 
field  of  the  y-polarized  el  and  e2  modes  will  be  enhanced.  To  obtain  enhancement  of  the  x  - 
polarized  ol  and  o2  modes,  the  slot  should  be  designed  with  Wx  <Wy.  In  the  following  discussion, 
we  focus  on  the  case  where  Wy  <  Wx. 


Figure  5-2  (a)  Diagram  of  the  Slot-Suzuki-phase  hybrid  lattice,  (b)  Normalized 
transmission  spectra.  Red  line  for  x  -polarization,  blue  fory-polarization.  (c)  Hz- field 
profile  of  se2  resonanee.  d)  field  profile  of  se2  resonanee. 
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In  Figure  5-2  (b),  we  plot  the  transmission  speetrum  of  the  SSP  lattiee  at  normal  ineidenee 
for  slot  dimensions  w/a  =  0.9  and  Wy/a  =  0.16.  The  other  struetural  and  material  parameters  are 
identieal  with  the  above-mentioned  SP  lattiee.  We  label  the  four  Fano  transmission  dips  as  sol, 
so2,  sel,  and  se2.  All  four  modes  exhibit  narrower  lineshapes  than  those  in  the  SP  lattiee  (Figure 
5-1  (b)).  Note  that  the  width  of  the  frequeney  window  is  the  same  for  Figure  5-1  (b)  and  Figure 
5-2  (b).  Thus,  higher  Q  faetors  ean  be  obtained  whieh  are  inversely  proportional  to  the  resonanee 
linewidths.  We  will  foeus  our  diseussion  of  field  loealization  on  the  se2  mode.  Figures  5-2  (e) 
and  (d)  plot  and  for  the  se2  mode  in  the  z  =  0  plane  of  the  SSP  lattiee.  The  eleetrie  field  of 
the  se2  mode  is  well  eonfined  around  the  eenter  of  the  slot  in  the  XY  plane.  In  faet,  the  modal 
volume  for  the  se2  mode  in  one  unit  eell  is  Fe//=  2.4x10  (Wn/)  ,  where  Xr  is  the  resonant 
wavelength  and  n/  is  the  refraetive  index  of  the  fluid.  The  mode  volume  per  unit  eell  is 
eomparable  to  that  in  PhC  mieroeavities  eontaining  slot  features  [25],  indieating  strong 
eonfinement  of  the  field  by  the  slot  within  eaeh  unit  eell  of  the  SSP  lattiee.  The  deeay  length  of 
the  field  intensity,  for  whieh  the  value  of  falls  to  1/e  of  the  value  at  the  slot  eenter,  are  0.28a  in 
the  X  -direetion,  0.09a  in  the  y-direetion,  and  0.3a  in  the  z  -direetion  (0.05a  above  the  PhC  slab). 
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Figure  5-3  (a)  Evolution  of  the  Q  and  wavelength  Xr  as  a  function  of  slot  width  Wy.  b) 
Evolution  of  the  Q  and  wavelength  Xr  as  a  function  of  slot  length  Wx.  Red  color  indicates  Q; 
blue  color  indicates  wavelength.  Dots  correspond  to  calculated  values;  lines  represent  a 

guide  for  the  eye. 


Adjusting  the  slot  dimensions  provides  a  flexible  strategy  for  tuning  the  Q  factor  and 
resonant  wavelength  Xr  of  the  GRM’s.  Varying  the  slot  width  (w;,)  can  dramatically  increase  the 
Q  factor  of  the  mode.  In  Eigure  5-3  (a),  we  plot  the  dependence  of  Q  on  Wy  for  fixed  Wx  =  0.9a 
for  the  se2  mode.  In  the  limit  where  the  slot  vanishes,  i.e.,  the  conventional  SP  lattice  with  only 
circular  holes,  the  Q  factor  for  the  e2  mode  is  262.  By  adding  a  narrow  slot  of  Wy  =  0.08a  in  the 
center  of  each  unit  cell,  the  Q  factor  is  increased  to  1750.  The  Q  factor  is  further  enhanced  by 
increasing  Wy  and  reaches  a  peak  value  of  around  123,000  at  Wy  =  0.16a.  This  Q  value  is 
comparable  to  values  obtained  in  square  PhC  lattices  for  coupled  GRM’s  [26].  If  the  slot  width 
Wy  continues  to  increase  beyond  0.16a,  the  Q  factor  decreases. 

The  slot  width  also  affects  the  resonant  wavelength  Xr.  The  blue  line  in  Eigure  5-3  (a)  shows 
that  Xr  linearly  decreases  from  1605  nm  for  Wy  =  0.08a  to  1510  nm  for  Wy  =  0.24a  with  a  slope  of 
approximately  -1.2.  The  graph  is  plotted  assuming  a  fixed  lattice  constant  a  of  515  nm. 
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Changing  the  slot  length  {w^  also  affeets  the  Q  faetor  and  resonant  wavelength  Xr.  We  fix  Wy 
=  0.16a  and  adjust  Wx.  As  indieated  by  Figure  5-3  (b),  the  Q  faetor  deereases  from  123,000  when 
Wx  deviates  from  0.9a,  but  remains  above  20,000  for  Wx  between  0.8a  and  1.15a.  The  blue  line 
shows  that  the  wavelength  Xr  remains  relatively  stable  with  a  shift  of  less  than  6  nm  in  the  same 
Wx  range.  This  relatively  low  variation  ean  be  explained  by  the  mode  distribution  shown  in 
Figure  5-2  (d).  The  eleetrie  field  is  tightly  eonfined  in  the  slot  with  a  deeay  length  in  the  x 
direetion  of  0.28a,  whieh  is  less  than  half  of  the  slot  length  (wx/2  =  0.45a).  Thus,  the  resonanee 
profile  is  only  lightly  infiueneed  by  Wx  deereasing  from  1.15a  to  0.8a.  The  ehoiee  of  Wx  can  be 
eonsidered  as  an  approach  for  fine-tuning  of  the  dipole  mode  Q  faetor  and  wavelength. 

5.3  Optical  Forces 

When  SSP  lattiees  with  different  slot  dimensions  are  eompared,  the  Q  factor  changes  by 
orders  of  magnitude,  while  the  mode  profile  remains  similar.  The  eomponent  Hz  has  a  dipole 
distribution  and  the  eleetrie  field  intensity  is  eoneentrated  in  the  slot.  The  SSP  lattiee  thus  offers 
the  flexibility  to  design  a  broad  range  of  Q  factors  and,  therefore,  trapping  forees.  The  main 
objeetive  of  this  seetion  is  to  predict  the  trapping  eapabilities  of  the  SSP  lattiee  with  lattiee 
eonstant  a  =  515  nm  on  a  dieleetrie  particle  {npoiy  =  1.60).  The  trapping  forees  exerted  on  the 
particle  are  eomputed  by  integrating  the  Maxwell  stress  tensor  (MST)  over  a  elosed  surfaee 
surrounding  the  partiele.  The  forees  are  numerieally  ealeulated  by  3D  FDTD  simulations.  We 
take  a  reetangular  solid  with  a  surfaee  several  mesh  points  away  from  the  nearest  edge  of  the 
partiele  as  the  integration  surfaee.  Due  to  the  high  Q  faetor  of  the  GRM,  it  is  eonvenient  to 
perform  the  foree  ealeulation  in  the  time  domain  rather  than  the  frequeney  domain.  We  exeite  the 
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mode  using  a  dipole  souree  inside  the  slab,  record  the  instantaneous  electromagnetic  fields  for 
several  optical  periods,  and  use  them  to  calculate  the  time-dependent  force.  We  then  time- 
average  the  force  and  normalize  it  to  the  power  P  coupled  to  the  se2  mode. 

We  calculate  the  optical  force  on  a  particle  of  radius  varying  from  25  nm  to  100  nm.  The 
particle  is  placed  right  above  the  center  of  the  slot  (at  position  x  =  0,  y  =  0),  with  its  bottom  edge 
45  nm  above  the  top  surface  of  the  slab.  Due  to  symmetry,  the  transverse  force  (Fxy)  on  the 
particle  vanishes.  Figure  5-4  (a)  shows  the  vertical  force  Fz/P  above  the  SSP  lattice  for  slot 
dimensions  Wx  =  0.9a  =  464  nm  and  Wy  between  0.08a  (41  nm)  and  0.24a  (124  nm).  The  force  is 
negative,  meaning  that  it  is  directed  towards  the  slab.  For  all  radii,  the  force  magnitude  increases 
to  a  peak  value  and  then  decreases  with  increasing  slot  width  (wy),  following  a  similar  trend  as 
the  Q  factor.  In  the  optimum  case,  a  slot  with  Wy=  0.16a  =  82  nm  enhances  the  force  by  two 
orders  of  magnitude  compared  to  a  slot  with  Wy  =  0.08a  =  41  nm.  The  optical  force  increases 
with  particle  radius.  For  a  particle  of  radius  25  nm,  the  maximum  force  magnitude  reaches  46  pN 
for  1  mW  power  per  unit  cell.  For  a  100-nm-radius  particle,  the  force  increases  to  46  pNmW~\ 
an  order  of  magnitude  enhancement.  In  Figure  5-4  (b)  we  plot  the  dependence  of  force  on 
particle  radius  for  a  particle  which  is  at  {x  =  0,  y  =  0)  and  has  its  bottom  edge  45nm  above  the  top 
surface  of  the  slab.  The  force  magnitude  increases  linearly  with  particle  radius  between  25  nm 
and  100  nm. 
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Figure  5-4  (a)  Optical  force  Fz  as  a  function  of  slot  width  Wy  for  four  different  particle  radii. 

The  slot  length  Wx  is  fixed  to  464  nm.  (b)  Optical  force  Fz  as  a  function  of  particle  radius. 
Dots  represent  the  calculated  values.  Lines  represent  a  guide  for  the  eye.  In  both  (a)  and  (b), 
the  particle  is  at  (x  =  0,  y  =  0)  and  has  its  bottom  edge  45  nm  above  the  top  surface  of  the 

slab. 


To  investigate  the  dependence  of  the  force  on  particle  position,  we  set  the  slot  dimensions  to 
464  nm  x  82  nm  (0.9a  x0.16a),  for  which  Q  achieves  its  optimal  value  of  123,000.  We  consider 
a  particle  of  radius  25  nm.  The  lower  edge  of  the  particle  is  45  nm  above  the  top  surface  of  the 
slab  (see  Figure  5-5  (a)).  We  calculate  the  vertical  force  (Fz)  and  transverse  force  (Fxy)  on  the 
particle  for  each  position  in  the  X¥ plane.  The  results  are  shown  in  Figures  5-5  (b)  and  (c).  The 
maximum  magnitude  of  Fz  is  approximately  46  pNmW~\  which  is  achieved  when  the  particle  is 
above  the  center  (x  =  0,  y  =  0).  For  all  positions  in  the  XT  plane,  the  vertical  force  is  directed 
towards  the  slab.  The  magnitude  of  the  in-plane  force  Fxy  is  shown  by  the  colormap  in  Figure  5-5 
(c),  and  the  force  direction  is  indicated  by  blue  arrows.  At  the  center  (x=  0,  y=  0),  the  in-plane 
force  vanishes.  The  strongest  in-plane  forces  point  toward  the  center.  The  maximum  in-plane 
force  for  the  particle  is  14  pNmW”^  and  is  weaker  than  the  maximum  vertical  force  of  46 
pNmW~\  The  ability  to  stably  trap  particles  is  thus  limited  by  the  in-plane  values.  Given  the 
spatial  map  of  the  in-plane  forces  Fxy  in  Figure  5-5  (c),  we  can  calculate  a  potential  map  in  the 
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XTplane.  We  calculate  the  potential  depth  A17by  integrating  the  in-plane  force  from  a  reference 
position  to  each  point  in  the  XY plane.  Here,  the  reference  position  is  taken  to  be  Sy/T).  The 
choice  of  the  reference  position  does  not  affect  the  relative  potential  depth.  The  potential  map  is 
shown  in  Figure  5-5  (d).  The  depth  of  the  trapping  potential  is  larger  than  377  KbT  for  an 
incident  power  of  1  mW  per  unit  cell,  where  KbT  is  the  thermal  energy  at  temperature  300  K.  A 
figure  of  merit  in  the  context  of  optical  trapping  is  the  stability  factor  S=IsXJ/KbT,  which  reaches 
377  mW~^  per  unit  cell  in  our  device.  In  order  to  achieve  a  stability  factor  equal  to  10,  generally 
considered  sufficient  to  achieve  stable  trapping  despite  the  presence  of  Brownian  motion  [27], 
the  power  needed  in  the  SSP  lattice  is  27  pW  per  unit  cell.  The  power  required  to  realize  stable 
trapping  in  the  SSP  lattice  is  reduced  by  about  40  times  compared  to  the  power  requirement  in 
the  square  lattice  proposed  in  our  previous  work  [12].  At  the  equilibrium  point,  the  trap  stiffness 
is  defined  as  -d  Fx/d  x  along  the  x  axis  and  -d  Fx/d  y  along  the  y  axis.  From  Figure  5-5  (c),  we 
estimate  trap  stiffnesses  of  0.12  and  0.28  pNnm~^mW~'  for  the  25  nm  particle  along  the  x  andy 
axes,  respectively,  corresponding  to  a  radial  trap  stiffness  of  0.08  pNnm~'mW~^  in  the  XY  plane. 
Greater  trap  stiffness  results  in  lower  uncertainty  in  the  trapping  position.  The  stability  factor  and 
the  radial  trap  stiffness  we  obtained  exceed  the  values  reported  for  dielectric  trapping  structures 
in  previous  works  [5]. 


92 


Ill 


X  (nm)  X  (nm) 


Figure  5-5  (a)  Diagram  of  the  Slot  Suzuki-phase  (SSP)  lattice  with  a  particle  above,  (b) 
Vertical  force  as  a  function  of  particle  position  in  the  XY  plane,  (c)  In-plane  force  Fxy  as 
a  function  of  position  in  the  XY  plane.  The  force  magnitude  is  indicated  by  the  colormap, 
and  the  force  direction  is  shown  by  the  blue  arrows,  (d)  Potential  map  in  the  XY  plane.  To 
obtain  the  results  shown  in  (b),  (c),  and  (d),  we  assume  that  the  particle  has  a  radius  of  25 
nm  and  is  placed  such  that  its  bottom  edge  is  45  nm  above  the  top  surface  of  the  slab. 


For  particle  diameters  less  than  the  slot  width  Wy,  particles  can  be  trapped  inside  the  slot.  The 
particle  is  physically  confined  by  the  slot  in  the  y  direction.  We  calculate  the  force  as  a  function 
of  position  in  the  XZ  plane  for  y  =  0.  Results  are  shown  in  Figure  5-6  (a).  The  force  Fy,  which  is 
normal  to  the  XZ  plane,  vanishes  due  to  symmetry,  and  thus,  we  only  plot  the  in-plane  force  F^z. 
The  forces  point  to  the  center  of  the  slot  at  (x  =  0,  z  =  0)  and  have  a  maximum  magnitude  of  153 
pNmW~\  The  corresponding  potential  map  shown  in  Figure  5-6  (b)  indicates  that  a  strong 
trapping  potential  depth  of  more  than  3500  KgT  per  milliwatt  per  unit  cell  is  achieved  within  the 
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slot.  For  a  particle  within  the  slot,  a  power  as  low  as  3  pW  per  unit  cell  is  required  for  stable 
optical  trapping.  The  radial  trapping  stiffness  in  the  XZ plane  is  -0.34  pNnm~'mW~*.  Within  the 
slot,  the  stability  and  trapping  stiffness  are  higher  than  outside  the  slot  and  comparable  to 
reported  values  for  particle  trapping  in  other  slot  structures  [10,  11]. 

The  presence  of  a  particle  in  the  slot  can  affect  the  resonance  wavelength.  Our  simulations 
show  that  when  every  unit  eell  of  the  structure  eontains  one  partiele  of  radius  25  nm  in  the  eenter 
of  the  slot,  the  resonance  wavelength  is  shifted  by  0.3  nm.  Nevertheless,  in  a  real  experiment,  the 
particles  are  likely  to  be  trapped  one  by  one,  with  the  number  inereasing  gradually  over  time.  In 
this  case,  any  given  particle  will  create  a  negligible  perturbation  on  the  mode.  As  larger  numbers 
of  particles  are  trapped,  a  gradual  shift  of  the  resonanee  wavelength  will  take  plaee,  requiring  a 
slow  adjustment  of  the  excitation  laser  (AXr  <0.3  nm). 


a)  FJP  (pN/mW)  b)  AU/KgT/P  (mW'') 


Figure  5-6  (a)  In-plane  optieal  force  F^z  for  a  partiele  with  radius  of  25  nm.  (b)  Optieal 

potential  map  for  the  particle. 
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5.4  Summary 

We  have  proposed  a  new  PhC  lattiee  for  optieal  trapping  of  two-dimensional  arrays  of 
nanoparticles.  Our  structure  is  created  by  using  the  2D  Suzuki-phase  PhC  lattice  as  a  base  and 
introducing  a  slot  into  each  unit  cell  to  localize  the  electromagnetic  field.  Optimizing  the  slot 
dimensions  increases  the  Q  factor  of  the  resonance  by  orders  of  magnitude.  Optical  power  values 
as  low  as  27  pW  per  unit  cell  are  predicted  for  trapping  of  25  nm  radius  beads,  a  reduction  of 
power  by  about  40  times  relative  to  our  previous  work.  Once  the  particle  is  on  the  slot,  optical 
power  of  3  pW  per  unit  cell  is  required  for  the  stable  optical  trapping.  Our  Slot-Suzuki-Phase 
lattice  is  a  promising  candidate  for  carrying  out  light-assisted  templated  self  assembly  processes. 
The  low  power  requirements  for  trapping  suggest  to  use  of  active  materials  for  this  purpose.  Slot 
photonic  crystal  microcavity  lasers  have  been  demonstrated  with  output  optical  power  as  high  as 
150  pW,  and  some  evidence  suggests  possible  optical  trapping  effects  in  such  structures  [28,  29]. 
The  Slot- Suzuki-phase  structure  we  propose  here  provides  a  way  of  effectively  combining 
multiple  slot  PhC  microcavities  into  a  high-Q  structure  with  extended  area.  Further  improvement 
of  the  design  can  be  achieved  by  band  engineering  techniques,  as  well  as  by  combining  the 
photonic  crystal  with  a  bottom  Bragg  reflector  [30,  31].  In  such  a  device,  we  expect  that  the  laser 
may  self  adapt,  adjusting  its  own  lasing  wavelength  in  response  to  the  resonance  shift  induced  by 
trapped  particles. 
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Chapter  6:  Characterization  of  Slot-graphite  Photonic  Crystal 
Lattice 

6.1  Background 

Photonic  crystal  structures  have  distinguished  points  with  zero  group  veloeity  along  the 
Brillouin  boundary.  The  zero-group-velocity  modes  have  been  applied  in  many  interesting  fields, 
including  laser  emission  [1-3],  quantum  optics  [4],  and  sensing  [5].  One  interesting  point  with 
zero  group  veloeity  along  the  Brillouin  boundary  is  /"-point.  The  /"-point  modes  are  above  the 
light  line.  The  modes  can  be  coupled  or  uncoupled  to  an  external  plane  wave  depending  on  the 
mode  symmetry  [6].  When  light  is  normally  incident  upon  the  periodically  patterned  slab,  it  can 
excite  /"-point  modes  or  as  so  called,  guided  resonance  modes,  whieh  propagate  in  the  plane  of 
the  slab  while  leaking  partially  to  the  surrounding  air  [7]. 

In  the  last  Section,  we  have  discussed  utilizing  /"-point  modes  in  a  Slot-Suzuki-phase 
photonic  crystal  lattice  for  light-assisted,  template  assembly  (LATS)  of  nanoparticles.  We  have 
pointed  out  that,  for  fixed  input  power,  the  optieal  forees  can  be  enhanced  by  using  a  high- 
quality-factor  and  well-confined  guided  resonance  modes.  In  other  words,  the  assembly 
capability  is  dependent  on  optieal  quality  factors.  In  this  Section,  we  experimentally 
demonstrated  the  high-Q  guided  resonance  modes  in  Slot-graphite  lattiees.  The  lattice  is  ereated 
by  using  regular  graphite  lattice  [8]  as  the  base  and  placing  a  rectangular  slot  into  the  center  of 
each  unit  cell.  Quality  factors  over  4x10*  are  predieted  by  theory.  The  Slot-graphite  lattices  have 
been  fabricated  on  silieon-on-insulator  platforms  and  optically  characterized  by  transmission 
spectroseopy.  We  have  observed  the  dependenee  of  resonant  wavelengths  and  Q  factors  on  the 
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slot  dimension.  Note  that  by  fine  eontrol  of  slot  geometry,  we  ean  tune  the  Q  faetors  by  orders  of 
magnitude.  The  modes  have  a  strong  field  eonfinement  within  the  slots,  proposing  potential 
applieations  in  LATS. 

6.2  Structure  Design 

The  Slot-graphite  photonie  erystal  lattiee  is  ereated  by  using  regular  graphite  lattiee  [8]  as  the 
base  (See  Figure  6-1  (a)).  In  the  eenter  of  eaeh  unit  eell,  we  plaee  a  reetangular  slot  oriented 
towards  the  middle  point  of  two  elosest  holes  (see  Figure  6-1  (b)).  This  novel  strueture  uses  the 
same  Bravais  lattiee  as  the  graphite  lattiee.  Nevertheless  the  slot  reduees  the  lattiee  symmetry 
from  Cfc  to  C2v  [9],  and  the  irredueible  Brillouin  zone  ehanges  from  1/12  of  the  first  Brillouin 
zone  (See  Figure  6-2  (a)  inset)  to  1/4  (Figure  6-2  (b)  inset).  The  irredueible  Brillouin  zone  of 
graphite  lattiee  is  marked  by  T-M-K-  F  in  Figure  6-2  (a)  inset,  while  slot-graphite  lattiee  is  F-M- 
K-Ml-Kl-  F  in  Figure  6-2  (b)  inset. 

Figure  6-2  (a)  shows  the  3D  photonie  band  strueture  of  a  graphite  lattiee.  The  bands  are 
ealeulated  using  guide  mode  expansion  method  [10].  We  assume  a  lattiee  eonstant  a=800  nm,  a 
hole  radius  r=0.157a  and  slab  thickness  of  250  nm  laying  on  top  of  silica  substrate.  The  device  is 
immersed  in  a  material  of  dielectric  constant  £6=1.77. 
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Figure  6-1  (a)  Diagram  of  the  graphite  photonic  crystal  lattice,  (b)  Diagram  of  the 
Slot-graphite  photonic  crystal  lattice.  The  arrows  represent  the  primitive  vectors 

of  the  Bravais  lattice. 


In  Figure  6-2  (a),  we  show  the  band  structure  of  the  graphite  lattice.  At  F-points  as 
highlighted  by  dashed  circles,  the  fourth  and  fifth  bands  are  degenerate.  The  corresponding  mode 
profiles  are  shown  by  Figure  6-2  (c).  The  /A  components  of  two  degenerate  modes  have  dipole- 
like  distribution  with  different  orientation.  And  we  call  them  dipole  modes.  The  electric  field  of 
the  fourth  (fifth)  modes  is  x  (y)  polarized.  Note  that  the  bands  are  the  same  for  the  two  T  points 
in  the  M-T-K  boundary  and  Ml-T-Kl  boundary.  The  reason  has  been  given  in  the  last  paragraph; 
in  the  graphite  lattice,  the  M-T-K  Brillouin  zone  is  irreducible  and  contains  all  band  information. 

In  Figure  6-2  (b),  we  show  the  band  structure  of  our  slot-graphite  lattice.  The  slot  dimension 

is  Wx=0.615a  and  Wy=0.09a.  The  other  geometries  are  identical  with  the  graphite  lattice.  The 

degeneration  of  the  fourth  and  fifth  band  of  the  graphite  lattice  at  the  T  -point  is  removed  due  to 

the  breaking  symmetry  created  by  the  slot.  By  adding  slot  in  each  unit  cell,  a  sharp  discontinuity 

of  the  dielectric  constant  is  created.  Due  to  the  Maxwell  continuity  condition,  the  normal 

components  of  eE  are  conserved  along  the  slot  boundary  [11,  12],  so  that  the  Ey  field  of  the  fifth 

band  is  enhanced  within  the  slot  as  shown  by  the  upper  part  of  Figure  6-2  (d).  The  light 
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confinement  within  slot  is  not  obvious  for  the  fourth  band,  sinee  this  band  is  polarized  along  x 
direction.  In  the  slot  device,  we  are  going  to  refer  to  the  fifth  band  as  slot  mode  whereas  the 
fourth  band  as  dipole  mode.  Note  that  the  slot-mode  band  around  the  F-point  is  fiat,  which  can 
be  very  important  for  specific  application  [13]. 


Figure  6-2  (a)  Photonic  band  structure  of  the  Graphite  photonic  crystal  lattice,  (b)  Photonie 
band  structure  of  the  Slot-graphite  photonic  crystal  lattice.  The  highlighted  region 
represents  the  targeted  bands,  (c)  //z-field  profile  (left)  and  (right)  profile  of  the  fifth 
band  (down)  and  fourth  band  (up)  of  the  graphite  photonie  crystal  lattice,  (d)  //z-field 
profile  (left)  and  E  (right)  profile  of  the  fourth  band  (down)  and  fifth  band  (up)  of  the 

Slot-graphite  photonic  crystal  lattice. 
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We  calculated  the  temporal  confinement  of  the  slot  mode  at  the  F-point.  Figure  6-3  (a)  shows 
the  dependence  of  the  Q  factor  (red  curve)  and  wavelength  (black  line)  on  the  slot  height  Wy, 
with  fixed  width  Wx  =0.675a.  The  simulation  is  performed  by  solving  the  Maxwell  equations 
using  a  3D  finite  difference  time-domain  (FDTD)  code,  Lumerical.  Interestingly,  the  Q  factor 
increases  from  1970  for  Wy  =0.06a  until  reaching  4.3x10^  for  Wy=0.0925a.  As  we  further  increase 
the  slot  width,  the  Q  factor  decreases  to  3100  for  Wy  =0. 1225a.  On  the  other  hand,  the  mode 
wavelength  decreases  as  Wy  increases,  following  a  linear  trend. 

Figure  6-3  (b)  shows  the  evolution  of  the  Q  factor  (red  curve)  and  wavelength  (black  line) 
with  slot  width  Wx,  for  fixed  slot  height  Wy=0.0925a.  The  Q  factor  decreases  as  Wx  is  increasing. 
Note  that  the  Q  factor  is  less  sensitive  to  w^than  Wy.  For  example,  a  change  of  Wx  of  20  nm  with 
respect  to  the  peak  preserves  a  Q-factor  over  85000  (20%  of  the  peak  value).  While  changing  20 
nm  of  Wy  from  the  peak  value  decreases  Q  from  4.3x10^  to  4300  (1%  of  the  initial  value).  The 
evolution  of  the  wavelength  with  Wx  has  two  regions  clearly  differentiated;  one  with  Wx  <0.675a 
and  another  one  with  Wx  >0.675a.  In  the  first  region,  the  wavelength  change  is  the  ~  1  nm  for  a 
76  nm  decrement  in  the  slot  length.  In  the  second  region,  there  is  a  change  of  ~  9  nm  for  a  78  nm 
increment  in  the  slot  length. 

After  evaluating  the  temporal  confinement,  we  evaluate  the  spatial  confinement  by  calculating 
modal  volume  per  unit  cell.  For  the  slot  mode  with  the  Q  factor  4.3x10^,  we  obtain  a  value  of 
V^j^^=0.036{X  /  2n  i)^ .  This  value  of  the  spatial  confinement  compares  well  with  the  previous 
results  obtained  in  PhC  microcavities  [14,  15]. 
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Figure  6-3  (a):  Evolution  of  the  Q  and  wavelength  as  function  of  Wy.  (b):  Evolution  of  the 
Q  and  wavelength  as  function  of  Wx.  Red  color  for  Q.  Black  for  wavelength.  Dots; 
Corresponds  to  the  calculate  values.  Eines  represents  a  guide  for  the  eye. 


6.3  Fabrication  and  Optical  Characterization 

We  fabricate  a  group  of  slot-graphite  photonic  crystal  devices  in  order  to  characterize  the 
evolution  of  the  slot  modes  with  varying  slot  height  Wy.  The  devices  are  fabricated  with  fixed 
lattice  constant  a=820  nm,  radius  of  air  holes  r/a=0.155,  and  slab  thickness  t=250  nm.  We  adjust 
the  device  parameters  so  that  we  have  a  relatively  low  Q<10^  for  easier  optical  characterization. 
We  expect  to  see  the  resonant  wavelengths  of  slot  modes  fall  in  our  tunable  laser  range,  and  the 
Q  factor  will  show  an  increasing-then-decreasing  trend  with  increasing  Wy.  The  computed  Q 
factor  is  not  far  away  of  the  highest  values  (Q-factor  over  9800)  we  measured  previously  in  the 
same  set-up  for  uncoupled  guided  resonant  modes  [16]. 
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Figure  6-4  SEM  image  of  one  fabricated  Slot-graphite  device. 


The  devices  have  been  fabricated  using  e-beam  lithography  and  inductively-coupled  plasma 
reactive  ion  etching  (ICP-RIE).  We  use  a  250  nm  thick  silicon  layer  on  top  of  a  3  pm  silica  layer 
(SOITEC).  The  sample  was  spin-coated  with  PMMA-A4  950K.  Elsing  Raithl50  e-beam  system, 
a  50  pm  diameter  Slot-graphite  pattern  has  been  exposed  with  an  acceleration  voltage  of  30KV. 
The  pattern  is  transferred  from  the  resist  to  the  silicon  layer  by  an  ICP-RIE  machine  using  a  gas 
mixture  SE6/C4E8.  Einally  the  remaining  resist  was  removed  using  O2  plasma  etching  and 
acetone.  Figure  6-4  shows  a  Scanning-electron  beam  (SEM)  image  of  one  of  the  fabricated 
devices.  We  analyze  the  SEM  images  and  extract  geometrical  parameters  [17].  The  slot  width  Wx 
is  fixed  to  550  nm,  and  Wy  is  increased  gradually  from  58  nm  to  93  nm  at  3.3  nm  by  each  step. 
Note  that  the  SEM  inspection  can  introduce  large  systematic  errors  as  much  as  5-10%  of  the 
absolute  dimension  [18]. 
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We  characterize  the  transmission  spectra  of  our  fabricated  samples  using  cross  polarization 
configuration.  The  sample  was  immersed  in  deionized  water.  We  measure  the  wavelength  and  Q 
factor  of  the  slot  modes  using  our  experimental  configuration  as  shown  by  Figure  6-5  (a). 
Incident  light  of  a  tunable  laser  (TL)  in  near-infrared  range  between  1500  nm  and  1620  nm  is 
collimated  from  fiber  to  free  space  by  a  lens  (LS).  The  light  is  focused  onto  the  sample  on 
translation  stage  (TS)  using  a  low  numerical  aperture  microscope  objective  (MO).  The  light 
transmitted  through  the  sample  is  collected  using  another  microscope  objective  and  focused  onto 
a  photodetector  (PD).  Note  that  two  crossed  polarizers  (PC)  are  positioned  before  and  after  the 
microscope  objectives,  which  can  cancel  out  the  Fabry-Perot  background  in  the  guided- 
resonance  transmission.  With  the  cross-polarization  setup,  we  measured  slot  modes  as  peaks  in 
the  transmission  spectrum  rather  than  dips  on  the  Fabry-Perot  background  as  shown  by  Figure  6- 
5  (b)  in  the  previous  section.  Figure  6-5  (b)  shows  a  transmission  spectrum  measured  in  a  slot- 
graphite  device.  The  slot  mode  corresponding  to  the  fifth  band  and  the  dipole  mode 
corresponding  to  the  fourth  band  are  clearly  visible. 
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Figure  6-5  (a)  Diagram  of  experimental  setup  used  to  characterize  the  optical  properties  of 
the  guided  resonant  device.  Two  cross  polarizers  (PC)  are  applied  before  and  after  the 
device  to  cancel  out  Fabry-Perot  noise.  TL;  tunable  laser;  MO;  microscope  objective;  LS; 

lens;  TS;  translation  stage;  PD;  photodiode,  (b)  One  measured  spectrum.  The  right 
transmission  peak  correspond  to  the  slot  mode,  and  the  left  peak  shows  the  dipole  mode 

with  lower  Q-factor. 


The  tested  Q-factor  and  wavelength  of  the  slot  modes  are  plotted  by  black  spheres  in  Figure 
6-6.  We  did  not  plot  the  evolution  of  the  dipole  modes  for  simplicity.  We  compared  the 
experimental  data  with  the  simulation  shown  by  gray  dashed  curves  in  Figure  6-6.  In  FDTD 
simulation  we  take  into  account  water  absorption  in  the  near-infrared  range. 

As  shown  by  Figure  6-6  (a),  within  our  expectation,  the  measured  Q-factor  first  increased 
and  then  decreased  with  increasing  Wy.  The  measured  data  are  missing  for  devices  with  Wy 
between  65  nm  and  71  nm.  The  possible  reason  is  that  the  slot  modes  in  these  devices  have  high 
Q  factors.  For  high  Q  modes,  it  is  more  challenging  for  external  light  to  couple  into  the  mode. 
We  see  the  intensity  of  the  transmission  peaks  for  guided  resonance  comparable  with  the 
background  noises.  Due  to  the  low  signal-to-noise  ratio,  we  did  not  obtain  corresponding  Q 
factors  and  wavelengths.  The  highest  Q  we  measured  is  around  5300.  Note  that  compared  with 
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simulation  (gray  dashed  curve),  the  tested  Q  is  lower  by  a  factor  of  2.  This  offset  may  be  due  to 
fabrication  imperfection  and  surface  roughness. 

Figure  6-6  (b)  shows  the  tested  wavelength  evolution  (black  spheres)  and  the  corresponding 
simulation  (gray  dashed  curve).  The  measured  wavelength  monotonically  decreases  from  1613 
nm  to  1560  nm  with  Wy  increasing  from  58  nm  to  93  nm.  On  average,  the  wavelength  decreases 
by  1 .5  nm  if  Wy  is  increased  by  1  nm. 
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Figure  6-6  (a):  Evolution  of  the  Q  factor  as  function  of  Wy.  (b):  Evolution  of  the 
wavelength  as  function  Wy.  Black  spheres  for  experimental  data.  Gray  dashed  curves  for 
simulation.  Blue  stars  for  experimental  data  when  Wy  is  increased  by  10%  by  considering 

SEM  system  errors. 


The  overall  wavelength  and  Q-factor  trends  follow  our  simulation  curves.  However,  the 
tested  wavelength  is  smaller  than  simulated  value  by  about  15  nm.  One  reason  is  that  the 
wavelength  is  very  sensitive  to  slot  height  Wy,  while  the  Wy  value  we  obtained  through  the  SEM 
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images  can  have  offset  as  much  as  10%  of  the  feature  dimension  [18].  We  use  the  blue  stars  in 
Figure  6-6  to  show  the  experimental  data  if  the  real  slot  height  Wy  is  increased  by  10%  of  the 
SEM  size.  In  this  case,  the  tested  and  simulated  wavelengths  are  in  good  agreement. 

6.4  Summary 

A  new  photonic  crystal  lattice  has  been  proposed  as  slot-graphite  lattice.  The  lattice  is 
created  by  using  the  graphite  lattice  as  the  base  and  placing  a  slot  in  the  center  of  unit  cell. 
Quality  factors  over  4x10^  are  predicted  for  guided  resonance  modes  with  light  well  confined 
within  slots.  The  structures  have  been  fabricated  on  SOI  wafer  and  immersed  in  D.I.  water  for 
measurement.  The  modes  have  been  characterized  by  cross-polarization  transmission 
spectroscopy.  The  evolution  of  the  Q  factors  and  wavelengths  follows  the  trend  predicted  by  our 
simulations.  Quality  factors  up  to  5300  have  been  measured.  Due  to  highly  tunable  Q  factor  and 
strong  light  confinement,  the  lattice  is  promising  for  lower-power  light  assisted,  templated 
assembly. 
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Chapter  7:  Radiation  Loss  in  Photonic  Crystal  CROW 

7.1  Background 

Coupled  resonator  optieal  waveguides  (CROWs)  [1],  waveguides  in  whieh  light  propagates 
by  “hopping”  between  loealized  resonator  modes,  have  been  widely  studied  for  applieation  to 
on-ehip  optieal  delays  [2],  CROWs  provide  both  slow  group  veloeity  and  low  dispersion  [3], 
essential  features  for  distortion-free  signal  delay.  In  dynamieally-tuned  CROWs,  light  pulses 
may  be  stopped,  stored,  and  released  on  demand  [4,  5].  CROWs  based  on  various  resonator  types 
have  been  investigated,  ineluding  mieroring  resonators  [6-8]  and  photonie-erystal  mieroeavities 
[9-15].  We  foeus  here  on  photonie-erystal  mieroeavity  CROWs,  whieh  provide  both  eompaet 
resonator  size  and  high  spatial  eonfinement;  the  latter  is  useful,  for  example,  for  aehieving 
nonlinear  effeets  at  low  power  levels  [16]. 

Minimization  of  propagation  loss  is  an  important  aspeet  of  CROW  design.  Various 
approaehes  have  been  used  to  design  low-loss  CROWs,  ineluding  eoupled-mode  theory  [17]  and 
the  finite-element  method  [18].  CROWs  formed  by  mieroeavities  in  photonie-erystal  slabs  are 
intrinsieally  lossy,  beeause  the  CROW  mode  lies  above  the  light  line  of  air,  the  mode  leaks  light 
vertieally  as  it  propagates  [19].  Moreover,  if  the  loss  is  uneven  aeross  the  CROW  band  for 
different  waveveetors,  the  effeetive  transmission  bandwidth  is  redueed,  and  the  transmitted 
signal  is  distorted.  The  intrinsie  loss  provides  a  lower  bound  on  the  total  loss,  whieh  may  also 
inelude  eontributions  due  to  fabrieation  imperfeetions  and  seattering. 

Previous  theoretieal  work  has  shown  that  intrinsie  CROW  loss  ean  depend  strongly  on 
frequeney  and  be  up  to  an  order  of  magnitude  larger  or  smaller  than  that  of  an  individual 
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resonator  [20,  21],  At  the  same  time,  experiments  have  demonstrated  particular  CROW  designs 
with  low  loss  comparable  to  that  of  an  individual  resonator  [22],  Given  these  considerations,  it  is 
important  to  understand  theoretically  how  the  CROW  and  resonator  losses  are  related,  as  well  as 
how  to  achieve  low,  flat  loss  across  the  entire  CROW  bandwidth.  In  this  paper,  we  use  a  model 
based  on  the  tight-binding  approximation  (TBA)  to  answer  these  questions. 

TBA  approaches  are  advantageous,  as  they  allow  the  calculation  of  the  frequency-dependent 
loss  of  a  CROW  from  a  single,  full-field  simulation  of  an  individual  cavity.  Previously,  Fussell 
and  Dignam  have  used  a  TBA  model  based  on  the  complex  dispersion  relation  of  the  CROW  to 
calculate  losses  [21].  Another  approach  is  to  use  the  TBA  to  relate  the  electromagnetic  fields  of 
the  CROW  to  those  of  the  isolated  cavity  [15].  In  this  case,  the  spatial  power  spectrum  of  the 
CROW  fields  may  be  obtained  from  that  of  the  isolated  cavity  by  taking  discrete  samples  in 
Fourier  space,  as  experimentally  demonstrated  in  [15].  Reference  [15]  has  suggested  that  the 
losses  may  then  be  calculated  from  the  power  spectrum  using  a  light-cone  picture  [23].  In  this 
paper,  we  explicitly  demonstrate  the  validity  of  this  approach  by  comparing  the  predictions  of 
the  TBA  light-cone  model  with  direct  calculations.  We  then  show  that  within  a  particular  family 
of  photonic -crystal  cavities,  a  decrease  in  cavity  loss  also  results  in  a  decrease  in  CROW  loss. 
While  the  CROW  loss  generally  varies  across  the  CROW  band,  we  show  that  the  TBA  predicts 
particular  cavity  separations  for  which  the  loss  is  flat. 

We  expect  our  results  to  be  useful  for  a  range  of  applications.  In  the  context  of  light- assisted 
self  assembly  above  photonic-crystal  slabs  [24,  25],  our  method  may  facilitate  the  design  of 
large-area,  coupled-cavity  modes  with  high  quality  factor  that  nevertheless  couple  to  normally- 
incident  radiation.  For  slow-  and  stopped-light  applications,  the  fidelity  of  signal  transmission  at 
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the  end  of  a  delay  line  depends  on  aehieving  low,  flat  loss  aeross  the  signal  bandwidth.  In 
eoupled-eavity  lasers,  the  ability  to  tailor  the  Q{k)  loss  profde  eould  provide  a  useful  method  to 
control  the  effects  of  mode  competition  and  tune  the  directional  output  of  the  laser. 

7.2  Theoretical  Formulation  in  the  Tight-binding  Approximation 

For  concreteness,  we  first  consider  the  CROW  structure  shown  in  Figure  7-1  (a).  The 
photonic  crystal  slab  is  a  triangular  lattice  of  air  holes  with  lattice  constant  a  and  hole  radius  in  a 
suspended,  dielectric  slab  with  height  and  refractive  index  1.55pm).  The  CROW  waveguide 
consists  of  periodically  spaced  defects.  Each  defect  is  a  micro-cavity  with  two  adjacent  missing 
holes,  known  as  an  L2  cavity.  The  center-to-center  spacing  of  adjacent  defects  is  5a.  The 
separation  width  S,  the  width  of  the  region  between  defects  that  does  not  contain  missing  holes, 
is  3a.  We  denote  the  CROW  as  LmSn,  with  m  being  the  number  of  missing  holes  of  the 
constituent  cavity  and  n  the  number  of  separation  holes.  We  calculate  the  band  structure  by  the 
guided-mode  expansion  method  [26].  Figure  7-1  (b)  shows  the  dispersion  relations  for  L2S2, 
L2S3,  and  L2S4  CROWs,  respectively.  The  band  of  each  CROW  mode  displays  a  cosinusoidal 
shape  [1,  3],  with  the  center  frequency  close  to  the  fundamental  mode  of  the  isolated  L2  cavity, 
0.26992(27ic/a).  The  bandwidth  of  a  CROW  mode  is  proportional  to  the  coupling  strength 
between  adjacent  defects.  Figure  7-1  (b)  shows  that  the  coupling  strength  decreases  as  the 
separation  is  increased. 
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Figure  7-1  (a)  Schematic  diagram  of  CROW  structure  L2S3.  (b)  Dispersion  curve  of 
CROW  mode  in  L2S2  (red),  L2S3  (blue)  and  L2S4  (magenta).  Dashed  line  shows  the 
fundamental  mode  of  the  isolated  L2  cavity. 


For  all  three  CROWs,  the  CROW  band  lies  above  the  light  line  of  air  and  is  thus  intrinsically 
lossy.  The  light  line  of  air  is  given  by  (O-ck,  where  k  is  the  in-plane  CROW  wavevector.  At  the 
edge  of  the  Brillouin  zone,  k  =  k^^=nlD  and  the  light-line  frequency  (u=0.125(27ic/a), 

0.1(27ic/a),  or  0.083(27ic/a)  for  D  =  Aa,  D  =  5a  ox  D  =  6a,  respectively.  In  all  three  cases,  the 
upper  edge  of  the  light  line  lies  well  below  the  frequency  range  shown  in  Figure  7- 1(b). 

We  define  a  waveguide  quality  factor  Q{k)  to  describe  the  leakage  of  the  CROW  mode; 

Q{k)^(0{k)T,^  (7-1) 

where  co{k)  is  the  dispersion  relation  of  the  CROW  mode,  and  ii/e  is  the  time  in  which  the  power 
in  the  waveguide  mode  decays  by  1/e.  Q{k)  can  be  converted  to  a  1/e  power-decay  length  in  the 
waveguide  as  Ly^{k)  =  Vg{k)Q(k)/  (0{k) ,  where  v^(A:)  is  the  group  velocity.  Alternately,  the  loss 
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in  dB/length  is  equal  to  -4.34dB/Zi/e.  While  the  loss  per  unit  length  diverges  at  A:  =  0  and 
0.5(27i/Z))  where  Vg=0,  Q{k)  is  well  defined  over  the  entire  Brillouin  zone. 


7.2.1  Light  Line  Picture 

For  isolated  mieroeavities,  it  has  been  shown  [23]  that  the  eavity  Q  ean  be  related  to  the 
Fourier  transform  of  the  eleetromagnetie  fields  above  the  surfaee  of  the  slab.  A  similar  relation  is 
true  for  the  waveguide  quality  faetor  Q{k). 

We  ean  write: 


Q{k)  =  (o{k) 


{Ujk)) 

{P(k)) 


(7-2) 


where  (U (k)'j  is  the  time-averaged  eleetromagnetie  field  energy: 


(U(k))  =  \\l\dr(e 


CROW 


{x,y,z)  Ei^{x,y,z) 


+  Po 


H, 


(7-3) 


where  £crow  is  the  dieleetrie  eonstant  of  the  CROW,  and  (x,  y,  z)  and  EI^  {x,  y,  z)  are 

the  eleetrie  and  magnetie  fields  of  the  CROW  mode  at  wave  veetor  k.  {P{k)')  is  the  time- 
averaged  emitted  power, 

{P{k))  =  (^  6/5  ~  Re  {e^{x,  y,  z)  x  Hl(x,  y,  z) ),  (7-4) 

where  the  integral  is  over  a  elosed  surfaee  that  eneloses  the  photonie  erystal  slab.  For  an 
infinitely  long  CROW,  both  (y{k)'j  and  {P{k)'j  are  infinite;  however,  the  ratio  (U(k)'j  /  (^P{k)'^ 
is  finite  and  equal  to  the  value  of  the  ratio  for  a  single  unit  eell. 
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{P(k)'j  can  be  rewritten  in  terms  of  the  two-dimensional  (2D)  spatial  power  spectra  of  the 
free-space  fields,  on  any  surfaee  above  the  strueture  [23]  ; 


{P{k))  = 


1  IKff 

2]l  2n  2k 


(7-5) 


where  -Wx  ^kii^x’^y)  Hi^.{q^,qy)  are  the  2D  Fourier  transforms  of  the 

fth  component  of  Ej^{x,y,z)  and  Hj^{x,y,z)  ,  and  (ie  x,y) : 


^k,.  i^x  ’^y)  =  \\  ^  y)  e 

H k,i  (qx  ’  )  =  Jj dxdyH^  .  ( :v;  y)  e (7-6) 
The  Fourier-space  integral  includes  all  q^^  within  the  light  eone  <  m/c).  Only  these  so-called 
“leaky  components”  contribute  to  radiation  loss.  We  will  refer  to  the  quantities  \E,^ Xqx^qy)\  ’ 

\H^  ^{q^,q^)^  appearing  in  Equation  (7-5)  as  the  E^,  Ey,  H^, 

and  Ely  power  spectra,  respectively. 

Due  to  the  periodicity  of  the  CROW  waveguide  in  the  x-direction,  Fourier  transforms  of  the 
field  components  are  only  nonzero  for  diserete  values  of  q^ .  Using  Bloeh’s  theorem,  the  CROW 

mode  ean  be  written  as  Ej^  (x,  y,  z)  =  (x,  y,  z) ,  where  (x  -I-  D,  y,  z)  =  U^.  (x,  y,  z) .  An 

analogous  statement  is  true  for  Hi^{x,y,z) .  It  follows  that  the  Fourier  transforms 
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i ’^y)-  ^kii^x’^ly)-^  except  when  q^=k  +  Ikh  / D  .  The  Fourier  transforms  are  nonzero 
over  a  continuous  range  of  values. 


7.2,2  Tight-binding  Approximation 

The  tight-binding  approximation  can  be  used  to  write  the  CROW  quality  factor  in  terms  of 
the  fields  of  the  isolated  cavity. 

First,  we  write  the  CROW  fields  in  terms  of  the  isolated  cavity  fields  E^(x,y,z)  and 
H^{x,y,z): 

Ek  T)  =  -£'0  S  y,  z) 

f1=—oo 

H,  {X,  y)  =  H,f^  (x  -  nD,  y,  z)  (7-7) 

n=— 00 

where  and  are  normalization  constants  [1].  Taking  into  account  that 


=2n'^5{x-27tl) 


(7-8) 


It  follows  that 


EkMx’^ly) 


^  =\eS 


(2^'\ 

\  D  y 


2nn 


S  X  ^1  - —  \KMx’<ly) 


WkMx^qy) 


=  i/n 


^  27t^ 
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E  n=- 


<lx-k- 


2nn 
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HaMx^^yi\  ('7-9) 
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where  the  sum  over  m  ean  be  interpreted  as  a  sum  over  unit  cells.  This  relation  is  schematically 
depicted  in  Figure  7-2(a).  The  tight-binding  approximation  predicts  that  for  fixed  qy,  the  power 
spectrum  of  each  field  component  of  the  CROW  mode  is  given  by  a  set  of  delta-function  spikes 
weighted  by  the  power  spectrum  of  the  isolated  cavity.  The  spikes  are  shifted  by  k,  and  the 
separation  between  neighboring  spikes  isItt  /  D  . 


Power  spectra 


(a) 


I  \ 


Figure  7-2  (a)  Schematic  picture  of  tight-binding  prediction  for  the  power  spectrum  of  a 
CROW  mode,  (b)  CROW  power  spectrum  giving  a  large  value  of  Q{k)  compared  to  the 
isolated  cavity,  (c)  CROW  Power  spectrum  giving  a  small  value  of  Q{k)  compared  to  the 

isolated  cavity. 


Within  the  tight-binding  approximation,  {P(k)'j  can  be  found  by  substituting  Equation  (7-9) 
into  Equation  (7-5).  (U (k))  can  be  found  by  substituting  Equation  (7-7)  into  Equation  (7-3): 

{U(k))=  ^lljdV{e,,„(x,y,z)Ei  y  j 

^  y  m=—°° 

+  X  Z  (x  -  nD,  y,  z)H^  (x  -  mD,  y,  z)| 

J  (7-10) 
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The  m  =  n  terms  will  dominate,  leading  to  the  approximation  that 


1  OO  r 

((/)  —  —  ^  E^{x,y,z) 


+  l^oHl  H^(x,y,z) 


j’ 


(7-11) 


independent  of  k.  For  weak  eoupling  between  defeets,  we  substitute  £n,  the  dieleetrie  funetion  of 
the  isolated  eavity,  for  scrow,  and 


i  T,  z)El  \e^(x,  y,  z) 

^  m=~oo  y 

-liv. 


+  ^qHI  H^{x,y,z) 


n 

j 


(7-12) 


where  Uq,  is  the  field  energy  of  the  isolated  cavity. 


Using  the  expressions  for  (/^(k))  and  {U(k)'j ,  the  final  result  is  that 


Q(k)^co^{U, 


1  1^0  ff  dq,  dq^ 


2Ve„ 


ff  y  g 
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(7-13) 


Here,  we  have  assumed  that  the  variation  in  (o  across  the  CROW  band.  Am,  is  small  compared  to 
the  isolated  cavity  frequency  mo.  In  this  case,  (o{k)  can  be  approximated  by  mo  in  Equation  (7- 
13).  For  the  numerical  examples  presented  below.  Am/mo  is  on  the  order  of  10'  . 

Equation  (7-13)  implies  that  within  the  tight-binding  approximation,  increasing  Qo  of  the 
isolated  cavity  will  tend  to  increase  the  Q{k)  of  the  corresponding  waveguide.  This  is  because  an 
increase  in  Qo  corresponds  to  a  decrease  in  the  power  spectra  values  of  the  isolated  cavity  within 
the  light  cone.  However,  due  to  the  discrete  nature  of  the  sum,  Q(k)  can  have  a  strong 
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dependence  on  the  wave  vector.  This  is  illustrated  schematically  in  Figures  7-2  (b)  and  (c).  For  a 
particular  value  of  k,  it  may  happen  that  the  spikes  in  the  CROW  power  spectrum  line  up  with 
minima  (Figure  7-2  (b))  or  maxima  (Figure  7-2  (c))  of  the  isolated  cavity  power  spectrum.  The 
Qo  of  the  isolated  cavity  is  inversely  proportional  to  the  area  under  the  dashed  line,  while  Q{k)  of 
the  CROW  mode  is  inversely  proportional  to  the  red  or  blue  dashed  areas  shown.  It  is  apparent 
that  the  dashed  area  can  either  be  smaller  (Figure  7-2  (b))  or  larger  (Figure  7-2  (c))  than  the  area 
under  the  curve.  As  a  result,  Q{k)  can  be  either  larger  or  smaller  than  Qq. 

In  general,  the  validity  of  the  TBA  improves  as  the  separation  between  cavities  increases.  For 
decreasing  separation,  as  the  fields  of  neighboring  cavities  interact  more  strongly,  it  becomes 
less  accurate  to  write  the  CROW  fields  in  terms  of  the  isolated  cavity  mode  (Equation  (7-7)).  We 
note  that  in  the  numerical  examples  below,  we  consider  cavities  in  which  the  modes  are  well 
isolated  and  the  spacing  between  modes  is  large  compared  to  the  width  of  the  CROW  band.  In 
some  crow’s,  the  simple  model  of  a  well  isolated  mode  is  not  sufficient,  and  coupling  with 
other  modes  of  the  isolated  cavity  must  be  considered.  In  this  case.  Equation  (7-7)  would  be 
modified  to  include  a  sum  over  multiple  cavity  modes.  The  expansion  coefficients  in  the  sum  can 
be  found  as  in  Reference  [27],  and  Equations  (7-10)  and  (7-13)  would  be  similarly  modified.  In 
cases  where  nearest  neighbour  or  next-nearest  neighbor  coupling  are  significant,  the  m-n^O 
terms  should  be  included  in  the  calculation  of  (U (k))  (Equation  (7-10)). 


7.3  Numerical  Validation  of  the  Tight  Binding  Approximation 

To  test  the  TBA  model  of  the  CROW  quality  factor,  we  perform  three-dimensional  finite- 
difference  time-domain  (EDTD)  simulations  [28]  of  Q{k)  for  CROW  structures  and  compare  the 
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results  with  the  TBA  prediction  of  Equation  (7-13).  We  start  with  CROW  structures  based  on  the 
L2  cavity,  such  as  the  one  shown  above  in  Figure  7-1  (a). 


Real  part  of  Ey  in  L2 
cavity 


/  \  Ey  power  spectrum 
VCj  for  L2S3  with /f=0 


Ey  power  spectrum 
for  L2  cavity 


1 


■0.5 

^0 


Q^yfalln 


Figure  7-3  (a)  Mode  distribution  of  the  real  part  of  Ey  on  the  x-y  plane  0.2a  above  the 
surface  of  the  isolated  F2  cavity,  (b)  Power  spectrum  of  F2  cavity.  The  blue  circle  is  the 
light  cone,  (c)  Power  spectrum  of  CROW  mode  in  F2S3  structure  with  k=  0.  (d)  Visual 
test  of  tight-binding  approximation  via  comparison  of  CROW  power  spectra  to  that  of 
isolated  F2  cavity.  The  blue  dashed  lines  indicate  light  lines. 
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For  the  isolated  cavity  simulation,  a  computational  cell  of  35axl0^/3ax4a  was  used. 
Perfectly-matched  layer  (PML)  boundaries  were  used  in  all  directions,  and  the  grid  resolution 
was  20  grid  points  per  a.  In  Figure  7-3  (a),  we  plot  Ey  for  the  fundamental  mode  of  the  L2  cavity 
on  a  plane  0.2a  above  the  slab  surface.  The  mode  is  antisymmetric  with  respect  to  both  the  x-z 

plane  andy-z  planes.  Figure  7-3  (b)  shows  the  Ey  power  spectrum,  \E^^y{q^,qy)\  for  the  isolated 

cavity.  We  note  that  the  power  spectrum  is  zero  along  =  0,  due  to  the  fact  that  E^^  is  odd 

with  respect  to  x  =  0 .  In  Figures  7-3  (a)  and  (b),  the  field  magnitude  and  power  spectrum  are 
each  normalized  to  their  maximum  values. 

We  next  calculate  the  full  CROW  fields  and  corresponding  power  spectrum.  We  use  the 
FDTD  method  to  simulate  a  single  unit  cell  of  the  L2S3  CROW.  The  boundary  conditions  are 
fixed  to  be  Bloch-periodic  in  the  x-direction  with  fixed  wave  vector  k.  In  the  y-direction,  the 

length  of  the  computational  cell  is  chosen  to  be  large  enough  (10V3a)  for  Q{E)  to  be 
independent  of  length.  Perfectly-matched  layer  boundaries  are  applied  in  the  y  and  z  directions. 
The  grid  resolution  was  20  points  per  a.  We  divide  the  calculated  field  components  by  6*"^  to 
obtain  the  periodic  field  envelopes.  We  then  Fourier  transform  the  envelopes  and  shift  them  by  k 

in  the  ^x-direction,  a  procedure  that  can  easily  be  shown  to  yield  \Ej^  and  i{q^,qy)^  ■ 

In  Figure  7-3  (c),  we  plot  the  Ey  power  spectrum  of  the  L2S3  CROW  mode  for  A:  =  0.  The 
CROW  power  spectrum  is  only  nonzero  at  discrete  values  of  qx  separated  by  a  distance 
A^x  “  (1/  5)(2;r  !  a)  =  Q.2{2n  I  a) .  In  this  case,  only  three  values  of  qx  fall  within  the  light  cone, 

^^=0  and  =±0.2(2;r/a) .  The  power  spectrum  is  zero  for  q^=Q,  since  E^  is  odd  with 


respect  to  x  =  0 . 
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Within  the  TBA,  the  CROW  power  speetrum  ean  be  found  by  using  the  Dirac  delta  function 
to  sample  the  power  spectrum  of  the  isolated  cavity  mode  (Equation  (7-9)).  To  test  this 
approximation  quantitatively,  we  plot  the  values  of  the  L2  isolated  cavity  power  spectrum  (red 
curve)  and  the  L2S3  CROW  power  spectra  (black  dots)  along  =  0  in  Figure  7-3  (d).  Multiple 

values  of  k  are  shown  on  the  same  plot,  and  the  blue  dashed  lines  indicate  the  position  of  the 
light  lines.  The  overall  agreement  is  quite  good. 

As  a  further  check  of  the  TBA  model,  we  quantitatively  compare  the  results  for  Q{k)  to  those 
obtained  by  direct  simulation  of  the  full  CROW  structure.  Figure  7-4  shows  results  for  CROW 
waveguides  with  varying  separation  between  cavities:  Figures  7-4  (a-c)  correspond  to  F2S2, 
F2S3,  and  F2S4  CROWs,  respectively.  Qo  for  the  isolated  cavity  is  shown  by  the  red  lines  and 
has  a  value  of  1440. 

The  results  labeled  “FDTD”  and  “FDTD-FC”  represent  two  different,  direct  methods  for 
calculating  Q{k)  from  the  full  CROW  structure.  In  both  cases,  a  Bloch  periodic  unit  cell  with 
fixed  k  was  used,  as  described  above.  A  narrowband,  pulsed  source  was  used  to  excite  the  mode 
of  interest.  In  the  first  case  (“FDTD,”  black  lines)  Q{k)  was  calculated  from  the  field  decay  time 
after  source  turn-off  In  the  second  case  (“FDTD-FC”,  magenta  lines),  Q{k)  was  instead 
calculated  from  an  integral  over  the  CROW  fields  via  the  light-cone  approach  (Equations  (7-2)- 
(7-5)).  The  agreement  between  the  two  direct  methods  is  good,  indicating  the  validity  of  the 
light-cone  formulation  for  CROW  modes.  For  high-Q  modes,  one  might  expect  the  “FDTD” 
method  to  be  more  accurate  than  the  “FDTD-FC”  method,  which  depends  on  an  integral  over 
small  field  components  and  may  therefore  be  more  susceptible  to  numerical  error  than  the 
extraction  of  a  decay  time. 
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The  results  labeled  “TBA”  are  obtained  by  using  the  tight-binding  approximation  to  ealeulate 
Q{k)  from  the  isolated  cavity  mode,  via  Equation  (7-13)  (blue  lines).  The  TBA  is  shown  to 
predict  the  Q{k)  well.  The  maximum  difference  between  the  TBA  method  and  the  FDTD  method 
is  12%  for  L2S2,  7%  for  L2S3,  and  3%  for  L2S4.  As  expected,  the  accuracy  of  the  TBA 
improves  with  increasing  cavity  separation.  Importantly,  using  the  TBA  allows  the  prediction  of 
the  entire  Q{k)  curve  for  all  three  CROWs  from  only  a  single  calculation  of  the  isolated  cavity 
mode. 
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Figure  7-4  Q{k)  for  F2  CROWs  with  varying  cavity  separations:  (a)  F2S2,  (b)  F2S3,  and 
(c)  F2S4.  Dots  represent  calculated  values.  Fines  represent  guides  for  the  eye. 


From  Figure  7-4,  we  notice  that  the  A:-dependence  of  the  CROW  quality  factor,  Q{k),  is  not 
universal.  Q  may  increase  or  decrease  with  k,  depending  on  cavity  separation.  For  all  three  cases 
shown  (F2S2,  F2S3,  and  F2S4),  the  minimum  and  maximum  values  of  Q{k)  occur  at  either 
k-Q  or  k  =  kg^=  Q.5{2n / D) .  Depending  on  k,  the  CROW  mode  can  have  either  a  higher  or 

lower  quality  factor  than  the  isolated  cavity.  In  the  middle  of  the  Brillouin  zone,  the  quality 
factors  are  similar.  At  A:  =  0.25(2;r  /  D),  the  difference  between  Q{k)  of  the  CROW  and  Qo  of  the 
isolated  cavity  is  below  10%  for  all  three  cavity  separations  shown. 
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Interestingly,  while  the  L2  eavity  mode  eannot  emit  light  vertieally  due  to  symmetry 
eonsiderations  [29],  the  CROW  ean.  For  the  isolated  eavity,  symmetry  prohibits  eoupling  to  a 
plane  wave  propagating  in  the  z-direetion.  However,  leakage  oeeurs  in  off-normal  directions, 
yielding  a  finite  quality  factor.  For  the  CROW  mode  at  Q{k  =  0) ,  symmetry  again  prohibits 
coupling  to  a  z-propagating  plane  wave.  However,  translational  symmetry  requires  that  any 
leakage  be  in  the  z-direction.  The  finite  value  of  Q{k  =  0)  indicates  that  vertical  leakage  does 
occur;  in  this  case,  it  arises  from  interference  between  off-normal  plane  waves  radiated  from  the 
individual  cavities. 

7.4  Factors  Influencing  the  CROW  Quality  Factor 

The  results  above  verity  that  the  TBA  model  is  a  good  predictor  of  Q{k)  in  CROWs  based  on 
the  L2  cavity.  Using  the  TBA,  we  can  gain  insight  into  how  the  CROW  quality  factor  depends 
on  wave  vector,  cavity  separation,  and  the  quality  factor  of  the  isolated  cavity. 

In  Figure  7-5,  we  plot  Q{k)  for  L2  CROWs  as  a  function  of  cavity  separation  (black  lines). 
Triangles  indicate  Q{0),  and  circles  indicate  Qiksz)-  All  values  are  calculated  using  the  TBA.  We 
see  from  the  figure  that  Q{0)  and  Qikgz)  oscillate  around  the  Qo  of  the  isolated  cavity,  shown  by 
the  dashed  line.  Interestingly,  Q{0)  and  Qiksz)  are  always  on  opposite  sides  of  Qo.  As  a  result,  it 
is  never  the  case  that  Q{k)  is  larger  or  smaller  than  Qo  for  all  k.  As  the  number  of  separation 
holes  is  increased,  both  Q{0)  and  Qiksz)  approach  Qo.  This  trend  is  intuitive;  in  the  limit  of  large 
separation,  the  CROW  can  be  viewed  as  a  chain  of  isolated  cavities. 

We  next  increase  the  length  of  the  cavity.  Results  for  L3  CROWs  are  shown  by  red  lines  in 

Figure  7-5  (a).  It  can  be  seen  from  the  figure  that  Q{0)  and  Qikgz)  again  oscillate  around  the 
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isolated  cavity  value.  Moreover,  the  increase  in  the  isolated  cavity  Qo  between  L2  and  L3 
corresponds  to  an  increase  in  ^(A:).  Similarly,  we  can  form  CROWs  from  L4  or  L5  cavities  (blue 
and  magenta  lines  in  Figure  7-5  (a)).  Within  the  family  of  L2  to  L5  cavities,  the  CROW  quality 
factor  increases  with  the  quality  factor  of  the  isolated  cavity.  This  result  corresponds  to  the 
intuitive  argument  presented  above  in  (Section  7.2).  As  the  quality  factor  of  the  isolated  cavity 
increases,  the  components  of  its  power  spectrum  inside  the  light  line  decrease.  This  also  tends  to 
decrease  the  leaky  components  of  CROW  power  spectrum,  which  may  be  obtained  from  the 
cavity  power  spectrum  via  Fourier  component  sampling. 

For  each  cavity  length  (or  data  color)  shown  in  Figure  7-5  (a),  we  notice  that  there  are  nodes 
where  Q{0)  and  Qiksz)  cross.  At  these  points,  we  expect  Q{k)  to  be  flat.  In  Figure  7-5  (b),  we  plot 
Q{k)  for  CROWs  designed  with  a  number  of  separation  holes  equal  to  the  first  node  in  Figure  7-5 
(a).  In  all  cases,  the  variation  in  Q{k)  is  small  compared  to  Figure  7-5  (a).  For  L3  and  L4,  the 
number  of  separation  holes  is  equal  to  3.  We  thus  obtain  designs  for  CROWs  whose  loss  is 
similar  to  the  isolated  cavity  across  the  entire  band.  For  L2  and  L5,  the  ideal  separation  is  non¬ 
integer.  While  loss  values  can  be  calculated  within  the  TBA,  it  is  not  obvious  how  to  achieve 
non-integer  separations  in  a  real  CROW  structure.  However,  we  have  verified  via  direct  FDTD 
simulations  that  for  the  L2  cavity,  simply  using  a  CROW  unit  cell  of  4.4a  and  allowing  a  slight 
overlap  of  circular  holes  at  each  edge  of  the  cell  yields  a  similarly  flat  Q{k)  trend  to  that  pictured, 
though  at  somewhat  higher  values.  We  have  also  verified  via  FDTD  simulations  (not  shown)  that 
in  L3,  L4  and  L5  cavities,  the  TBA  method  gives  accurate  results  for  Q{k)  provided  that  the 
number  of  separation  holes  is  greater  than  2. 
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Figure  7-5  (a)  Quality  factor  for  CROWs  based  on  different  constituent  cavities.  Triangles 
indicate  Q(0);  circles  indicate  Qiksz)-  Dashed  lines  show  Q  for  single  cavity.  Calculations 
are  performed  using  the  TBA  method,  (b)  Q{k)  for  CROWs  designed  with  a  number  of 
separation  holes  equal  to  the  first  node  in  (a). 


7.5  Summary 

In  conclusion,  we  have  applied  the  tight-binding  approximation  (TBA)  to  predict  the 

radiation  loss  of  CROW  waveguides  in  photonic  crystal  slabs.  Our  formulation  shows  that  as  the 

quality  factor  of  the  isolated  cavity  increases,  the  CROW  quality  factor  tends  to  increase  as  well. 

This  trend  can  be  understood  within  the  light-cone  picture  of  radiation  loss.  For  both  isolated 

cavities  and  CROWs,  the  quality  factor  is  related  to  the  spatial  power  spectrum  of  the  fields 

above  the  photonic  crystal  slab.  As  the  magnitude  of  the  components  inside  the  light  cone 

decrease,  the  quality  factor  increases.  Within  the  TBA,  the  CROW  power  spectrum  can  be 

approximated  by  sampling  the  isolated  cavity  power  spectrum  in  Fourier  space.  The  loss  of  the 

CROW  thus  tends  to  follow  that  of  the  isolated  cavity.  We  have  verified  the  validity  of  the  TBA 

by  comparing  predicted  values  of  Q{k)  for  L2  CROWs  with  those  obtained  by  direct  simulation. 
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Importantly,  our  method  allows  the  calculation  of  the  wave-vector  dependent  CROW  quality 
factor  for  a  range  of  cavity  separations  from  a  single  simulation  of  the  isolated  cavity  fields.  As  a 
result,  it  is  possible  to  scan  various  CROW  configurations  with  minimal  computational  power.  In 
particular,  we  demonstrate  that  it  is  possible  to  design  waveguides  with  low,  fiat  losses  across  the 
entire  Brillouin  zone,  even  at  small  intercavity  separations. 

We  expect  our  results  to  be  useful  for  a  range  of  applications.  For  slow-  and  stopped-light 
applications,  the  fidelity  of  signal  transmission  at  the  end  of  a  delay  line  depends  on  achieving 
low,  fiat  loss  across  the  signal  bandwidth.  In  coupled-cavity  lasers,  the  ability  to  tailor  the  Q{k) 
loss  profile  could  provide  a  useful  method  to  control  the  effects  of  mode  competition  and  tune 
the  directional  output  of  the  laser.  Lastly,  in  the  context  of  light-assisted  self  assembly  above 
photonic-crystal  slabs  [24,  25],  our  method  may  facilitate  the  design  of  large-area,  coupled- 
cavity  modes  with  high  quality  factor  that  nevertheless  couple  to  normally-incident  radiation. 
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Appendix: 

Force  Calculation  in  Time  Domain 

The  Maxwell  stress  tensor  (MST)  is  a  general  formalism  whieh  allows  to  ealeulate  the  foree 
exerted  on  an  object  for  a  given  electromagnetic  field  distribution.  To  evaluate  optical  force 
exerting  on  a  component,  the  MST  is  numerically  integrated  over  a  closed  surface  surrounding 
the  object: 

where  a  and  (5  indicate  direction  x,  y  or  z,  n  is  the  outward  normal  to  the  closed  surface  S,  and 
the  stress  tensor, 


-  E  +  /iH  ■  H) . 

To  calculate  the  stress  tensor  Tap,  full  electromagnetic  field  distributions  are  needed,  which 
can  be  obtained  from  an  appropriate  full-vectorial  electromagnetic  solver.  We  show  an  example 
of  using  a  FDTD  solver  and  the  MST  method  to  solve  for  optical  force  on  polystyrene  particle  by 
the  enhanced  light  of  guided  resonance  mode  in  Slot-Suzuki-phase  lattice.  An  artificial  surface  is 
taken  to  enclose  the  polystyrene  nanoparticle.  As  shown  by  Figure  A-1,  here  an  integration  cube 
is  taken  with  20  nm  away  from  the  particle  (refractive  index-'1.60).  It  is  highly  recommended  to 
check  if  the  calculated  force  was  significantly  varied  with  different  integration  surface.  A  force 
dependence  on  surface  indicates  questionable  modeling. 
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Figure  A-1  Schematic  of  one  unit  cell  of  Slot-Suzuki-phase  lattice.  A  cube  is  taken 
enclosing  the  particle  as  marked  by  blue  dashed  lines  for  integrating  the  Maxwell  Stress 

tensor. 

Optical  forces  are  correlated  with  optical  performance.  As  a  starting  point,  the  device  without 
particles  are  simulated  to  get  an  idea  of  critical  optical  performance  including  resonance 
wavelength,  quality  factor  and  mode  profde.  A  test  of  mesh  size  should  be  conducted  to  check 
how  parameters  converge.  Generally  speaking,  finer  meshes  are  favorable  for  more  accurate 
calculation  in  the  cost  of  computational  resource.  By  convergence  test,  a  reasonable  mesh  size 
could  be  decided  to  achieve  acceptable  results  without  occupying  too  much  resource. 

For  force  calculation,  a  particle  is  included  in  the  computational  cell.  Finer  meshes  are 
required  in  the  region  of  the  particle  and  small  features  like  slots  in  the  device.  In  this  case,  it  is 
advantageous  to  use  a  FDTD  package  which  can  define  un-uniform  meshes.  One  FDTD 
implementation,  Lumerical,  is  applied  in  our  modeling  for  force  calculation  of  particles  with 
radius  of  25  nm. 

Light  distribution  varies  with  wavelength,  and  thus,  optical  force  is  dependent  on  wavelength. 
Using  frequency-domain  solver,  a  group  of  calculations  are  needed  to  scan  over  the  wavelength 
of  interest.  By  using  time-domain  solver  like  FDTD,  we  can  obtain  the  entire  frequency  spectrum 
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of  force  in  a  single  simulation,  by  Fourier- transforming  the  response  to  a  short  pulse.  The 
frequeney  resolution  is  deeided  by  the  total  running  time.  For  narrow-bandwidth  features  like 
high-Q  resonanee,  FDTD  simulation  needs  to  be  run  by  long  time  to  aeeurately  compute  the 
fields  and  forees.  If  the  resonanee  is  exeited  by  aehromatie  light  souree  with  different 
frequencies,  foree  at  eaeh  frequeney  eontributes  to  the  effeetive  foree. 

For  high-Q  resonanee  mode,  it  is  eonvenient  to  perform  the  foree  ealeulation  in  the  time 
domain  rather  than  the  frequeney  domain.  In  FDTD,  we  excite  the  mode  using  a  dipole  souree 
inside  the  slab.  For  effieient  exeitation,  the  bandwidth  of  the  dipole  souree  is  larger  than  the 
resonanee  linewidth.  A  eouple  of  dipole  sourees  may  be  positioned  at  positions  with  expected 
strongest  intensity.  To  have  an  idea  of  the  mode  profile,  an  eigenmode  solver,  e.g.  MPB  or  GME, 
ean  be  used  before  FDTD  ealeulation.  After  the  souree  is  tuned  off,  the  fields  in  the  resonator 
start  to  deeay.  The  deeaying  rate  ean  be  deseribed  by  quality  faetor  Q.  The  instantaneous 
eleetromagnetie  fields  on  the  integration  cube  are  reeorded  for  several  optieal  periods.  Equation 
(A-1)  is  utilized  to  ealeulate  the  time-dependent  foree.  Typieally,  the  time  seale  for  meehanieal 
response  is  longer  than  the  optieal  period,  and  the  time-averaged  foree  over  the  optieal  periods 
may  be  eonsidered.  No  long  running  time  is  needed.  Another  advantage  of  time-domain  method 
is  that  all  frequeney  eomponents  are  taken  into  aeeount  for  foree  ealeulation. 

The  foree  ean  be  normalized  over  eleetromagnetie  field  energy  U.  After  the  souree  is  tuned 
off,  the  instantaneous  eleetromagnetie  fields  in  the  unit  cell  are  recorded  for  several  optieal 
periods.  The  time-dependent  field  energy  is  ealeulated  by  taking  the  volume  integration  over  the 
unit  eell: 

^  y’  \H{x,  y,  zf  )  (A-2) 

^  V 
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The  time-averaged  energy  over  the  optieal  periods  can  be  calculated,  and  the  corresponding  F/U 
is  obtained. 

The  force  can  also  be  normalized  by  emission  power  P=coUIQ. 


F_F  Q 
P~  U  CO 


(A-2) 


In  steady  state,  the  emission  power  is  equal  to  incident  power.  The  FIP  can  be  calculated  by 
the  value  of  FlU  and  Equation  (A-3). 


Manipulating  Nano  wires  by  Optical  Forces 

In  the  thesis,  different  configurations  have  been  discussed  to  achieve  optical  forces  to  either 
reposition  waveguides,  resonators  and  nanoparticles.  It  will  be  also  interesting  to  manipulate 
nanowire  arrays  and  tune  optical  performance  by  optomechanical  effects. 

Nano  wire  arrays,  as  shown  by  Figure  A-2,  is  a  2D  grating  structure.  It  is  known  that  when 
light  is  normally  incident  upon  a  ID-  or  2D-periodically  patterned  structure,  it  can  excite  guided 
resonance  modes  with  increased  electromagnetic  field  intensity  within  the  grating.  It  is  expected 
that  that  the  enhanced  field  intensity  between  two  closely-positioned  nanowires  will  give  rise  to 
an  optical  force  on  the  nano  wires. 
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Figure  A-2  Schematic  of  square-lattice  nanowire  arrays  standing  on  substrate.  In 
each  unit,  there  are  two  closely-positioned  nano  wires.  Assume  the  material  is  InP. 

Light  is  normally  incident  into  the  array  from  free  space. 

An  example  design  is  shown  in  Figure  A-3  (a).  For  simplicity,  one  unit  cell  is  plotted.  An  air 

-5 

gap  of  68  nm  is  between  two  nanowires  with  dimensions  of  135x135x2700  nm  .  Assume  the 
material  is  InP.  The  transmission  spectrum  (black  curve  in  Figure  A-3  (b))  is  calculated  by  there- 
dimensional  FDTD  method.  Guided  resonance  modes  are  identified  by  transmission  dips 
imposed  on  Fabry-perot  fringes.  The  optical  force  as  function  of  wavelength  is  plotted  by  the 
red  curve  in  Figure  A-3  (b).  The  force  is  calculated  by  Maxwell  Stress  Tensor  method.  It  is 
notable  that  force  magnitude  is  significantly  enhanced  on  guided  resonance,  especially  for  the 
resonance  mode  at  1550  nm.  The  negative  sign  means  an  attractive  force.  Figure  A-3  (c)  shows  a 
side  view  of  the  ^x-component  distribution  at  wavelength  1550  nm.  Note  that  the  field  intensity 
is  enhanced  in  the  air  gap.  It  is  expected  the  field  confinement  will  generate  large  optical  force 
on  both  nanowires. 
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Figure  A-3  (a)  Schematic  of  one  unit  cell.  Dimensions  are  as  given.  The  air  gap 
between  two  nanowires  is  68  nm.  (b)  Transmission  (black  curve)  for  normal 
incidence.  Guided  resonance  modes  are  identified  by  transmission  dips  imposed 
on  Fabry-perot  fringes.  Force  spectrum  (red  curve)  in  the  same  wavelength 
window,  (c)  Side  view.  The  blue  and  red  color  are  for  the  component. 


If  1550-nm  laser  light  is  normally  incident  with  power  of  1  mW  per  unit  cell,  the  optical 
force  is  around  -2  nN  on  each  nanowire,  causing  a  maximum  bending  of  2  nm  at  the  nanowire 
top.  The  transmission  spectrum  can  be  modified  by  adjusting  air  gap.  With  the  air  gap  reduced 
by  2  nm,  the  resonance  wavelength  will  be  red  shifted  by  around  0.3  nm,  one  third  of  the 
resonance  linewidth.  Larger  change  in  the  air  gap  is  desirable  for  easier  optical  probing  by 
reading  transmission  shift.  For  fixed  input  power,  larger  displacements  can  be  obtained  by 
modulating  the  incident  optical  power  at  the  mechanical  resonance  frequency. 
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(a)  (b) 


Figure  A-4  SEM  image  of  top-down  fabricated  silicon  (a)  square  lattice  of 
nanowires,  and  (b)  graphite  lattice  of  nanowires. 

Nanowires  can  be  fabricated  by  top-down  method.  An  array  of  600-nm-tall  silicon 
nanowires  have  been  fabricated  using  e-beam  lithography  and  inductively-coupled  plasma 
reactive  ion  etching  (ICP-RIE).  We  use  a  2-pm  thick  silicon  layer  on  top  of  a  2  pm  silica  layer 
(Ultrasil).  The  sample  was  spin-coated  at  3000  rpm  with  HSQ  (XR-1541),  a  negative  e-beam 
resist.  Using  RaithlSO  e-beam  system,  a  square-lattice  pattern  has  been  exposed  with  an 
acceleration  voltage  of  30  KV  and  dose  of  2000  uA»s/cm  .  The  sample  is  developed  in  NaOH  (1 
wt%)  /  NaCl  (4  wt%).  The  pattern  is  transferred  from  the  resist  to  the  silicon  layer  by  an  ICP- 
RIE  machine  using  a  gas  mixture  SE6/C4Eg.  HSQ  is  very  sensitive  to  temperature  and  moisture.  It 
is  highly  recommended  to  prepare,  exposure  and  develop  quickly.  HSQ  requires  longer  e-beam 
exposure  time  than  PMMA.  Under  our  etching  conditions,  the  etching  selectivity  of  HSQ  over  Si 
is  1:8.5.  Eigure  A-4  (a)  and  (b)  respectively  show  a  square  lattice  and  graphite  lattice  of  nanowire 
arrays. 
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